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I. INTRODUCTION 


Flight simulation has been an important computer graphics application, 
embracing a range of systems from a $32.00 program for a personal computer 
(Ref. 1] to special purpose machines costing millions of dollars [Ref. 2]. The 
capabilities of these systems are spread across a range nearly as wide as their 
costs, with great variances in speed (frames displayed per second). realism, 
flexibility, and area of flight. We present here a system that is relatively 
inexpensive, vet still fast enough to present a real-time three-dimensional view of 
digitized terrain. We built this system on a commercially available, high- 
performance graphics workstation, the Silicon Graphics, Incorporated IRIS-2400 
Turbo. The IRIS system was selected because of its local availability and its 
performance capabilities. The flight simulator presented here does not use the 
natural color and shape of individual terrain elements (in order to achieve real- 
time performance), but it is sufficiently realistic to provide the feeling of flight 


and allow identification of the displayed terrain and targets. 


A. FOG-M 
1. Background 
The project presented here was built in response to the United States 
Army Combat Developments Experimentation Center’s need to simulate the 
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operation of the Fiber-Optically Guided Missile (FOG-M) [Ref. 3], but this missile 
is also being considered for use by the United States Marine Corps [Ref. 4]. 
Simulation is necessary in order to test and evaluate the tactics, doctrine and 
training requirements associated with the missile without the expense and danger 
of actual firings during simulated combat field trials. The FOG-M is a generic 
family of remotely-piloted, video-guided munitions, but for the purpose of this 
prototype simulator, the weapons are all logically equivalent, and the entire 


family is referred to as “‘the missile.”’ 


In order to avoid security constraints, the 
parameters and operational characteristics used in this work were not taken from 
exact FOG-M specifications. The parameters and technical specifications are all 
estimates, based on reasonableness and consistency with general, unclassified 
descriptions of the FOG-M. 
2. Description 

The actual FOG-M missile is six inches in diameter, five and one-half feet 
high, weighs eighty-three pounds, and costs about $20,000 [Ref. 4]. It has a video 
camera mounted in its nose, which transmits a black-and-white picture back to 
the operator’s console (which consists of a television screen, a computer, and a 
joystick) over the fiber-optic link. (The simulator display offers the user the choice 
of either color or black-and-white; color is the default for the simulator despite the 
operator view of the missile being black-and-white. The color compensates for 


some of the loss in realism and identifiability inherent in a _ polygonal 


representation of natural objects). Before launch, in normal operation, the missile 
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is given a general direction to a target and the altitude of the highest point within 
its range. The simulator allows values in excess of FOG-M operational 
capabilities for speed, range, and altitude above ground level (AGL), but the 
default values of two hundred knots, ten kilometers, and one thousand meters are 
characteristic of this type of missile. As soon as ie missile is in position, it begins 
transmitting video images. When launched, the missile rises to approximately 
two hundred feet above the highest terrain point, and then levels off in horizontal 
flight in the targeted direction. The operator controls the pan and tilt angle of 
the camera with the joystick, and can dial in changes to the heading and altitude 
of the missile. The operator also has the capability to zoom the camera’s field of 
view from eight degrees to fifty-five degrees, and to designate (‘‘lock-on”’ to) a 


target for automatic homing by the missile. 


B. ASPECTS OF FLIGHT SIMULATION 

There are many aspects to flight simulation. Modern commercial simulators 
provide sophisticated mock-ups of cockpits and controls and highly detailed out 
the window views. By mounting the simulator on a moving platform, a true sense 
of the physical feelings of acceleration and roll can be achieved. These systems 
also cost millions of dollars. 

One of the first decisions that must be made when designing a flight simulator 
is, ‘‘For what purpose will the simulator be used?’’ The answer to this question 


drives most of the design decisions that have to be made. Since the purpose of 
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this project is to provide a simulation of the FOG-M missile as viewed from its 
operator’s console, it is felt that the most important items to model are the 
simulated video display of the terrain and the actual operator controls. The 
terrain should appear realistic enough that its major features are recognizable to a 
person familiar with the area. The controls should allow for the same 
functionality as the actual console. The simulator must, of course, also provide a 
feeling that the missile is in motion over the terrain. The effectiveness of the 
feeling of motion provided by a flight simulator can be largely measured by two 
criteria: the realism of the displayed scene and the update rate of the display. 
1. Realism 

Many factors contribute to the perceived realism of a displayed natural 
scene. Early attempts to quantitatively measure realism consisted of counting the 
number of ‘‘edges’’ or lines that a scene contained. This later gave way to 
counting the number of “‘faces’’ or polygons in a scene. Since each polygon was 
colored in a single shade, it was felt that each polygon represented a single ‘‘bit”’ 
of information in the scene. Therefore, the more polygons the scene contained, 
the more “‘realistic’’ it was felt to be [Ref. 5:pp. 27-28]. 

The latest advances in computer graphics have also made this measure of 
effectiveness obsolete. With the introduction of systems that are able to fill 
polygons with textured patterns, a single polygon can now contain thousands of 
‘bits’? of information. As a result, a scene drawn with a few textured polygons 


can appear more realistic than one with an order of magnitude more untextured 


13 


ones. Early textures consisted of superimposing things such as mathematical 
noise functions or stripes on the polygons. More recent advances have allowed the 
texture to be derived from digital photographs of a similar scene. For example, 
polygons representing a part of terrain covering by meadow could be filled with a 
digital texture derived from an aerial photograph of a meadow [Ref. 5: p. 28]. 

Since most currently available graphics workstations do not support the 
use of texture filled polygons, the use of texture was deemed to be outside the 
scope of the current project. Rather, the project’s work concentrated on 
determining how realistically a scene could be rendered in real-time incorporating 
only the use of lighting and shading models along with terrain hidden-surface 
algorithms. These topics are covered in more detail in Chapter V. 

2. Frame Update Speed 

Another important aspect of a flight simulator’s performance is the speed 
at which it is capable of displaying successive frames in a scene. The faster the 
update rate, the more continuous the motion appears. As a reference, standard 
motion picture film is projected at a rate of twenty-four frames per second. 
Although the IRIS workstation is capable of displaying up to sixty frames per 
second, the amount of computation that must be done between successive frames 
in the simulation has limited the update rate to an average of three frames per 
second. While this presents a less than smooth motion, it is felt to be adequate 


for the purposes of the prototype. 
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C. ORGANIZATION 

The above sections of this chapter have provided background on flight 
simulation in general, and the missile whose flight is specifically being simulated. 
Chapter II provides an overview of the hardware used in running the simulation. 
The structure anal content of the Defense Mapping Agency (DMA) Digital 
Terrain Elevation Data (DTED) are discussed in Chapter III. Chapter IV 
discusses the motivation behind and creation of the two-dimensional contour map 
displays. Chapter V covers the storage and use of the DMA DTED to produce a 
lighted and shaded three-dimensional polygonal terrain display. The mathematics 
and process involved in simulating flight over the terrain are detailed in Chapter 
VI. Chapter VII discusses the creation, insertion, animation, and designation of 
targets. Chapter VIII covers the creation and drawing of cultural features. 
Chapter IX contains a user’s guide for operation of the FOG-M simulator. 
Chapter X concludes with a discussion of limitations, future extensions and 
research topics, and summarizes the research conducted. Narrative descriptions of 
the modules and listings of the program source code for each of the modules are 


included as Appendices A and B respectively. 
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Il. COMPUTER SYSTEM 


AS discussed in Chapter I, flight simulators are nothing new. The significance 
of this work lies in the speed with which it was developed, the display rate 
achieved, and the realism and fidelity of the display in comparison to the cost of 
the system that supports it. This project was technically feasible only because of 
the capabilities and high performance of the IRIS (Integrated Raster Imaging 
System) Turbo 2400 Graphics Workstation, manufactured by Silicon Graphics, 
Incorporated. Others have also used the IRIS as a base on which to build flight 
simulators [Ref. 6]. This low-cost ($50,000 to $100,00) three-dimensional display 


system is summarized in Figure 2.1 and is discussed more fully below. 


A. HARDWARE AND SYSTEM OVERVIEW 

The IRIS has a conventional Von Neumann type computer architecture but 
adds custom-built special purpose VLSI circuits and a pipelined design to provide 
the graphics functions that are implemented in software on other comparably- 
priced workstations. Conceptually, there three pipelined components in the IRIS 
hardware: the applications/graphics processor, the Geometry Pipeline. and the 
raster subsystem |Ref. 7:p. 1-1]. The applications/graphics processor is a 
conventional Motorola MC68020 processor running at 16.7 MHz. This processor 
runs the applications program(s) within a UNIX System V operating system. 


16 


ETHERNET to Vax and other IRIS 


af 


32 bit 16.7 MHz Motorola MC68020 CPU 





6 Megabytes of CPU Memory 

32 1024x768 bitplanes of Display Memory 

Hardware matrix multiplier & floating point accelerator 

Hardware Gouraud shading, depth cueing & backface polygon removal 
12 pipelined custom VLSI Geometry Engines 

16-bit Z-buffer for Hidden Surface Elimination 


2 72 Megabyte Winchester Disk Drives 
60 Hz non-interlaced 19" RGB high resolution monitor 


83 key up-down encoded keyboard 


3 button mouse 


32-button and 8-dial valuator boxes 
Unix System V 
Ethernet to VAX’s 


IRIS Graphics Library 


Features of the IRIS Turbo 2400 Graphics Workstation 
Figure 2.1 
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Applications either issue graphics commands in tmmediate mode, in which case 
they are sent through the Geometry Pipeline immediately as individual graphics 
primitives, or compile graphics commands into graphical objects, in which case 
they are sent through the Geometry Pipeline as a single geometric entity when 
explicitly called at some later point in time. 

The Geometry Pipeline takes commands in terms of the user’s world 
coordinates, performs specified matrix transformations on them using the matrix 
multiplier and floating point accelerator built into the hardware, and then clips 
and scales the transformed coordinates into screen coordinates. The raster 
subsystem takes the screen coordinates output by the Geometry Pipeline and 
updates the bitplanes (display memory) to contain the lines, polygons, or 
characters specified by the input coordinates. The raster subsystem also performs 
polygon fill, shading, depth-cueing, and hidden surface removal. A conventional 
video controller uses the values in the bitplanes and the color table to produce an 


image on the monitor. 


B. SOFTWARE 

The C programming language is native to UNIX and is the language used for 
all ot the IRIS system software. The IRIS Graphics Library. which provides both 
high- and low-level graphics and utility commands, can be called in C, 
FORTRAN, Pascal, or LISP. However, due to the built-in bias of UNIX and the 


IRIS, plus the local pool of knowledge, the C programming language was the 
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pro forma choice for programming all parts of the FOG-M simulator. The IRIS 


User’s Guide |Ref. 7] breaks the Graphics Library commands into the following 


twelve categories: 


Global State commands initialize the hardware and control global variables, 
and are used mostly in FOG-M’s init_irts routine. 


Drawing Primitives are used throughout FOG-M. They create points, lines, 
polygons, circles, arcs, and text strings. 


Coordinate Transformations specify mappings within and between user- 
defined world coordinates and screen coordinates. These are used to move 
targets and to simulate flight. 


Drawing Attribute commands specify textures and fonts. Although texture 
would greatly improve the appearance of the terrain, the IRIS provided 
textures are applied in the screen coordinate system, so they do not scale or 
tilt to conform to the terrain, and produce a very artificial result. 


Display Mode and Color commands determine how the bitplanes are used 
and what colors appear on the screen. These include the commands that set 
double-buffering, establish writemasks, and define the color table. 


Input / Output commands initialize and read the dials and mouse. 


Object Creation and Editing commands allow manipulation of complex 
displays as a single entity. They are used in all FOG-M displays. 


Picking and Selecting commands are not used in FOG-M. 
Geometry Pipeline Feedback commands are not used in FOG-M. 


Curve and Surface commands draw complex curves and smooth surfaces. 
Experiments with these produced more realistic terrain images, but not even 
close to real-time, making flight animation impossible. 


Shading and Depth—cueing commands provide Gouraud shading of polygons 
and intensities that vary with distance from the viewer. 


Teztport commands define an area of the screen for text. They are not used 
in FOG-M. 


Also available on the system, and used by FOG-M, are the math library with 


sine, 


cosine, arctangent, hypotenuse, and exponentiation functions, and routines 


that access the system clock in order to determine elapsed time. 
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II. DIGITAL ELEVATION TERRAIN DATA 


A. INTRODUCTION 

Unlike other flight simulation systems, which may rely on manual creation of 
the terrain [Ref. 8], the source data for the terrain in the FOG-M simulation is a 
Defense Mapping Agency (DMA) digital terrain elevation database (DTED) for 
Fort Hunter-Liggett. California. The database is not Level 1 DTED, but rather a 
DMA special product produced about 1980 at a higher resolution than normal 
Level 1 DTED [Ref. 9]. Level 1 DMA data contains elevation points spaced at 
three arc-second intervals, or approximately every one hundred meters. The Fort 
Hunter-Liggett special data contains points at twelve and one-half meter spacing, 


which is eight times the resolution of Level 1 data. 


B. COVERAGE 

The area covered by the database is thirty-six kilometers wide and thirty-five 
kilometers high, with 6400 data points per square kilometer. This area includes 
most of Fort Hunter-Liggett plus some surrounding land, and is bounded by 
latitudes 26° 05° 00°’ (to the north) and 35° 50’ 00°’ (south) and longitudes 
121° 04° 30°’ (east) and 121° 20° 30°’ (west). In terms of Universal Transverse 
Mercator (UTM) coordinates, the area has easting (X) of 10SFQ41000 to 
10SFQ77000 and northing (Y) of 10SFQ60000 to 10SFQ95000. The database 
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appears to be based on DMA forty foot interval contour map products, because 
peaks tend to have flattened tops. This was confirmed both by a comparison of - 
surveyed instrumentation sites on or near peaks with their digital terrain values 


[Ref. 10: pp. 1-2], and by a Bezier surface patch image of the data created locally. 


mo lL RUCTURE 

The data is stored in an unformatted sequential file that is organized as a 
stream of integers. Each integer (sixteen bits) represents both the vegetation code 
and bald terrain elevation in feet at one sampling point, as illustrated in Figure 


3.1 below. 


Veg. Code Bald Terrain Elevation 
Ditcselioel4 lout 10°98 7 6 «6 424 3 2 1 (iO 


Figure 3.1 DTED Data Encoding 


The thirteen low-order (rightmost) bits contain the elevation, allowing a range 
irom zero to 8191 teet. although the highest point in the database is 3744 feet. 
The three high-order (leftmost) bits specify one of eight vegetation codes, which 
are given in Table 3.1 below. Vegetation codes are only available for points 


within the boundaries of Fort Hunter-Liggett proper. The file is written one 
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TABLE 3.1 DTED VEGETATION CODES 


Description 


Less than one meter 
One to four meters 
Four to eight meters 

Eight to twelve meters 
Twelve to twenty meters 
Greater than twenty meters 
No data available 

Unused 





square kilometer at a time, beginning with the lower left one kilometer grid square 
(41,60), proceeding up the column to the upper left grid square (41,94), then 
doing the next column from bottom to top (42,60 to 42,94) and so on; the upper 
right one kilometer grid square (76,94) is the last one written. Within each one 
kilometer grid square, the individual data points are written in the same pattern, 
beginning with the lower left, doing each column from bottom to top, and doing 
the columns from left to right. This file layout is summarized in Figure 3.2. The 
position in the file of the elevation for a point expressed in five digit local UTM X 


and Y coordinates is found as shown in Equation 3.1. 


position = 35 * (integer(X/1000) — 41) + (integer( Y/1000) — 59) (3.1) 


DA EOCATION 
The compiete DTED fle occupies 16,128,000 bytes or storage. Due to a local 
shortage of available disk space, this file must permanently reside on the UNIX 


VAX 11/785 system rather than on the IRIS system. The FOG-M simulator 
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Figure 3.2 DTED File Layout 
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presently operates on a ten kilometer square extract from this database. A 
program on the VAX called make—database—e allows interactive specification 
of the area and resolution desired, and produces an extract. This extract is sent 
over the Ethernet to the IRIS to serve as the input for a FOG-M run. However, if 
the data is sent directly, it is received with each pair of bytes swapped, so another 
program, swapdima, is run on the VAX before transmittal. This program swaps 
the low- and high-order bytes of each integer so that the swapping during 


transmission is cancelled. 
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IV. TWO-DIMENSIONAL TERRAIN MAP PORTRAYAL 


The two-dimensional representation of the terrain was begun as the first 
graphics portion of the system, in order to gain familiarity with the IRIS graphics 
workstation and the Defense Mapping Agency (DMA) digital terrain elevation 
data (DTED). Contour maps are the traditional approach to two-dimensional 
terrain portrayal, and thus were the basis for the two-dimensional images of the 
terrain generated here (Figure 4.1). Although these two-dimensional images are 
not true contour maps, they are still referred to as such in hie study because of 
their close relation and common origin. The algorithms for determining and 
drawing the forty foot contour lines found on a normal contour map seemed non- 
trivial, so a simpler alternative was chosen. Hach elevation datum is represented 


by a tile, with the implicit X and Z (easting and northing, respectively) 


coordinates of the elevation datum being the center of the tile. 


A. COLORS 

The color of a tile is determined by its vegetation code, and its intensity (or 
shading) by its elevation. The intent was to use green for tiles with vegetation 
and brown for tiles without vegetation. However, the DTED vegetation codes 
lump together both “no vegetation”’ and ‘“‘vegetation less than one meter high.”’ 
Brown tiles thus include both unvegetated areas (e.g. rock slabs, areas above the 
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treeline) and grasslands or meadows. This is significant in the Fort Hunter- 
Liggett area, because most of the valleys are covered in grass, and all of the high 
ground is below the treeline. The result is a map with brown valleys and green 
ridgelines. While this was readily accepted as natural by most viewers, pilots 
with a background in low-level flight found it awkward, and contrary to their 
expectations (from flight charts) of green valleys and brown mountains. While 
this might be significant in other flight simulation applications (particularly those 
designed for pilots), the initial representation was deemed most appropriate for 
the target audience of the FOG-M simulator. 

A similar initial, intuitive choice was made for the elevation-keyed shading. 
High intensity (light) colors were used for higher elevations, and low intensity 
(dark) colors for lower elevations. This was accepted as natural by almost all 
viewers. The optimum number of intensities (shadings) to use in the map was 
experimentally determined to be sixteen. A small power of two was desirable due 
to the nature of the writemasks used to improve display speed. A large number of 
colors provides greater elevation definition and prevents large masses of the same 
color in areas where elevations change gradually. However, having too many 
colors destroys the contour-map effect, since adjacent colors are so close that no 
boundary is distinguishable between them. Hight shades each of green and brown 
were used initially. The shift to sixteen shades of each produces a better looking 
map. Due to the RGB (red, green, blue) nature of color creation on the IRIS, the 


greens were still barely differentiable at thirty-two shades, but the browns (a 
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combination of mostly red, some green, and, in some shades, a trace of blue) 
began to blend together. 

To determine the elevations at which color shades should change (in order to 
use the full range of shades), the maximum and minimum elevations of the 
terrain section in use must be known. Rather than preprocess the data before each 
run, these values are coded as constants in a header file. The equation for which 
color index to use is straightforward (see Equation 4.1) but takes significant time 


when repeated ten thousand times. 


elevation—MIN 
index = base index + —~———————. * #_ of shades (4.1) 
7 MAX-—MIN — 


Therefore, the fifteen points at which the shade changes are precalculated and 
stored in an array so that no calculations are needed at each point, just an array 


lookup. 


B. DRAWING 

The map can then be produced by determining the color and shade for each 
tile, and drawing it as a filled square. However, an increase in speed can be gained 
by exploiting the structure of the data and the line drawing hardware of the IRIS. 
The data is stiil processed a point at a time within each one kilometer column, 
but no drawing is done until an elevation/shading breakpoint is reached. Then a 
single line of one tile’s width is drawn to color all tiles since the previous elevation 


breakpoint. 
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C. WRITEMASKS 

A more significant speed improvement (on the order of fifty per cent more 
frames per second) was achieved with writemasks. Writemasks are a relatively 
low-level hardware feature that can be used for many purposes. In the FOG-M 
simulator, they are used to prevent the contour map from being overwritten. 
This allows the map to be drawn only once into the bitplanes, and have it remain 
on the screen without being re-drawn during each frame update. In order to 
understand how writemasks work, one must understand the layout and use of the 
IRIS’s color table and bitplanes. 

1. Color Table 

The color table associates a particular binary number with a color. 

When the display system asks what color some number is, the color table replies 
with the intensities for the red, green and blue color guns that will produce the 
color defined for the input number. This input number is referred to as a 
colorinder. Thus the color displayed on the screen depends on the colorindex 
associated with a given pixel, and the color associated with that colorindex in the 
color table. Table 4.1 gives the color table entries that are the defaults on the 
IRIS workstation. 

2. Bitplanes 

The colorindex that is associated with each pixel is stored in the display 

memory, which is composed of bitplanes. Each bitplane has one bit for each pixel 


on the display screen, so a bitplane is 1024 bits wide, 768 bits high and one bit 
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TABLE 4.1 IRIS DEFAULT COLORINDEX DEFINITIONS 



















Colorindex | 

Binary 

| Black | 0 ~~ _9000000000000000 | 
Red 1 OOO0000000000001 
Green 2 0000000000000010 
Yellow 3 0000000000000011 
Blue 4 0000000000000100 
Magenta 3 | OQO000000000000101 
Cyan 6 0000000000000110 
White 7 O0000000000000111 


deep. When used in doudle—buffer mode (as in FOG-M), the IRIS uses sixteen 
bitplanes (numbered 0 to 15) for each buffer. The frontbuffer is the one whose 
binary contents define the image being displayed. While the frontbuffer is being 
displayed, the next image is created by issuing drawing commands which affect 
only the backbuffer. Once a new image is completed in the backbuffer, the 
buffers are swapped, so the backbuffer becomes the frontbuffer and is displayed. 
The old frontbuffer becomes the backbuffer, and is then available for update. 
3. Writemask Example 

Consider the pixel at location (0,0) — the lower left corner of the screen. 
The colorindex of that pixel 1s determined by sixteen bits: one from the lower left 
corner of each bitplane. The display system reads those sixteen bits as a binary 
number (the colorindex), and uses the color table to determine what color to 


make that pixel. For example, using the default colors defined in Table 4.1 above, 


that pixel will be colored black if all sixteen bitplanes have zeroes in their lower- 
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left corners, since the value of the sixteen bit binary number 0000000000000000, is 
zero. If the current color is set to magenta (color five, whose binary value has ones 
in bits zero and two) and a drawing command is issued, bitplanes zero and two 
are set to one, and all other bitplanes are set to zero, for every pixel covered by 
the drawing command. These pixels will now be displayed as magenta, because 
the colorindex constructed from the sixteen bitplanes will be 0000000000000101, 
(Fn) and the color table tells the display system that color 5,, is magenta. 

The previous example showed that a drawing command works by 
placing ones in certain bitplanes, and zeroes in all of the rest, with the current 
color specifying which bitplanes get which. A writemask tells each bitplane to 
either allow or ignore the changes a drawing command says to make. In normal 
double-buffered usage, the writemask is 1111111111111111,, meaning all sixteen 
bitplanes should allow updates. Now suppose there is an image on the screen 
which uses just the default eight colors. Bitplanes three through fifteen are all 
zeroes, because all of the colors have colorindices with three or less binary digits, 
which will be in bitplanes zero, one, and two. If the writemask is changed to 
1111111111111000, after drawing the image, those lower three bitplanes are 
‘frozen’? and will not be changed by any drawing command. Setting the color to 
Diack and clearing the screen will not change anything. The upper bitp!anes will 
be set to all zeroes, which they already were. The lower three bitplanes will be 


told to reset to zero, but will not do it because they are protected by the 


writemask. 
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Now suppose you want to draw a grey line on top of the image. The line 
only needs one color, so it can be drawn in one bitplane. (Two bitplanes will allow 
three more colors on top of the map, three bitplanes allow seven, etc.) The first 
‘‘free’’ bitplane is number three. Turning on a bit in this plane (and not turning 
on any bits in higher planes) requires a colorindex in the range 1000, to 1111, (8,, 
to oe) Defining color eight in the color table as grey, making color eight the 
current color, and then drawing the line is sufficient to get the image into the 
bitplanes, but the display will not show the desired effect. If the image 
underneath the line is black (i.e. bitplanes zero through two are all zeroes form 
some pixels), the line will appear grey, as intended, for those pixels. However, if 
the image underneath the line is red (i.e. the lower bitplanes contain 001,), the 
composite colorindex retrieved by the display system is 0000000000001001, or 9,,) 
and since color nine is not defined in the color table, it appears as black. Thus 
every colorindex that has bit three (because the line is in bitplane 3) set to one 
(i.e. colorindices 1000, to 1111,, or 8,, to 15,,) must be defined as grey in order to 
produce the desired image. 

4. Writemasks in FOG-M 

The map image used in FOG-M is stored in the first six bitplanes 
(numbered 0 through 5) of both butfers, which means sixty-four colors are 
available: eight are the IRIS defaults, sixteen are shades of brown, sixteen are 
shades of green, and twenty-four are unused. The writemask defined as 


SAVEMAP (CO,, or 0000000011000000,) allows things to be drawn on top of the 
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map in bitplanes six and seven. Colorindices 64 through 127 are all defined as 
blue in the color table, so anything drawn in bitplane six appears on top of the 
map in blue. Similarly, bitplane seven is used for red, with colorindices 128 
through 255 all correspondingly defined to be red. 

The speed improvement due to writemasks in FOG-M comes from not 
having to re-draw the map each time the screen is updated. The cost is the use of 
Many more indices in the color table, which limits the number of colors available 
for use on top of the map. For our simulation system, with only two colors 
needed on top of the map, there is plenty of room in the color table. Therefore, 


the gain in speed comes at no real cost. 
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V. THREE-DIMENSIONAL TERRAIN CONSTRUCTION 


A. REPRESENTATION DECISIONS 
1. Polygons versus Patches 

Early experiments in the study involved attempting to display the 
terrain using parametric bi-cubic surface patches. A surface patch is simply a 
smooth curved surface fitted to a set of data points. A discussion of the theory 
and use of surface patches can be found in the IRIS User’s Guide [Ref. 7:sec. 11-3] 
and Hearn and Baker [Ref. 1l:pp. 193-205]. It was quickly determined that it 
would not be possible to use surface patches to represent the terrain and still 
maintain a real-time update of the terrain during flight. 

An alternate method of displaying a three-dimensional object is through 
the use of a set of planar polygon surfaces that join at common edges to form the 
terrain object. This method has the advantage of being much simpler, and 
therefore faster, to generate and display. For this reason it was chosen for use in 
the project. 

Figure 5.1 shows the method of constructing the terrain surface as a set 
of triangles. The term gridsquare is used in the remainder of the chapter to refer 
to a set of two triangles with a common hypotenuse that form a square of the 


terrain grid. 


34 


"GRID- 
SQUARE" 





eee: East 


View from 


above looking down on the terrain. 
See aia 1) 


ter. Onn 


elevation points are connected 
triangular polygons with common 
edges. 


pIcinewomderolygonal lérrain, Construction 


35 


2. Resolution 

The special DMA data file used in this project contains elevation data 
that is spaced at a twelve and one-half meter interval. One of the first questions 
which had to be answered concerning the three-dimensional! portrayal of this data 
was, ‘In how fine a resolution can the data be displayed, while still allowing for a 
sufficient frame update speed?” Early test runs showed that using the full twelve 
and one-half meter resolution would be much too slow, although it provided an 
excellent representation of the terrain. An adequate frame update rate 
(approximately three to four frames per second) was achieved with a seventy-five 
meter resolution or every sixth data point. Since this was an earlv test. displaying 
terrain without any targets or cultural features, a one cence meter resolution 
was decided upon for use in the remainder of the project. This allowed for an 
adequate “‘cushion”’ of processing time to complete the additional computations 
that would be needed in the final product, while still providing an adequate 
degree of resolution. 

3. Elevation Scaling 

After viewing the early representations of the terrain, it appeared that 
the hills did not give an appropriate appearance of height. Although this was a 
subjective judgement, it was shared by most people who viewed the display and 
compared it to aerial photographs of the area. Because of this, it was decided to 
scale the elevations of the displayed points upward. Two approaches, linear 


scaling and exponential scaling, were examined. 
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In the linear scaling approach, each elevation point was simply 


multiplied by a scale factor as shown in Equation 5.1. 


IC eno elev |, foe) 


new 


Using this approach, it appeared that a scaling factor between 1.5 and 2.0 was 
necessary to achieve the desired effect. 
In the exponential approach, the elevation of each point was raised to a 


fixed power as shown in Equation 5.2. 


Elev. = Elev. fo) 


ne 


This approach has the effect of exaggerating the higher elevations to a greater 
degree than the lower ones. It was chosen as the approach for use in the project 
based on subjective observations of the displays produced by the two methods. 
The scaling factor, o, was chosen as 1.05. Using this factor produces the 
equivalent of a linear scaling of 1.5 for the maximum elevation and 1.4 for the 
minimum elevation contained in our area of interest. 

Subsequent to the decision to use an exaggerated elevation scale, 
research results were discovered which supported it. In a study conducted by the 
U.S. Army Research Institute for the Behavioral and Social Sciences, observers 
were asked to pick a computer generated line drawing that best matched actual 
terrain. The line drawings had different exaggerations of the vertical (elevation) 


scale. The overall ratios chosen by the four observers ranged from 1.25:1 to 


oT 


1.50:1. The drawings presented to the observers had exaggeration ratios ranging 
from Labo! 1. focles eiasb2) 
4. Shading and Texturing 
As explained above, each one hundred meter square of the terrain, a 
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‘‘osridsquare,”’ is represented by two triangles in three-space that share a commen 
diagonal edge. The process of applying colors to these polygons, shading, was the 
next area of research in the project. 
a. Elevation Based Shading 
Three different shading algorithms were investigated. The first was 
a simple algorithm where the shade of a polygon was a function of its elevation. 
Higher elevations are shaded in lighter shades of green while lower elevations 


receive darker shades. Equation 5.3 represents the assignment of a shade from the 


color map. 


elev— Min Elev 
color index = base index + ————__ * #_ of “shades ang) 


7 7 Maz Elev— Min Elev 


The darkest green is stored in the base index color map location and the lightest 
green in the baseinder + # of shades location. Although this approach works 
well for two-dimensional contour maps (see Chapter IV), and is currently used in 
another “low cost” simulator Ref. 6j, it did not appear to present a realistic view 
of the terrain. An advantage of this approach, however, is that the calculation of 


the color index is simple enough to be done with no preprocessing. 
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b. Lambert’s Cosine Law Shading 
The second method of determining the shade for a polygon involved 
the use of a point light source and Lambert’s cosine law |[Ref. 11:p. 278]. Let N 
be a unit normal vector to the polygon, and L be a unit vector in the direction of 
the light source. The angle between N and i ®, is the angle of incidence. 
Lambert’s Law states that the intensity of the light reflected from the polygon is 


proportional to cos ® (Equation 5.4). 
Tacos ® (5.4) 


In order to use this law, the normal vector (V), the light source vector cay. and 
the angle between them (®) must be known. N can be determined by taking the 
cross product of vl and v2, where vi is a unit vector in the direction from vertex 
B to vertex C of the polygon, and v2 is a unit vector in the direction from vertex 


B to vertex A of the polygon (Equation 5.5 and Figure 5.2). 
N = v1 x v2 (5.5) 


With N and L available, cos ® can be computed as their dot product (Equation 


oak 
cose @=N-L (om 


Since the intensity is proportional to cos ®, the appropriate coior index to use can 


be computed as 


color index = min indez + (# shades*cos ©) (5.8) 
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Figure 5.2 Lambert’s Cosine Law 
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where min tnder is the color index of the lowest intensity green and 
min index + # shades is the color index of the highest intensity green. 
c. Gouraud Shading 

The final shading model investigated involved the use of Gouraud 
shading. The purpose of Gouraud ead ine is to provide a continuous transition of 
shades across a polygon so that the shades at the edges of adjoining polygons 
match. This in effect eliminates the visible boundary between polygons and 
provides a smooth continuous surface. The Gouraud algorithm involves 
interpolating to determine the intensity to be used at each pixel along a scan line, 
and is illustrated in Figure 5.3 as reproduced from Hearn and Baker [Ref. 11:p. 
290]. To use the algorithm, intensity values for each vertex of the polygon must 
be known. In the project’s implementation, the intensity at each vertex was 
computed as the average of the intensity values for all the polygons meeting at 
that vertex, where the individual polygon’s intensity values were calculated using 
Lambert’s cosine law. 

The use of this model posed two problems. First, even though the 
IRIS supports Gouraud shading in its graphics library, its use increased the time 
between frames to an unacceptable rate (approximately one and one-half to three 
seconds between frames). Second, the smoothing of the algorithm worked too 
well, resulting in terrain displays that lacked the necessary position cues to detect 
motion. This second problem could be alleviated by adding artificial texture to 


the terrain but in light of the speed problem it was not pursued further. 
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For interpolated shading, 
at point 4 is determined from intensity values 


intensity at point 6 is 


at points 1 and 2, 


determined from values at points 2 and 3, 
along 


and 


intensities at other points (sucasace) 


the scan line are interpolated between the 


values at points 4 and 6. 


Figure 5.3 The Gouraud Shading Algorithm 
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d. Adding Texture 

Lambert’s cosine law was chosen as the shading model for use in the 
project, providing the most realistic display within the allowed computation time 
constraints. However, a problem with its use is that the flat valleys, with little 
variance in the surface normals of their polygons, produce large geographic areas 
having a near constant shade. This results in a lack of motion cues in these areas 
similar to that experienced with the Gouraud shading model. To remedy this 
situation, a simple artificial texture, in the form of a checker board, was imposed 
on the terrain. The checker board effect was implemented as follows. First, the 
shades for the two triangles in each gridsquare were averaged, a this average 
shade was used for both of them. This of course causes the visible boundary 
between the triangles to disappear leaving a square shaded in a single color. 
Second, two slightly offset color ramps were used with adjacent grid squares using 
different ramps to compute their shades. One ramp is composed of green 
intensities ranging from 255 to 50, while the other uses intensities ranging from 
245 to 40.* This causes the shades for two adjacent gridsquares with identical 


surface normals to vary, providing the necessary texturing. 


*A value of 255 is the highest intensity green obtainable, a value of zero indicates the absence 
of the color green. 
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B. INTERNAL DATA STRUCTURES 

Two global arrays are maintained which store the information necessary to 
display the terrain. The first 1s a five-dimensional array, savetriangle, that stores 
the values of the coordinates for each triangle making up the terrain structure. 
The second is a two-dimensional array savecolor that stores the color map indices 
for each of the terrain’s grid squares. The purpose and range of each of 
savetriangle’s indices is shown in- Table 5.1. For example, 
savetriangle([3][5](1][1][2] would contain the value of the Y coordinate (fifth 
dimension = 2), of the second vertex (fourth dimension = 1), of the northern 
triangle (third dimension = 1), of the grid square with X index five and Z index 


three (second dimension = five and first dimension = three). 


TABLE 5.1 LAYOUT OF THE SAVETRIANGLE ARRAY 
Index Range 


Dimension Purpose 
' 


' First Grid square index in the Z direction. 0 
is the southern most square, 98 is the 
northern most. 





Second Grid square index in the X direction. 0 
is the western most, 98 is the eastern 
most. 

Third Wi 0. ih 1 I Triangle identifier within a grid square. 


| 0 is the southern triangle. 1 is the! 
| northern. | 


! 





| yu | ae Vertex number of the triangle. 0 is the 
first vertex, 2 is the last. 





Coordinate identifier of the vertex. O is 


the X coordinate, 1 the Y coordinate 
and 2 the Z coordinate. 
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Table 5.2 lists the purpose and ranges of each of savecolor’s indices. For 
example, savecolor(30|{10] contains the color map index to be used for the grid 


square with a Z index of thirty and an X index of ten. 


TABLE 5.2 LAYOUT OF THE SAVECOLOR ARRAY 


Index Range | 


| Dimension 


[First | Grid square index in the Z direction. 0 
is the southern most square, 98 is the | 
northern most. | 


Purpose 


Second Grid square index in the X direction. 
is the western most, 98 is the eastern 
most. 





These two arrays contain all the information necessary to constrict an image 
of the terrain. The following chapter provides the details of using their data to 
create a real-time, updated image of the terrain as it is seen from the FOG-M’s 


camera. 
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VI. FLIGHT SIMULATION 


A. OVERVIEW 

The previous chapter discussed the methodology of constructing the three- 
dimensional terrain from the provided elevation data. This chapter’s purpose is 
to explain the details of displaying this terrain in real time as it is seen through 
the missile’s camera. 

The high level pseudocode for the main program’s terrain display loop is 
shown in Figure 6.1. Chapter VII explains the details of step two. The details of 
steps one and six are explained in Appendix B under the procedures readcontrols 
(for step one) and edit navboz and edit indboz (for step two). The remainder of 


this chapter discusses the details, considerations, and results of implementing 


steps three through five. 


B. UPDATING THE MISSILE’S POSITION 


Determining the missile’s new position can be broken into two cases: 


[1] the missile is under operator control and its new position is a function of the 
old position, the commanded direction of flight. the commanded altitude. 
and the commanded speed. 


2] the missile is locked onto a target and its new position is a function of its oid 
position, the position of the desired target, and the commanded speed. 


In both cases, a very large simplifying assumption is made to ignore the 


dynamics of the missile’s flight. This means that the missile is able to 
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While missile is flying do 
1) Read the values from the operator’s controls 
2) Determine new positions for all the targets 
3) Determine the new position for the missile 
4) Determine the position of where the camera is looking 
5) Display the terrain as seen by the camera 
6) Update the operator's control indicators 
End while 


Figure 6.1 Main Display Loop Pseudocode 


instantaneously change heading, speed, and altitude. This assumption was made 
only because of development time constraints. It is felt that the computations 
necessary to more realistically model the dynamics of the flight can be done 
without a serious degradation of the simulator’s performance. 
1. Case 1 - Operator Control 
Under this case the missile’s X,Y, and Z coordinates are computed as 


shown below. 


ADist = Speed*A Time (6.1) 
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Where 


- ADist is the distance traveled over the ground since the last position was 
calculated. 


- Speed is the missile’s speed in feet per second and 


- ATime is the elapsed time since the last position was calculated ~ 
Having calculated the distance the missile must move during this frame the 


missile’s new coordinates (MX,MY,MZ) can be calculated as 


MX,,,,, = MX,,,+[cos(Dir,_4) *ADist| (6.2) 

MZ_,,., = MZ,,,—|sin( Dir. ,) *A Dist] (6.3) 

VAY = Alien | (6.4) 
Where 


- Dir, 8 the commanded heading in radians 
- Alt. is the commanded altitude in feet 


- a is the altitude scaling factor (see Chapter V, Section A.3). 
2. Case 2 - Locked Onto a Target 
In the case where the missile is locked onto a target, the missile’s new 
position is computed as follows. A Dist is computed as in Equation 6.1. Next the 
missile’s heading is computed so as to steer it directly toward the target’s 


position: 


Dir, , = arctan2(-—(TZ—MZ|,{| TX—X)) (6.5) 


gt 


48 


Where 


E Dir, is the direction from the missile’s position to the target’s position 


TX is the X coordinate of the target’s position 


TZ is the Z coordinate of the target’s position 


- MX is the X coordinate of the missile’s position 


- MZ is the Z coordinate of the missile’s position 
° a * 
- arctan2(a,6) is a function which returns the arctan} — {| in the range 
b 
0 to 2II, based on the sign of a and b. 


Once Dir, is known, the missile’s new X and Z coordinates can be calculated as 


MX, .y = MX, )4t+|¢0s(Dir,.,) *A Dist] (6.6) 
MZ, = MZ,,,—|stn(Dir,,,) *A Dist] (6.7) 


Next the missile’s altitude (MY) is adjusted a proportion of the total altitude 
difference between it and the target, based on the ratio of ADist to the total 


distance (along the horizontal plane) to the target. 








Dist,,,=V (TX-MX)°+(TZ-M2Z) (6.8) 
f A Dist 

MY ee! * oe eet) : (6.9) 
Dist, 





Where 


- Dist, is the distance to the target measured along a horizontal plane. 


- MY and TY are the Y (altitude) coordinates of the missile and target, 
respectively. 
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C. DETERMINING THE LINE OF SIGHT 

Once the new position of the missile has been calculated, the next step in 
displaying the terrain is to determine another point along the camera’s line of 
sight: the look-at position. This calculation is also broken into two cases based 
on whether the missile is or is not locked onto a target (see Figure 6.2). 

The case where the missile is locked on is trivial, the look-at position is 


simply set to the coordinates of the locked-on target. 


LX=TX (6.10) 
LY=TY (6.11) 
LZ=TZ (6.12) 


Where LX, LY, and LZ are the X, Y, and Z coordinates of the look-at position. 
This centers the target in the displayed three-dimensional scene. 

When the missile is not locked onto a target, the camera’s look-at position is 
a function of the missile’s position, the missile’s heading, and the pan and tilt 


angles of the camera. It is determined as follows 


Dir, , = Head, ,+ Pan (6.13) 
LX = MX+|cos(Dir, ,)* Dist, 1] (6.14) 
LZ = MZ—|\sin( Dine oe (6.15) 

LY = MY+([Dist,, *tan(Tilt)] (6.16) 
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Case 2 - Missile Not Locked on Target 


Figure 6.2 Determining the Camera’s Look-at Position 
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Where 


- Dir), is the direction the camera is looking 
- Pan is the pan angle of the camera 
- Tilt is the tilt angle ot the camera 


- Dist,,,, is an arbitrary distance over the ground that the camera looks ahead. 
Since the only purpose of LX, LY, and LZ is to determine a point along the 
camera’s line of sight, any positive number will be acceptable. A value of five 
kilometers is currently used. 


DD, DISPLAYING FOE SCR ME 
Once a line of sight has been determined, the next steps are to apply the 
appropriate viewing transformations, draw the filled polygons that make up the 
terrain, and add other items to the scene such as targets and roads. 
1. Viewing Transformations 
It is possible to project a three-dimensional object onto a two 
dimensional viewing surface in two basic ways. In one method, the parallel 
projection all the points of the object are projected along parallel lines. This has 
the advantage of preserving the relative dimensions and angles within an object 
and is used when accurate views of various sides of an object are needed such as 
in architectural drawings. In the other method, the perspective projection, all 
the points of an object are projected along lines that converge at a single point 
called the Center of Projection. In this method, relative dimensions are not 
preserved. Lines closer to the projection plane appear larger than those that are 


more distant. The perspective projection provides a view of three-dimensional 
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objects that is more realistic, similar to that provided by the human eye or a 
camera. Both these projections are illustrated in Figure 6.3. [Ref. l1l:pp. 235-241] 

Because of its more realistic presentation of the scene, a perspective 
projection was used for the project’s three-dimensional scenes. The IRIS’s 
graphics library provides a procedure called perspective which constructs the 
necessary transformation matrix* to obtain a perspective projection. The matrix 


is defined as [Ref. 7:p. C-2] 


Perspective (fovy,aspect ,near,far) = 








fovy 
cot ( ) 
Z 
0 0 0 
aspect 
fovy 
0 cot ( 0 0 (6.17) 
2 
far+near 
0 0 —_——————_ -l 
far—near 
2x far x near 
0 0 SS 


far—near 


Where 


- fovy is the field of view angle 


- aspect is the aspect ratio, a ratio of the distance a viewer sees in the X 
direction to the distance he sees in the Y direction. It is generaily set to be 
the same as the ratio of the width to the height of the viewport. 

- near and far are the distances from the viewer to the near and far clipping 
planes. 


*A knowledge of using transformation matrices to perform graphical operations is assumed 
here. Hearn and Baker [Ref. 11:chaps. 11-12] provides excellent coverage of the subject. 
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Figure 6.3 Parallel and Perspective Projections 
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The perspective projection forms a view frustum as shown in Figure 6.4. 
Any object within the frustum between the near and far clipping planes will be 
displayed in the scene. Objects outside this view volume are clipped and 
discarded. 

Next, the frustum formed by the perspective projection must be 
positioned along the camera’s line of sight. This is accomplished by another 
transformation matrix constructed via a graphics library procedure named lookat. 


The lookat procedure takes the following inputs: 


. M3 and V_: the X, Y, and Z coordinates of the center of projection. 
ae ee and P_: the X, Y, and Z coordinates of the look-at position. 


Twist, a right handed rotation of the scene about the line of sight. 
The transformation matrix formed by lookat is actually the result of multiplying 


four other transformation matrices [Ref. 7:p. C-2] 


Lookat ( ay Vy Ee,» Twist) = 


| 6.18 
Trans(- ie a a. V.) * Rot (©) x Rot (®) x Rot (— Twist) | | 


Phere Trans(—V_,— Vi V)= (6.19) 
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Clipping 


Planes 





The perspective command defines a near and 


far clipping plane, a field of View ame 


Figure 6.4 The Perspective Command 
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Rot, () = (6.20) 


1 0 0 0; 
0 cos(®) sin(®) O as 
Rot,(®) = 0 —sin(®) cos(®) 0 ean 


0 0 0 1 


cos(—Twist) sin(—Twist) 0 0 


—sin(—Twist) cos(—Twist) 0 0 


Rot,(— Twist) = (6.22) 
0 O 1 O 
0 O O 1 
La ve 
And © = sin (6.23) 
Ce ees aN) 
VSP 
. =i yy 
® = sin (6.24) 





As can be seen, this transformation simply translates the center of projection to 
the origin, then rotates the view frustum about X and Y axes to align with the 


line of sight. Finally the twist angle is added with a rotation about the Z axis. 
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In the flight simulation, the twist angle is analogous to the “roll’’ angle of an 
aircraft or missile. A value of zero is currently used, but other values could be 
used if the roll of the missile during flight was added to the model. 

2. Determining Which Polygons to Draw 

After the correct viewing transformations have been applied, the 
polygons that comprise the scene must be drawn. Although the IRIS will ‘“‘clip” 
polygons which lie outside the perspective projection’s view volume, an increase in 
frame update speed can be achieved by not attempting to draw those that 
obviously lie outside. This is discussed further in the following section on 
simulator performance. 

The term view—bound is used to describe a north-south oriented 
bounding box around those parts of the scene that are sent to the graphics 
pipeline. The view-bound is described by the index of the northernmost, 
southernmost, easternmost, and westernmost gridsquare to be drawn. It is 
calculated by extending (if necessary) the line-of-sight vector until it intersects the 
horizontal plane Y = Min elev, where Min elev is the minimum elevation value 
of the terrain. The view—bound is calculated as being 20 gridsquares to the 
north, south, east, and west of this intersection point. If the missile’s X and Z 
coordinates are not within the calculated view—bound, the bounds are extended to 


include them. Figure 6.5 illustrates this construction. 
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East View-bound 


1) Line of sight vector is extended down 
to intersect the minimum elevation plane. 


2) View bound extends 20 gridsquares north, 
SOutm, east and west of the intersection: 


3) Bound is extended, if necessary to 
include the missile’s position. 


Figure 6.5 Construction of the View-bound 
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3. Hidden Surface Removal 

A final detail that must be taken care of is the removal of hidden 
surfaces from the scene. A hidden surface is simply a part of the scene that is 
obscured by some object in the foreground, such as a valley that it hidden behind 
a large hill. 

The IRIS supports a method in hardware called Z-Buffering. In this 
method, a buffer is maintained for each pixel position on the monitor and 
contains the ‘‘depth” (transformed Z coordinate) of the part of the scene that 
generated that pixel. Before drawing is started, the buffer is initialized to the 
maximum depth value (the value of the far clipping plane) for each pixel position. 
Before each new pixel is drawn, its depth is compared to the depth stored in the 
buffer. If its depth is greater than the stored depth it is not drawn. If it is less 
than the stored depth, it is drawn and its depth value replaces the value in the 
buffer. This method could not be used in the project. for two reasons. First, with 
comparisons having to be made on a pixel-by-pixel basis, it slows down the frame 
update rate to an unacceptable level. Second, the IRIS does not allow the use of 
Z-buffering and double-buffering at the same time. Double-buffering is necessary 
to implement the animation of the scenes. 

Another common method of hidden surface removal is the painter’s 
algorithm. It derives its name from the way a painter would draw a scene on 
canvas, drawing in all the background and then adding foreground objects by 


painting over the background objects they obscure. Implementing this algorithm 
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in computer graphics means drawing the scene in an ordered fashion, such that 
the most distant objects from the viewer are drawn first and those closest to the 
viewer are drawn last. Since the gridsquares comprising the terrain form well 
defined rows and columns, an efficient implementation of this algorithm is 
possible. That implementation is described below. 

The implementation can be thought of on a conceptual level as follows. 
A line, perpendicular to the line-of-sight, is constructed to serve as a pseudo- 
scanline. Gridsquares within the view-bound are drawn as they are intersected by 
this scanline. The scanline is first positioned along the line-of-sight vector so that 
it intersects the far corner gridsquare of the view-bound. After all the gridsquares 
along the scanline have been drawn, it is moved one gridsquare closer to the view 
position, along the line-of-sight vector, and the process is repeated. This 
continues until all the gridsquares within the view-bound have been drawn. 
Figure 6.6 illustrates this process. 

From Figure 6.6, notice that each scanline passes through three 
gridsquares in a column, shifts over a column, then passes through three 
gridsquares in the next column. The number of gridsquares drawn in a column 
(or row) before advancing to the next column (or row) can be determined by 
computing the tangent or the scanline’s direction. If the magnitude of the 
tangent is greater than 1.0, scanlines will run and shift along columns of 
gridsquares. If it is less than 1.0, scanlines will run and shift along rows of 


gridsquares. The term threshold is used in the remainder of the algorithm to 
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Figure 6.6 The Scanline Hidden Surface Algorithm 
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describe the number of gridsquares drawn before a shift of column (or row) takes 


place. It is computed as 


I 
nearest integeritan(Dir.___ ) if jtan({Dir _)| 21.0 














threshold = (6.25) 
(tan(Dir,_)) , if |tan( Dir 150 


nearest integer mm) 














The pseudocode for implementing the algorithm is shown in Figure 6.7. 


The case shown is for a line-of-sight direction that is in the first octant (between 


O and — radians). The algorithm for the other seven octants is similar, the 
4 


difference being the direction the scan line advances, and the direction it shifts 
when the threshold is reached. Table 6.1 summarizes these parameters for all 
eight octants. 


TABLE 6.1 VARYING PARAMETERS FOR THE SCANLINE ALGORITHM 
BASED ON THE OCTANT OF THE LOOK DIRECTION 


Look Directions |[ Scan Line Advances |] _ 


When Threshold is Reached 





ro | W/4 [Noh [~ South || Shift one column Bast 
Hija |_W/2 [| Bast_ [West || Shift one row North 
| 1/2 | 3it/2 | West | Bast || Shift one row North 


2 
3 
es | See eeConthiee| PShifronc column West 
a ee ome ae Nothel fieShiftonecolumn West | 
6 | | 


bidy eS | atl) 2 | West | Kast | Shirt one row South 
















| mot ew “aT? ihe Pease | West | Shift one row South 


8 M4 | 2M |] South | North {| __Shift one column East__ 


Notice the step draw gridsquare/z indez//xz indez] in the algorithm. 





Since a gridsquare contains terrain, and can also contain roads and targets, an 
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Calculate the threshold value 
count < 0 


start x index <— west view bound 
start 2 index — north view bound 


While start z index > south view bound do 
z index <— start z_index 
x index <— start x index 


while (x_index < east view bound) and (z index 2 south view bound) do 
{ traverse a scanline } 
draw gridsquare|z_ index||x index| 
z index — 2 index - 1 {move it one gridsquare south} 
count < count + l 
if count = threshold then 
x index — x index + 1 {move it one gridsquare east} 


count + 0 { reset count} 


endif 
end while 
{move on to next scanline: start it one gridsquare to the west} 
start x < start x- 1 
count <« 0 
if (start x < west view bound) then 
start x < west view bound 


start z + start 2 - threshold 


endif 


endwhile 


Figure 6.7 Pseudocode for the First Octant Scanline Algorithm 


ordering of these parts of the gridsquare must also take place. The two triangles 


forming the terrain are drawn first, next any roads are drawn, and finally any 
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targets are drawn. The details of integrating the targets and roads into the scene 
are covered in the following two chapters. 

The resulting scene is shown in Figure 6.8, a photograph of the IRIS 
monitor during the flight simulation. Note how the hidden surface removal allows 
the foreground hills to naturally obscure the valleys behind them. Also note the 


effect of the lighting model and texturing described in Chapter V. 


E. SIMULATOR PERFORMANCE 

Data collected while running the simulator shows that the average frame 
update rate is approximately four frames per second. The Unix profile utility 
was used to determine which procedures accounted for the majority of the 


simulator’s time usage. Table 6.2 shows the results for the top four routines. 


TABLE 6.2 FOG-M ROUTINES USING THE MOST CPU TIME 
% CPU Time | Routine Name —— Purpose 


[16.9 | polf | Iris graphics library filled polygon routine. 
display terrain | Output 3-D scene with hidden surface removal. 


8.7 malloc C language built in routine for dynamic 
memory allocation. 





eal gl findhash Low level Iris graphics library routine, used for 
| the hash tables associated with graphical 
objects (Not user accessible). | 


The top two entries in Table 6.2 are directly involved with outputting polygons to 
build the terrain image. It is therefore reasonable to believe that the frame 
update rate depends heavily on the number of polygons that are passed to the 


geometry engines. 
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Figure 6.9 is a scatterplot showing the frame update speed achieved when 
various numbers of polygons were attempted to be drawn. The data was 
generated by reading the system clock before each frame update and calculating 
the number of polygons based on the vwiew—bound that was used during that 
frame. The graph clearly shows the effect the view-bound has on the frame 
update rate. The next two entries, malloc and gl findhash, are traceable to the 
making and deleting of the graphical objects that store the targets (this process is 
explained in Chapter VII). As an experiment, the construction and deletion of 
the targets’ objects was removed from the simulation and the targets were simply 
displayed in stationary positions. The profile results from the simulator run in 
this configuration is shown in Table 6.3. Figure 6.10 is another scatterplot, 
generated in the same manner as Figure 6.9, except that the simulator was run in 
the stationary target configuration. Eliminating the dynamic memory 
Management associated with the target’s graphical objects increased the average 
frame update rate from 2.99 to 3.90 frames per second. Also, the maximum frame 
update rate achieved doubled from 7.5 to 15.0 frames per second. This would 
suggest that an area for further research is an improved algorithm for target 
updating that does not involve dynamically allocating memory. 

The tact that the frame update rate is so heavily dependent on the number otf 
polygons passed to the geometry engine suggests that a more sophisticated 
method of determining the view-bound may pay off in increased performance. 


For example, the present method does not take into account the field of view 
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angle. It should be possible to bound the line-of-sight intersection point with less 
than twenty grid squares when the field of view angle is small. However, any new 
algorithm developed can not be so sophisticated that it negates the performance 


increase by requiring intensive computations. 


TABLE 6.3 FOG-M ROUTINES USING THE MOST CPU TIME 
oe STATIONARY TARGETS) 


% CPU Time 
Beso) mmo) CC Ti Iris inca eines eniicdm als eonioutinen nal library filled polygon routine. 


_ 5 Hae ae terrain | Output 3-D scene with hidden surface removal. | 3-D scene with een eUNEEe removal. 
current drawing color. 

line intersect2 | Finds the intersection point of two lines. This 

routine is used exclusively during the road 


building process and therefore is not used at 
all in the display loop. 


poly Iris graphics library unfilled polygon routine. 
Used during the display loop to outline the 
road segments. 
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Vil. TARGET INTEGRATION 


A. GENERAL 

The primary targets of a FOG-M missile are tanks, helicopters, and 
reinforced ground installations. The simulator is designed to handle many types 
of targets, including various tanks and helicopters, but only a single type of tank 
is currently implemented. The prototype simulator provides an Ethernet 
networking capability to allow the input of actual target positions in real-time. 
This simulates the input that would be received by a production eralatoe during 
computerized mock combat field experiments. In its networking mode, the 
simulator receives target position and orientation data from an interactive 
program running on a different IRIS workstation. The target program, still in 
testing and not detailed in this study, provides the capability to dynamically 
insert and delete targets at any location, and to modify their speed and direction. 
In the simulator’s stand-alone mode, there are ten tanks defined by default that 
criss-cross the ten kilometer square terrain area. These tank targets move at a 
constant speed of fifteen knots and reverse direction when they reach one of she 


edges of the ten kilometer terrain square. No automated path planning is 


presently performed in either mode, so the tanks blithely traverse even the 
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steepest terrain. The default targets minimize this problem by traveling the length 


of the valleys for the most part. 


B. TARGET CREATION 
Target creation is simplified through the use of graphical objects. The actual 
image of a tank is defined initially by the tedious specification of the three 
coordinates of each vertex of each of the polygons that comprise the tank (Figure 
7.1). Using objects, this need only be done once, placed in an object, and then 
referred to by a single name within each target object. Thus each target is 
described by an object (the tank object) within another object (the target object). 
In addition to the tank object, the target object also contains the transformation 
commands that move the tank from the origin to its location on the terrain (a 
translation), and face it in the direction it is moving (a rotation). 
1. The System Matrix 
The rotation and translation commands work by modifying the system 
matriz. The system matrix is a global data structure that is used to transform 
coordinates from the three-dimensional world space into the two-dimensional 
screen space. Each transformation can be performed as a series of computations 
on individual Y.Y,. and Z coordinates, but the transformations can also be 
accomplished with a single matrix multiplication. The IRIS has a matrix 
multiplier built into its hardware, so matrix operations are very efficient. At least 


three transformations must be applied to every endpoint on the tank: a coordinate 
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scaling, a translation, and a rotation. Rather than do three separate matrix 
multiplications, the three transformation matrices can be combined, so that all of 
the transformations are accomplished in a single matrix multiplication. The 
matrices are combined by applying each of a to the system matrix. Each 
point is now completely transformed through a single multiplication with the 
system matrix. When a new transformation is needed, the system matrix must be 
reset by applying the inverses of the old transformations, or by copying the 
original contents back into the system matrix. Two commands are provided with 
the IRIS to support the latter method. Pushmatriz takes a copy of the system 
matrix’s current contents and saves it on the system stack. After the 
transformations have been applied, and the drawing that used _ those 
transformations has been completed, the system matrix is reset by calling 
popmatriz, which retrieves the copy placed on the stack by pushmatrix and 
restores the contents of the system matrix to the previously saved values. 
2. Target Transformations 
The tank is initially defined with its center interior at the origin 
(coordinates (0,0,0)). While it is not important which point on or in the tank is 
placed at the origin, it is crucial that the tank be defined somewhere around the 
origin in order for the rotation command to have the desired etfect. The originai 
direction of the tank is significant only to the extent that it must be known in 
order to calculate the appropriate rotation to achieve a specified heading. The 


tank in FOG-M faces to the right (zero radians mathematically, or a compass 
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heading of ninety degrees) initially. During target creation, dummy (zero valued) 
rotation and translation commands are placed in the target object, to be updated 
for display by a later editing of the object. Since all rotation and translation 
commands affect the system matrix (as previously described) and are cumulative, 
each target object must apply its transformations, be drawn, and then remove 
those transformations so that latter drawing commands are not distorted. Within 
each target object, the contents of the system matrix are saved with a pushmatrix 
call, the appropriate rotation and translation commands are applied to the system 
matrix (in reverse order, due to the nature of matrix multiplication), the target is 
drawn by calling the tank object, and then popmatrix is called to reset the system 


matrix. 


C. ANIMATION 

Animation of the targets is accomplished using the objects and 
transformations described above. The targets must be moved slightly before 
being redrawn in the next frame. This requires new (X,Y,Z ) coordinates, from the 
network or from local calculations. Then a global data structure is updated to 
indicate when in the display algorithm the target should be drawn, and the 
translation command in the target object is edited to provide the new coordinates. 
As each frame is displayed, targets appear in slightly shifted positions, and give 


the appearance of animated motion. 
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The calculation of new coordinates requires the maintenance of position, 
speed, and direction data for each target. The total distance traveled between 
screen updates is the product of the elapsed time (obtained from the IRIS’s real- 
time clock) and the target’s speed, scaled so the units match. In the networking 
version of the simulator this is done remotely; in the stand-alone version 
everything must be maintained locally. The target’s direction of travel is stored 
in radians, and is measured using the standard mathematical convention as 
opposed to a compass heading (Figure 7.2). This allows calculation of the the 


appropriate east/west (AX) and north/south (AZ) movement as follows: 


AX =cos(direction) “time “speed * scale factor (78 


AZ=—sin(direction) * time * speed * scale factor (72h 


The new target (X,Z) position is the sum of the old position and the offsets 
(AX,AZ) from Equations 7.1 and 7.2. Since all of the current targets are tanks, 
their Y coordinates (altitude) should be taken from the height of the terrain 
underneath the tank. This is obtained from the DTED interpolation routine 


gnd_ level, which is called with the new (X,Z) coordinates as input parameters. 


D. DISPLAY 

Chapter Five explained the exploitation of the structure of the data and the 
use of the painter’s algorithm to solve the polygon ordering problem without 
resorting to slower or more complicated schemes like Z-buffering or Binary Space 


Partitioning [Ref. 13]. Targets cannot merely be drawn after the terrain because 
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of the same ordering problem. Otherwise, targets appear in front of everything, 
and it is impossible to simulate a target moving out of sight into the distance or 
behind some terrain feature. The implementation of the target display algorithm 
is greatly facilitated by the use of objects. Objects allow the grouping of drawing 
commands into a subroutine-like package, which can be edited (effectively 
allowing parameterization) and then displayed with a single command. A two- 
dimensional array of object ‘“‘names”’ (the object—name-—array) is initialized so 
each element of the array represents the target object to be drawn in the one 
hundred meter square of terrain with the same indices. Since the C programming 
language recognizes the value integer zero as FALSE, and anything else as TRUE, 
this APRA does double duty as an array of booleans indicating the presence or 
absence of a target object in a particular one hundred meter grid square. (No 
target objects are given the ‘‘name” zero, which would indicate FALSE.) A list of 
targets is used to reset this array to all zeroes before each screen update (i.e. only 
those elements that contained targets need to be zeroed) so maintenance overhead 
of the array is minimized. The new target positions are received over the 
network, or are calculated, based on each target’s position, speed, and direction, 
plus the elapsed real-time since the last update. The appropriate object-name- 
array indices are calculated from the new target position and the object-name- 
array is updated. If this is the first (or only) target in the designated one hundred 
meter grid square, the update is accomplished by making a new object, and 


setting the object-name-array element equal to the new object’s integer ‘‘name.”’ 
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If the array shows that some other target is already in that particular piece of 
terrain (i.e. the object-name-array element is non-zero), the current target is just 
added to the object specified by the ‘‘name’”’ in the array. Once this has been 
done for each target, this array is available for the display terrain module. 
Display terrain checks the array as it draws each square of the terrain to see if 
any targets should be drawn. If so, it calls the indicated target object just after it 
has drawn the one hundred meter grid square on which the target(s) rests. Note 
that this causes the target(s) to be drawn at the correct time for the painter’s 
algorithm. The correct place to draw the target still must be specified by the 
transformation commands within the target object. 

In some cases it is necessary to draw a target more than once. Targets that 
straddle a one hundred meter grid square boundary must be drawn on top of both 
(or possibly all four) grid squares in order to avoid being partially obscured by 
whichever grid square is drawn last. (The target must be drawn immediately after 
the grid square on which it rests to ensure that the target will be obscured when 
it should be, by terrain drawn in the foreground.) Since the calculation of 
boundary intersection involves several trigonometric functions and an allowance 
for the distance between the center of the tank and its boundaries (which varies 
with the direction of the tank). a simplifying algorithm is used. If the tank is close 
enough to a boundary that the most distant part of the tank might cross the 
boundary (see tanks A and B in Figure 7.3), the target object is also drawn after 


the adjoining grid square(s). 
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The one hundred meter grid square is essentially divided into three areas: 
the middle, its sides, and its corners. In the middle, the tank cannot overlap any 
other grid square. On the sides, the tank may overlap one adjoining grid square, 
and in the corners, the tank may overlap three adjoining grid squares. The 
reference point on the tank (the position the X,Y, and Z coordinates refer to) is 
located at the very center of the tank. The tank is thirty feet long, so the most 
distant parts of the tank are within a fifteen foot radius of the tank’s reference 
point. The lines that mark the side and corner areas are thus fifteen feet inside the 
borders of the grid square. Once the tank’s reference point is within these areas, 
it is potentially obscured by the later drawing of the adjacent grid square(s). It 
might not be obscured if it is paralleling a side, for example, but the overhead of 
drawing it twice (or even four times) when it does not need to be is smaller than 
the overhead of the calculations to determine if the position and direction of the 
tank have it actually crossing one or more edges. 

The repeated drawing is accomplished by adding a ‘‘new”’ target to the array 
of target objects. The “‘new”’ target object is drawn at the exact same location in 
the three-dimensional terrain, but it is drawn after a different one hundred meter 
grid square. so it will have different target object array indices. and be in a 
separate target object, even though the two (or four) targets drawn will overwrite 


each other and produce a single image. 
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VIII. CULTURAL FEATURE INTEGRATION 


The addition of cultural features add much to the realism of the displayed 
scene. They also provide valuable landmarks from which a person observing the 
scene can geographically orient himself. This chapter covers the addition of one 
type of cultural feature, roads, fe the FOG-M simulation. Roads were chosen as 
the first feature to add because of the special problems associated with their 
implementation, the ease of extracting their locations from contour maps, and the 
visual impact added to all parts of the scene due to their wide-ranging locations. 
Three areas will be discussed: (1) the format of the external data file that contains 
the road’s locations, (2) the process of mapping the roads onto the existing 


terrain, and (3) the integration of the roads into the terrain display loop. 


A. EXTERNAL DATA FILE FORMAT 

The data being used in the simulation was obtained by manually extracting 
the roads’ positions from a DMA Topographic Center (DMATC) contour map of 
the area. Although this data is available in the DMA’s Digital Feature Analysis 
Data (DFAD) file, the software necessary to access it was not available. The road 
data file’s format is such that the DF AD data can be easily used when the access 


software is developed. 
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Figure 8.1 shows a segment of the file containing data for two roads along 
with a diagram showing their locations within the terrain. Each road entry is 
composed of three parts. The first part is the width of the road in feet. Next is 
an integer N, where N is the number of data points used to digitize the road. 
Third is a set N coordinate pairs, where each pair represents the location of a 
digitized point along the road’s centerline. The first coordinate of the pair is the 
east-west location of the point. It is measured in feet from the western terrain 
boundary. The second coordinate of the pair is the north-south location of the 
point, measured in feet from the southern terrain boundary. All the data is stored 
as ASCII text, which facilitates editing of the data using any text editor. The 
DFAD data file also contains road width information (in meters) and stores roads 
as a series of digitized points. The major difference is that DFAD’s points are 
stored as latitudes and longitudes, which need to be converted before they can be 


used in the simulation. [Ref. 9] 


B. CONSTRUCTION OF THE ROAD POLYGONS 

Knowing the width and centerline locations for the road, the next step is to 
construct the polygons which represent it. Although, this seems like a simple 
procedure, it is complicated by the fact that the road must follow the rise and fall 
of the terrain. Also, in order for hidden surface elimination to occur, the road 
must be divided at the gridsquare boundaries so that each piece can be drawn 


along with its corresponding gridsquare. The result is that the road must be 
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broken into many planar polygons, where each polygon is a portion of the road 
that overlays one of the terrain triangles within a gridsquare. Figure 8.2 
illustrates this division and defines some of the terms used in the description that 
follows. The high level pseudocode for processing the road data and constructing 
the planar polygons is shown in Figure 8.3. As the pseudocode shows, each road 


is processed a segment at a time. For each segment 


- The end points of the segment’s left and right side are calculated. A look- 
ahead to the next road segment is done, allowing the ends of adjacent 
segments to be calculated so that they meet cleanly. 


- A bounding box, which contains all the gridsquares intersected by the 
segment, is constructed. 


Next, for each gridsquare in the bounding box, the road segment is divided into 
the road-polygons at the gridtriangle boundaries. Note that all the vertices of the 


road-polygons fall into one of five types: 


- The intersection of a segment’s left side with the side of a gridtriangle. 
- The intersection of a segment’s right side with the side of a gridtriangle. 
- A gridsquare’s cornerpoint that is contained within the road segment. 

- An endpoint of the left side of the road. 

- An endpoint of the right side of the road. 


The road polygon is constructed by finding all the above vertices which exist, and 
ordering them counterclockwise. The counterclockwise ordering is necessary for 
backface polygon removal to take place. The intersections only derine the .Y and 
Z coordinates or the vertices. The Y (elevation) coordinate is found by 
interpolating between the terrain’s elevation at the three corners of the 


corresponding gridtriangle. 
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Figure 8.2 Constructing the Road Polygons 
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While more data in the road data file do 
read width of road 
read number of points 
read segment’s start coordinate pair (seg _start) 
read segment’s end coordinate pair (seg end) 
for i = 3 to number of points + 1 do 
if i < number of points then 
read the next segment’s end coordinate pair (next seg end) 
else 
next seg end x ¢ seg end x 
next seg end z *+ seg end 2 


endif 


calculate the start and end points for the segment’s left and right side 
(left start, left end, right start, right end) 


calculate a bounding box around the road segment 
for each gridsquare within the bounding box do 


Construct the polygon which overlays the gridsquare’s northern triangle 
Add the polygon to the road object associated with this gridsquare 


Construct the polygon which overlays the gridsquare’s southern triangle 
Add the polygon to the road object associated with this gridsquare 
right start < right end 


endwhile 


Figure 8.3 Pseudocode for Constructing Road Polygons 


C. INTERNAL ROAD-POLYGON STORAGE 

A global, two-dimensional array of graphicalobjects, named road, is used to 
store the road polygons. Each entry in the array corresponds to the pieces of road 
that lie within a gridsquare. An object is created when the first road-polygon is 
constructed for a gridsquare, with subsequent road-polygons being inserted into 


the already existing object. Since the roads are static in nature, the use of objects 
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does not present the dynamic memory allocation problems associated with their 
use in storing targets (see the Simulator Performance Section of Chapter VI). As 
each gridsquare of the terrain is drawn, a check is made to see if a road object 
exists for that square. If one does exist, the associated road-polygons are drawn 
immediately after the terrain. This insures that hidden surface elimination occurs 
for the roads as well as the terrain. A photograph of terrain which includes some 


sections of roads can be seen in Chapter VII, Figure 7.1). 
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IX. FOG-M SIMULATOR USER’S GUIDE 


A. OVERVIEW 

This section of the report is a user’s guide to running the FOG-M simulator. 
The simulator was built to be largely self documenting. Instructions are clearly 
displayed on the screen, including diagrams which serve as a reminder of the 
functions of the various controls. A knowledge of the logon procedure for the 
IRIS workstation and the basic commands of the UNIX operating system is 


assumed. 


B. STARTING THE SIMULATION 

To start the simulation, logon to the IRIS workstation and use the UNIX cd 
command to change to the directory containing the simulation. Currently the 
simulation is in the directory /work/terrain. Therefore issue the command: 

cd /work/terrain 

Next, start execution of the simulation by typing the command fogm. A 
welcome screen will appear on the display as shown in Figure 9.1. Pressing all 
three of the mouse buttons simultaneously will stop the program and return 
to the UNIX command level. This option of pressing all three buttons to exit is 
available at any time during the execution of the program. Pressing the middle 
mouse button advances the display to the next screen of instructions. When the 
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user has advanced through the welcome screen and the two instruction screens 
(Figures 9.2 and 9.3) he is presented with a display showing a two-dimensional 


contour map. This is the prelaunch phase of the simulation. 


C. PRELAUNCH CONTROLS 
The purpose of the prelaunch phase is to allow the user to designate a missile 
launch position and a suspected target location position. In effect, the user 
describes an initial flight path for the missile. 
1. The Prelaunch Display 
The prelaunch display is divided into three sections as shown in Figure 
9.4. The upper right corner of the dispiay contains an instruction box which 
summarizes the functions of the mouse buttons for this phase. The lower right 
corner contains a prelaunch statistics box. The meanings of the various items 
within the statistics box are explained below. The majority of the display is 
occupied by a two-dimensional contour map. Each of the square grids on the 
contour map represents a one square kilometer area. The colors on the map can 
be interpreted as follows. Green areas indicate terrain that is covered with 
vegetation that is greater than one meter high. Brown areas indicate terrain 
where the vegetation is less than one meter high. Within each of the color 
categories, the elevation of the terrain is indicated by the intensity of the color, 


with the brighter colors representing the higher elevations. 
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2. Selecting the Launch Position 

The launch position must be selected first. To select the launch position, 
use the mouse to move the red arrow cursor to the desired location on the contour 
map. As the cursor is moved, the UTM coordinates of the current cursor location 
are shown in the Launch Position field of the statistics box. These coordinates 
can be used when a more accurate selection of the launch position is required than 
is obtainable from the contour map alone. When the cursor is in the desired 
position, press the left mouse button to lock in that position. A blue circle will 
appear on the contour map showing the position selected and the workstation will 


> confirming the selection. The launch position can be changed any time 


“beep,” 
before the launching of the missile by simply moving to the new desired location 
and pressing the left mouse button. 
3. Selecting the Target Position 
The target position can only be selected after a launch position has been 


set. After the launch position has been selected, moving the cursor over the 


contour map produces the following effects: 


- The UTM coordinates of the current cursor position are shown in the Target 
Location field of the statistics box. 


- A “rubber band” line is drawn on the contour map from the launch position 
to the current cursor location. This line represents the flight path the missile 
would take if the current cursor position was selected as the target location. 


- The direction and length of the flight path represented by the above line are 
displayed in the statistics box in the Heading and Distance fields respectively. 


Once the cursor is at the desired target location, press the right mouse button 
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to lock in the position. A red circle will appear on the contour map showing the 
selected location and the workstation will ‘“‘beep,”’ confirming the selection. 

The missile is now ready for launch. The target location can be changed 
any time before launch by simply moving the cursor to the desired new location 
and pressing the right mouse button. 

4. Launching the Missile 

Launching can not take place until both a launch and target location 
have been selected. If the launch and target locations selected are acceptable, the 
missile is ‘‘launched”’ by pressing the middle mouse button. 

If this is the initial launch of this execution of the program, a several 
(three to four) minute delay will follow during which calculations are done to 
construct the upcoming three-dimensional scenes. Again, this delay only occurs 
during the first launch of any execution. Subsequent launches proceed with no 
delay. During this delay, a countdown will appear in the bottom of the statistics 


box. Launch occurs when the countdown reaches zero. 


D. IN-FLIGHT CONTROLS 
1. The In-Flight Display 
After the missile is launched. the display changes to the in-Hight display 


shown in Figure 9.5. The left side of the display contains: 


- A three-dimensional view of the terrain as seen from the missile’s camera. 


- Aslider bar scale along the bottom edge indicating the camera pan angle. 
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Figure 9.5 The In-Flight Display 


- A slider bar scale along the left hand edge indicating the camera tilt angle. 


- A box in the lower left corner containing either the word DESIGNATE or 
REJECT. The word DESIGNATE in this box indicates that the missile is 
not locked on to a target and is waiting for a command to designate one. 
The word REJECT indicates that the missile is locked on to a target and is 
waiting for a command to reject that target. 


- Cross hairs used to sight the camera onto a target. 
The upper right corner of the display contains a scaled copy of the contour map 
seen in the prelaunch phase. The red arrow superimposed on the contour map 
shows the missile’s current position (the tail of the arrow) and its direction of 
flight. The red rectangle on the map indicates that area of the terrain that is 

currently being shown in the three-dimensional display. 
The middle right section of the display contains four indicators which 


show the following: 


The speed of the missile in knots. 


The direction the missile is traveling in degrees. 
The height of the missile above ground level (AGL) in feet. 
The height of the missile above mean sea level (MSL) in feet. 


A slider bar indicating the zoom setting of the camera in degrees. 

The lower right section of the display contains a summary of the functions 
performed by the mouse and dials. These are explained further below. The in- 
flight phase continues until the missile impacts a designated target or all three 


mouse buttons are pressed simultaneously (to stop the execution of the 


simulation). 
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2. Controlling the Camera 
The ranges and initial values of the camera’s functions are shown in 


Table 9.1. All of the camera’s functions are controlled with the mouse. 


- To pan the camera, move the mouse left or right as needed. 
- To tilt the camera, move the mouse up or down as needed. 
- To zoom in to a tighter field of view, press the left mouse button. 


- To zoom out to a wider field of view, press the right mouse button. 
3. Controlling the Missile Flight 
The missile can be controlled by changing its direction, speed, and 
altitude. The ranges and initial values of each of the flight parameters is shown 
in Table 9.2. The missile flight parameters are controlled by using the dials on 
the IRIS’s button/dial box (see Figure 9.6). Dial zero (lower left) controls the 
missile’s direction, dial one (lower right) controls the missile’s altitude, and dial 


two (above dial zero) controls the missile’s speed. Refer to the display’s control 


TABLE 9.1 CAMERA CONTROL RANGES AND INITIAL VALUES 


Control =_ Initial Value 


25 degrees right 25 degrees left | O degrees 


25 degrees down | 15 degrees up 15 degrees down 


| | 55 degrees | 8 degrees | 95 degrees | 





TABLE 9.2 MISSILE CONTROL RANGES AND INITIAL VALUES 


Control sin Initial Value 


Altitude 10,000 MSL 200 AGL 200 AGL 
Speed 400 kts O kts 200 kts 
Direction | 359.9 degrees | O degrees From prelaunch 
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Figure 9.6 IRIS Dial Box Fuctions 
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summary box for a reminder of each dial’s purpose and location during flight. 


The controls are used as follows: 


- Direction of flight - Turning dial zero clockwise turns the missile to the 
right. Turning it counterclockwise turns it to the left. The missile will move 
freely through the 360 degree mark so that, for example, turning the missile 
right two degrees from a heading of 359 degrees will produce a heading of 
001. 


- Altitude - Turning dial one clockwise increases the missile’s altitude up to 
the maximum of 10,000 feet MSL. Turning the dial counterclockwise 
decreases the missile’s altitude. The simulator will not allow an altitude to 
be selected that is less than 200 feet above ground level. 


- Speed - Turning dial two clockwise increases the missile’s speed, while 
counterclockwise decreases the speed. 


4, Designating and Rejecting Targets 

The middle mouse button is used to designate (lock on to) and reject 
(release the lock on) targets. When the missile is not locked on to a target the 
word DESIGNATE will appear in the lower left corner of the display. To 
designate a target, center the target within the cross hairs and press the middle 
mouse button. In order for the missile to lock on, some portion of the target 
must be in the center of the cross hairs. If the designation is successful, the 
workstation will ‘beep’? and word REJECT will appear in place of the word 
DESIGNATE on the display. Once a target is designated the missile will 
automatically adjust its heading and altitude to home in on the selected target. 
An explosion is displayed after impact with the target occurs. The user is then 


returned to the prelaunch phase of the simulation to begin another launch. 
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A locked on target can be rejected and missile flight control returned to 
the user by pressing the middle mouse button any time before impact with the 
target occurs. The workstation will respond with a “beep” and the 
reject/designate box will again show the word DESIGNATE. The missile is now 


ready to accept the designation of a new target. 
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X. CONCLUSIONS AND RECOMMENDATIONS 


A. LIMITATIONS 

There are several limitations to the flight simulator presented in this study. 
First, a trade-off had to be made between resolution and frame update (display) 
speed. Even though data was available at a resolution of twelve and one-half 
meters, the simulator uses one hundred meter resolution in order to achieve an 
acceptable frame update rate. 

Second, the simulator’s flight is confined to a ten kilometer square area. Any 
ten kilometer square area of the DTED file can be used during a run of the 
simulation, but the simulator must be exited before switching to a new area. This 
limitation is not too restrictive for the current range of the FOG-M, but may be 
inadequate if the range of the missile is increased as planned. 

Third, road data is available in a format usable by the simulator for only one 
10 kilometer square area. Since access routines were not developed for the DFAD 
data file, roads must be digitized by hand. 

Fourth, the simulator does not model any of the missile’s flight dynamics. As 
stated earlier, this limitation was imposed only because of development time 
constraints. It is felt that the dynamics can be acceptably modeled without 


adversely affecting the performance. 
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B. FUTURE RESEARCH AREAS 

A follow-on to this project, which will provide more realistic targets and 
allow viewing of the scene as seen from inside any of them, is currently underway 
at the Naval Postgraduate School. The project's plans are to use the Ethernet to 
allow several workstations to take part in the simulation simultaneously. Each 
workstation will control one weapon (a target or the missile) and its monitor will 
display the scene as viewed from that weapon. 

Work is also underway at the Naval Postgraduate School in the use of 
digitized photographic images on the IRIS. This work could possibly be 
incorporated into the FOG-M project through the use of digitized target images, 
digitized cultural features, or digitized textures for the terrains. 

Another possible research area is the addition of various environmental effects 
into the simulation. These include clouds, smoke, and rain, which affect the 
camera’s view by reducing visibility, and also dust, which aids the missile 
operator in acquiring moving targets. 

Much work could be done in the area of the missile’s flight dynamics. The 
goal would be to provide an acceptably accurate model without too much of a 


sacrifice in speed. 


C. SUMMARY AND CONCLUSION>S 
The project has proven the practicality and feasibility of building a low-cost 


flight simulator with commercial, off-the-shelf hardware. With a relatively small 
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investment of time and funds, a simulator with significant capabilities was 
developed. As the speed and power of graphics hardware increases, even more 


realistic displays at faster update rates will be possible. 
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APPENDIX A — MODULE DESCRIPTIONS 


BUILD ROAD.C 
Input: None. 
Output: None. 


Side Effects: Modifies the global array road, an array of graphical objects, where 
each object contains the polygons representing the road in a 
particular gridsquare. 


Description: Build road reads the file road width and centerline information 
from the file Road.data and constructs polygons which represent 
the road. The polygons are stored in the array of graphical objects 
road. A more detailed discussion of building the roads is contained 
in Chapter VIII. 


BUILDTERRAIN.C 

Input: None. 

Output: None. 

Side Effects: Buildterrain modifies the global arrays savetriangle and gridcolor. 


Description: Buildterrain reads terrain height information from the global array 
gridpizel and constructs the terrain as a set of planar triangles. 
The details of constructing the triangles and the format of the 
savetriangle and gridcolor arrays can be found in Chapter VI. 


COLORRAMP.C 


Input: The inputs to colorramp are two booleans, greyscale and init. If 
greyscale is TRUE, the terrain, sky, and target colortable entries 
are defined in shades of grey to produce a black-and-white image. 
If greyscale is FALSE, the terrain colors are green, the sky is blue, 
and targets are brown. Init is set to TRUE when this routine is 
initially called, so that every entry in the colortable is defined, 
including those for terrain, sky, targets, and writemasked lines on 
top of the contour maps. Should the display be switched between 
color and black-and-white, only the terrain, sky, and target entries 
need to be redefined, which is what happens when init is FALSE. 


Output: None. 


Side Effects: Colorramp changes the system’s colortable, and thus determines 
the colors that appear on the display for the images drawn by 
other routines. 
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Description: 


Colorramp is called by the main program fogm as part of the 
initialization that takes place before the flying loop is entered. At 
that point, greyscale is set to its default value (usually FALSE, 
indicating color images) and init is TRUE. The readcontrols 
routine also calls colorramp to toggle the display image between 
color and black-and-white, based on the position of one of the 
dials. This call is made with the desired value for greyscale and 
with init FALSE. Colorramp uses the IRIS routine mapcolor to 
directly update the colortable for the contour map colors, and calls 
the user written routine gammaramp to define appropriately 
shaded ranges of the greens and browns (or greys) used for the 
terrain and targets. 


COMPASS.C 


Input: 


Output: 


Side Effects: 


Description: 


Compass takes as input a float, direction, which is an angle in 
radians. 


Compass returns a float which is the compass direction in degrees 
corresponding to the input direction. 


None. 


The function Compass converts an radian angle measured using 
the standard mathematical convention, and converts it to a degree 
angle measured using the standard navigational convention. 


DISP TERRAIN.C 


Input: 


Output: 
Side Effects: 


Description: 


Display terrain takes eleven inputs: the X, Y, and Z, coordinates 
of the missile position VX, VY, and VZ; the X, Y, and Z 
coordinates of the camera’s look-at position PX, PY, PZ; the field 
of view angle (camera zoom value), FOVY; and the X and Z 
ranges of gridsquares to be displayed, FIRST X, FIRST Z, 
LAST X and LAST Z. 7 7 


None. 
None. 


Disp terrain outputs a frame of the terrain scene to the monitor 
using a hidden surface algorithm. The scene contains terrain, 
roads, and targets. Details of the hidden surface algorithm can be 
found in Chapter VI. 
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DIST TO LOS.C 


Input: 


Output: 


Side Effects: 


Description: 


Dist to los takes seven inputs: the X, Y, and Z coordinates of the 
start of a line segment; the X, Y, and Z coordinates of the end 
point of a line segment; and three dimensional array, pt, which 
contains the coordinates of a point. 


Dist to los returns a float which is the perpendicular distance 
from the input point, pt, to the input line. 


None. 


Function which computes the perpendicular distance from a point 
to a line in three-space. 


DO BOUNDARY.C 


Input: 


Do boundary takes the following inputs: 
- An integer Bound type which is interpreted as: 

0 - a diagonal boundary 

1 - a horizontal boundary 

2 - a vertical boundary 
- An integer which triangle that is interpreted as: 

O - the lower triangle of the gridsquare. 

1 - the upper triangle of the gridsquare. 
- The indices, zgrid and zgrid, of the gridsquare for which the road 
is being constructed. 
- The coordinates of the start point of the boundary stored in a 
three dimensional array, bound start. 
- The coordinates of the end point of the boundary stored in a 
three dimensional array, bound end. 
- The coordinates of the start point of the left side of the road 
stored in a three dimensional array, left start. 
- The coordinates of the end point of left side of the road stored in 
a three dimensional array, left end. 
- The coordinates of the start point of the right side of the road 
stored in a three dimensional array, right start. 
- The coordinates of the end point of right side of the road stored 
in a three dimensional array, right end. 
- A boolean, start corner flag, which is TRUE if the gridsquare 
corner at the boundary’s start is ALREADY in the road polygon 
array, FALSE otherwise. 
- A boolean, end corner flag, which is TRUE if the gridsquare 
corner at the boundary’s end is ALREADY in the road polygon 
array, FALSE otherwise. 
- The partially complete road polygon array, road poly. 
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Output: 


Side Effects: 


Description: 


- An integer, vertex cnt, that is the number of vertices currently in 
the road poly array. 


Do boundary outputs the following: 

- start corner flag (see Inputs for a description) 

- end corner flag (see Inputs for a description) 

- road poly, the road polygon array with the vertices along this 
boundary added. 


- vertex cnt (see Inputs for a description) 
None. 


Do boundary’s purpose is to find all the intersections of the road’s 
left and right sides with the input boundary of a gridtriangle. As 
an intersection is found the point is put into a temporary array. 
After all the intersections are found for the boundary the points in 
the temporary array are sorted then added to the existing 
road poly array. The order of the sorting is such that the resulting 
road poly array will be ordered counterclockwise. See Chapter 
VIII for a detailed description of building the roads. 


EDIT INDBOX.C 


Input: 


Output: 
Side Effects: 


Description: 


The inputs to edit indboz are the name of the indicator object, the 
tags within that object for each of the indicators, and current 
values for the following missile parameters: X,Y, and Z position 
coordinates, pan, tilt, and zoom angles, and designate/reject 
status. 


None. 


Since edit indbor changes the indicator object, it has the side 
effect of changing the display when the indicator object is next 
called and displayed. 


The indicator object is edited between each display frame so that 
the heads-up display and the indicator box indicators show the 
current values for the missile’s speed, heading, altitude, camera 
pan angle, camera tilt angle, camera field of view (zoom), and 
designate/reject status. The input speed, heading, and MSL 
altitude (Y position coordinate) are converted to strings for 
display. AGL altitude is calculated as the difference between MSL 
altitude and the elevation of the ground directly below the missile 
as obtained from gnd level with the X and Z position coordinates 
as input. The boolean designate determines whether 
“DESIGNATE” or “REJECT” is printed in the lower left corner 
of the terrain display. Finally, the positions of the tilt, pan, and 
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zoom indicators are calculated from the missile parameters. The 
equations in the code have been simplified to avoid excess 
computation; the derivations are given below. 


The x screen coordinate of the zoom (field of view, or fov) indicator 
is fixed. The y screen coordinate varies from 200 (at 8 fov) to 70 
(at 55° fov). The input missile parameter zoom is in tenths of 
degrees, and thus ranges from 80 to 550. The y coordinate is 
determined from Equation A.1. 


zoom 
y = 200 — =: 
10 








eae 
ee 


oo — 8 
(A.1) 


zoom * —0.2766 + 222.128 


Likewise, the screen x coordinate of the tilt indicator is fixed, while 
the y coordinate varies from 680 (at +25° tilt) to 50 (at —25° tilt). 
The input missile parameter ¢ilt is in radians, and is converted to 
degrees by multiplying it with the RTOD (Radians TO Degrees) 
constant from the header file fogm.h. The y coordinate of the tilt 
indicator is calculated as shown in Equation A.2. 


ms 
| 


680 — 50 
50 {ft *DTOR) + 25 | ‘i a 


29 ~— —209 
(A.2) 


tilt * 721.92682 + 365 


The pan slider bar is horizontal, so the y coordinate is fixed, and 
the x coordinate ranges from 120 (at —25° pan) to 750 (at +25 
pan). Like tilt, the pan value is in radians and must be converted 
to degrees. The pan indicator x coordinate is given by Equation 


A.3. 
750 — 120 
i (5) | (pan * DTOR) + 25 | : as 


29 — —20 
(A.3) 


pan * —721.92682 + 435 
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EXPLOSION.C 


Input: 
Output: 


Side Effects: 


Description: 


FOGM.C 
Input: 


Output: 


Side Effects: 


Description: 


None. 
None. 
None. 


The explosion routine simulates the effect of a missile destroying a 
target by rapidly flashing a succession of red, black, and yellow 
screens. One buffer is kept black to pronounce the flash effect, and 
the other buffer is alternately cleared to red, yellow, red, yellow, 
and red. A short pause with a cleared, black screen is provided 
before the routine exits. 


Fogm is the name given to the main program in the simulator. It 
has no parameters, but gets data from its header files and through 
the readdata routine. Interactive input is also received vial the 
readcontrols routine. 


None. 
None. 


The fogm program consists of global variable declarations, local 
variable declarations, system initializations, an active loop, and 
some exit housekeeping. The initialization portion includes reading 
in the DMA elevation data, making network connection (if in use), 
setting the IRIS display configuration, defining the color table 
entries, building all of the graphical objects used in the displays, 
and computing the lighting and position of the polygons used to 
produce the terrain image. Within the active loop is some 
additional initializations and the flying loop. In the active loop 
initializations, the dial and mouse controls are reset to their initial 
defaults, and the display buffers are loaded with the images that 
remain unchanged during flight simulation (the contour map and 
the legend/instruction box). Control is then passed to the flying 
loop, which produces the flight simulation images until either a 
target is hit or the simulation exit command is received. If a target 
was hit, an explosion is displayed and the pre-launch phase of 
designating launch and target positions is re-entered. If all three 
mouse buttons have been pressed, the display is cleared and 
various system parameters are reset to provide a graceful exit from 
the simulator. 
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The flying loop contains the subroutine calls that produce the 
simulation of flight. First, the mouse and dials are checked for 
control input. Then the targets’, missile’s, and lookat reference 
point’s positions are all updated based on the elapsed time since 
the previous frame and the appropriate speeds. View bounds is 
called to determine which one kilometer grid squares are in view, 
and then the indicators are all updated to show the new control 
values, missile statistics, and view area. The main display routine 
then draws the appropriate sections of the terrain, plus cultural 
features and targets where appropriate. Finally, the updated 
indicator objects are drawn, and the display buffers are swapped to 
display the newly created image. 


GAMMARAMP.C 


Input: 


Output: 
Side Effects: 


Description: 


The inputs to gammaramp are a correction factor, a color table 
starting index, the number of color table entries (shades) to be 
defined, red, green, and blue intensities for the brightest color to be 
defined, and finally, red, green, and blue intensities for the darkest 
color to be defined. 


None. 


Gammaramp has the side effect of defining entries in the system 
color table. 


Displayed colors do not correspond linearly to the numeric red, 
green, and blue intensity values that are used to produce them. If a 
range of colors (0 .. #colors-1) is defined in the straightforward way 
with a uniform increment, the intensity of the n~ color (J) is 
given by Equation A.4, and the bright colors will appear more 
widely spaced than the dark colors. 

MazI — Minl 
=f * Sou (A.4) 


; # colors 


Gammaramp avoids this by using a power function to increase 
spacing between the dark colors’ intensity values and to decrease 
the intensity increment as the colors get brighter. The strength of 
the correction is determined by a value 7, which is constant for a 
given range, but must be experimentally determined for each range 
that differs in color or number of colors. FOG-M uses a ¥y value of 
1.5. The intensity of the n™ color in a gammaramp created table is 
given by Equation A.5. 
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l 


n 7 
a * (MazI — MinI) + MinI (A.5) 
#colors —1 


GET TGT POS.C 


Input: 


Output: 


Side Effects: 


Description: 


The input to get tgt pos is a socket number for Ethernet 
communication (if in use), a boolean indicating designate/reject 
status, the index of the currently designated target, and the 
‘‘name”’ of the tank object. 


Output is the new X,Y,Z position coordinates of the currently 
designated target. 


Get tgt pos updates several global data structures. It sets the 
number of target images, updates the target position arrays, and 
updates the array of target object names. 


The primary purpose of get tgt pos is to move the targets in the 
simulation. If the networking capability is in use, the target 
positions for the next frame are received over the network. When 
networking is not in use, targets are moved at a set speed of fifteen 
knots, and reverse course when they reach the boundaries of the 
ten kilometer square terrain area. As explained in Chapter VII, an 
array of graphical objects is defined to match one object per one 
hundred meter square of terrain, and this array is also used as 
booleans to indicate the presence or absence of targets in the one 
hundred meter grid square. Get tgt pos begins by removing each 
target from this array. New target positions are calculated or 
received over the network. If one of the targets has been ‘‘locked- 
onto,’ its new position is returned to be used as the current aim 
point for the missile. This is easily determined if networking is off 
because the designated target’s index remains the same and the 
new position can be directly accessed. The index correspondence is 
not guaranteed when networking, so the index of the new target 
whose coordinates are closest to the old targeted point is used. 


Targets that straddle a one hundred meter grid square boundary 
must be drawn on top of both (or possibly all four) grid squares in 
order to avoid being partially obscured by whichever square is 
drawn last. (The target must be drawn immediately after the grid 
square on which it rests to ensure that the target will be obscured 
when it should be by terrain drawn in the foreground.) Since the 
calculation of boundary intersection requires several trigonometric 
functions plus an allowance for the distance between the center of 
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the tank and its boundaries (which varies with the direction of the 
tank), a simplifying algorithm is used. If the tank is close enough 
to a boundary that the most distant part of the tank might cross 
the boundary, the target is also drawn after the adjoining grid 
square(s) (see Figure 7.3). This is done by adding a ‘“‘new”’ target 
to the array of target objects. The “‘new”’ target object is drawn at 
the exact same location in the three-dimensional terrain, but it is 
drawn after a different one hundred meter grid square, so it will 
have different target object array indices, and be in a separate 
target object. 


After all of the targets (originals and boundary copies) have 
updated positions and target object array indices, objects are 
added to the target object array as described in Chapter VII. This 
array is then used by the terrain display routine to actually draw 
the targets. 


GND LEVEL.C 


Input: 


Output: 
Side Effects: 


Description: 


Gnd level takes as inputs the X¥ and Z coordinates of the point for 
which the elevation is desired. 


Gnd level returns a float which is the elevation at point X and Z. 


None. 


Gnd level computes, through interpolation, the scaled elevation of 
any point within the terrain boundaries. A calculation is done to 
determine which gridtriangle contains the point. Then, using the 
known elevations at the vertices of the triangle, the elevation of the 
point is found. 


IN THIS POLY.C 


Input: 


Output: 


In this poly takes the following inputs: 

- An array of points, polygon, which define a polygon. (Note: only 
the X and Z coordinates of the points are used, the Y value is 
ignored). 

- An integer, num vertez, that is the number of vertices in 
polygon. 

- A point, pnt, that is to be tested. (Note: only the X and Z 
coordinates of the point is used, the Y value is ignored). 


In this poly returns a boolean which is TRUE if pnt is inside the 
polygon defined by polygon, FALSE otherwise. 
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Side Effects: 


Description: 


None. 


In this poly is a function which tests whether a point is inside a 
given polygon, where both the point and the polygon are in the XZ 
plane. The algorithm used constructs a bounding box around the 
polygon. If the point lies outside the bounding it obviously can 
not be inside the polygon. If the point lies inside the bounding box 
a further test is made. A line is constructed from a point outside 
the bounding box to the point to be tested. Each of the edges of 
the polygons are then tested to see if they intersect the constructed 
line and a count is kept of the number that do intersect. The 
point lies inside the polygon if and only if the constructed line 
intersects an odd number of the polygon’s edges. 


INIT CTRLS.C 


Input: 


Output: 


Side Effects: 


Init_ctrils takes as inputs the initial altitude of the missile, in feet; 
the initial heading of the missile in degrees; and a _ boolean, 
greyscale, which is TRUE if greyscaled terrain is to be displayed 
and FALSE if color terrain is to be displayed. 


Init ctrls has as outputs the initial pan angle of the camera in 
radians; the initial tlt angle of the camera in radians, and the 
initial zoom setting of the camera in tenths of a degree. 


The MOUSEX, MOUSEY, DIALO, DIAL1, DIAL2, and DIAL3 


valuators are set as a result of calling this routine. 


Description: Init ctris’s purpose is to initialize the mouse and dial valuators 
used for the operator controls. The initial altitude, heading, and 
greyscale valuator settings are passed in as inputs. The pan, tilt, 
and field of view settings are read from an "include" file and their 
values passed back as outputs. 

INIT IRIS.C 

Input: None. 

Output: None. 

Side Effects: Calling this routine sets the Iris attributes and configures the Iris. 

Description: Init irtis accomplishes the following: it puts the Iris into 


doublebuffer mode, sets the chunksize (the minimum memory 
increment used in objects), sets the monitor type to either NTSC 
or HZ60, and enables backface polygon removal. 
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INIT TGTS.C 


Input: 
Output: 
Side Effects: 


None. 
None. 


Init _tgts always initializes the global target object array to all 
zeros. If target data is not being received over the network, 
intt_tgts also defines ten targets by setting initial values in the 
global target counter, target position array, and target direction 
array. An auxiliary function init tgt is used to perform the actual 
update of the global arrays. 7 


INTERP ELEV.C 


Input: 


Output: 


Side Effects: 


Description: 


Interp elev takes three inputs, each an array of X, Y, and Z 
coordinates, representing a point. One array is the start point of a 
line, the second array is the end point of a line, and the third array 
is a point along the line. 


Interp elev returns a float that is the elevation value of the point 
along the line. 


None. 


Interp elev returns a float which is the linear interpolation of the 
Y (elevation) coordinate of the point along the line, based on the 
elevations at the start and end points of the line. 


LIGHT ORIENT.C 


input: 


Output: 


Side Effects: 


Description: 


Light ortent takes as inputs the following: 

- An array of coordinates for the polygon. 

- An integer, num coords, the number of coordinates in the 
polygon. 

- The X, Y, and Z coordinates of a point that is "behind" the 
polygon (an interior point). 

- The X, Y, and Z coordinates of a light source. 

- The minimum and maximum color map indices to be used for 
this polygon. 


Light ortent returns the color map index of the color to use in 
lighting this polygon. It also reorders the polygon array (if 
necessary) so that the points are ordered counterclockwise. 


None. 


Light orient computes a lighting for a polygon based on Lambert’s 
cosine law, which states that the intensity of the light reflected 
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from an object is proportional to the cos(®), where ® is the angle 
of incidence of the light ray. (see Figure 5.2). Light orient also 
orders the vertices of the polygon in a counterclockwise fashion so 
that backface polygon removal can take place (see the module 
description for npoly orient). 


LINE INTER2.C 


Input: 


Output: 


Side Effects: 


Description: 


Line inter2 takes the following inputs: 

- An array containing the X and Z coordinates of the start point of 
line one is ignored.) 

- An array containing the X and Z coordinates of the end of 
line one. (Note: a three element array is used, but the second, Y 
coordinate, element is ignored.) 

- An array containing the X, Y, and Z coordinates of the start of 
line two. (Note: a three element array is used, but the second, Y 
coordinate, element is ignored.) 

- An array containing the X, Y, and Z coordinates of the end of 
line two. (Note: a three element array is used, but the second, Y 
coordinate, element is ignored.) 


Line tnter2 returns as outputs: 
- An array containing the X and Z coordinates of the intersection 
of line one and line two. If the lines do not intersect these values 
are undefined not considered in the calculation). 
- An integer which can be interpreted as follows: 
O - the lines do not intersect. 
1 - the lines intersect, but the intersection uses an 
extension of at least one of the lines past its start or 
end points. 
2 - the lines intersect, and the intersection occurs 
between the input start and end points of both lines. 


None. 


Line tnter2 computes the point of intersection between two lines 
in the XZ plane. The type of intersection, as explained above in 
"Output" is also determined. Throughout the routine, three 
element arrays are used for compatibility with other routines. The 
second, Y, coordinate is not considered in any of the calculations. 


MAKEINDBOX.C 


Input: 


None. 
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Omtput: 


Side Effects: 


Description: 


Makeindbor returns a graphical object “‘name,”’ tags for editing the 
speed, direction, altitude, and designate/reject readouts, and tags 
for editing the zoom, pan, and tilt indicators. 


None. 


Makeindbor generates a graphical object that contains both the 
indicator box in the middle of the displays on the right side of the 
screen and the “heads-up” display that is superimposed on the 
terrain image (Figure 6.8). The object consists almost entirely of 
straightforward line and character string drawing commands, but 
there are two interesting points. First, within a single object, there 
are two different coordinate systems: one for the indicators 
superimposed on the terrain, and another for the separate indicator 
box. This is accomplished with an ortho2 call for each coordinate 
system, and by bracketing each ortho2 with pushmatriz and 
popmatriz commands. Note that the heads-up display is truly 
superimposed; it is specified in two-dimensional screen coordinates 
as opposed to the three-dimensional terrain coordinates. 

The second interesting aspect is the movement of the slider bar 
indicators. Drawing the indicators as polygons would require a 
sequence of pushmatriz, translate, and popmatriz calls for each 
indicator, with movement achieved by editing the translate call. To 
avoid all of this matrix movement and multiplication, the 
‘‘triangle’’ of the indicator is actually an overlapped line that 
‘fills’? the triangle by spiraling inwards. The line is drawn relative 
to the indicated point. with each segment of the line specified as 
offsets from that initial point, rather than as absolute coordinates 
(Figure A.1). Movement of an indicator triangle defined in this 
way is achieved by editing the parameters of a move2 call in the 
object, which sets the current graphics drawing position to the 
indicated point on the slider bar scale. Makeindbor is called once 
by fogm before the flying loop is entered, and then the object is 
edited (to update the indicator values) and called (to display it) 
every frame. 


MAKEINSTRBOX.C 


Input: 
Output: 
Side Effects: 


None. 
Makeinstrboz returns the name of an object to fogm. 


None. 
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Indicator Fill using Line Segments 
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Description: 


Makeinstrbor creates the object that produces the display in the 
lower right of the screen (Figure 6.8) during flight simulation. This 
display contains the legend for the FOG-M controls and the flight 
parameters they affect. Makeinstrboz is called once by fogm to 
create the object, and then the object is called twice per flight to 
put the image into each buffer. Note that writemasks are not 
necessary as they are with makemap and makenavbor, because 
nothing else writes to the instruction box portion of the screen 
during flight. The image thus remains undisturbed in the bitplanes 
despite the changes in other screen areas. 


MAKEMAP.C 


Input: 


Output: 


Side Effects: 


Description: 


The input to makemap is the globally defined array of elevation 
and vegetation values, gridpirel. 


’ 


The output from makemap is a graphical object ‘‘name,’’ which is 


returned to fogm. 
None. 


Makemap generates the object containing the contour map and 
grid that appear full screen during the pre-launch phase, and 
appear in the upper right corner of the screen during flight 
simulation (Figure 4.1). The map is produced using the 
methodology described in Chapter IV. Fogm calls the object 
returned by drawcontour twice, in order to place the map image in 
both buffers. The image is then protected from overwrite by a 
writemask. Fogm also passes the object name to prelaunch, which 
uses it in much the same way as fogm. 


MAKESCREENS.C 


Input: 
Output: 


Side Effects: 


Description: 


None. 


Makescreens returns an array of objects: instruction panel, 
statistics box, flight path between launch and target endpoints, 
and the three welcome screens, plus tags to update the statistics 
and flight path. 


None. 


Makescreens builds all of the objects (mostly screens of text) that 
are used by prelaunch. 
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MAKETANK.C 


Input: 
Output: 


Side Effects: 


Description: 


None. 


Maketank returns the name of an object containing a single tank, 
drawn around the origin. 


None. 


Maketank builds a object that consists solely of the drawing 
commands to produce a single tank. The tank is thirty-two feet 
long, ten feet high, and ten feet wide. Its center bottom is at the 
origin (coordinates 0,0,0), with its left side on the plane Z = -5, its 
back on the plane X = -15, its bottom on the plane Y = 0, and it 
faces to the right along the positive x axis. For each of the twenty 
polygon faces that make the tank, the X,Y, and Z coordinates of 
each polygon vertex are stored in an array, passed to lightorient, 
and then drawn with polf, the filled polygon drawing command. 
Lightorient ensures the vertices are ordered counter-clockwise in 
the array (with respect to an interior point) for backface polygon 
removal, and then calculates the appropriate color for the polygon 
using the same lighting model that is used for the terrain (see 


Chapter V). 


NEAREST TGT.C 


Input: 


Output: 


Side Effects: 


Description: 


Nearest tgt takes as inputs the X, Y, and Z coordinates of the 
missile position, and the X, Y, and Z coordinates of the camera’s 
look-at position. (The end points of the line of sight vector). 


Nearest tgt returns as output an integer, tgt tdz, which is the 
target index of the target that is closest to the line of sight vector. 


None. 


For each of the existing targets, nearest tgt computes the distance 
between the target and the line of sight vector. It returns the 
index of the target that was found to be closest. In the case of two 
targets which are the same distance apart, the highest index value 
will be returned. 


NPOLY ORIENT.C 


Input: 


Npoly ortent takes as input: 

- An integer, num coords, that is the number of vertices in the 
polygon. 

- An array containing the coordinates of the polygon. 

- The X, Y, and Z coordinates of a point that is "behind" the 


121 


Output: 


Side Effects: 


Description: 


polygon (an "interior" point). 


Npoly ortent returns as output an integer which is interpreted as: 
1 - the vertices of the polygon are ordered clockwise. 
2 - the vertices of the polygon are_ ordered 
counterclockwise. 


None. 


Npoly orient determines if the polygon is ordered clockwise or 
counterclockwise by computing two points: one along the normal 
vector and the other, the same distance from the polygon, but 
along the vector in the direction opposite the normal. Next the 
distance between these points and the "interior" point is 
computed. If the "interior" point is closer to the point along the 
normal vector, the polygon is ordered clockwise, otherwise the 
polygon is ordered counterclockwise. 


PRELAUNCH.C 


Input: 


Output: 


Side Effects: 


Description: 


The input to prelaunch is two arrays. The first contains objects, 
and the second contains tags for editing those objects. 


Prelaunch returns the .,Y, and Z coordinates of the missile’s 
designated launch position, and the initial direction of flight for the 
missile. This direction is returned in both radians and compass 
degrees (Figure 7.1). 

None. 


Prelaunch first provides three screens of introductory information. 
Each screen is an object defined by makescreens. After those, the 
user is presented with a full screen contour map of the ten 
kilometer by ten kilometer area available for overflight. Mouse- 
selected points define the missile’s initial position and direction of 
flight, and are displayed on top of the map. The map is writemask 
protected, so it is only drawn twice (once for each buffer) even 
though the flight path is repeatedly drawn and erased on top of the 
map. The flight path is made to act like a rubber band between 
the launch and cursor positions by repeatedly editing of the 
positions in the object containing the flight path line drawing 
commands. Once the flight path is confirmed, the launch position 
and heading are returned to the fogm program. 
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RANDNUM.C 


Input: 
Output: 
Side Effects: 


Description: 


Randnum uses the global random number seed. 
Randnum returns a floating point random number. 


The global seed value used by randnum is updated during every 
invocation. 


Randnum is a linear congruential pseudo-random number 
generator. The algorithm is a modified version of the one given by 
Sedgewick |Ref. 13]. It uses a a special piecewise multiplication 
routine mult to preserve the low-order digits of the newly 
generated seed even in case of overflow. The value returned is the 
new seed, scaled to fall between zero and one, inclusive. The 
random numbers are used in fogm to vary the point on the tank 
that the missile aims for. This simulates the variance in impact 
point that results from the optical homing of the real missile. 


RANDSEED.C 


Input: 
Output: 
Side Effects: 


Description: 


Randseed takes a long integer as input. 
None. 
Randseed updates the global random number seed value. 


The pseudo-random number generator implemented in randnum 
always returns the same string of numbers when it starts with a 
given seed value. Randseed provides the means to change that 
initial seed value so that different program runs will have different 
strings of ‘‘random’”’ numbers. 


READCONTROLS.C 


Input: 


Output: 


Side Effects: 


The inputs to readcontrols are the global X,Y, and Z random 
offset values for the aim point on the target, the current 
designate/reject status, and the black-and-white versus color 
boolean greyscale. 


All of the user-commanded control values are output from 
readcontrols: missile speed, heading and altitude, camera pan, tilt, 
and zoom angles, plus designate/reject status, greyscale status. 
Readcontrols also returns values for the booleans that control the 
active and flying loops. 


When a target is first designated, readcontrols calls randnum and 
updates the global target aim offsets randz, randy, and randz. 
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Description: 


Readcontrols checks the status of all of the valuators that provide 
input to the FOG-M simulator, and performs scaling, units 
conversion, and immediate processing, as appropriate. It 
determines whether to accept or reject a ‘“‘designate’’ command, 
based on the color index of the pixel at the center of the screen. (If 
a tank is in the crosshairs, the color index will be from the tank’s 
color ramp, and a designate command will be accepted. Otherwise, 
a designate command will be ignored.) 


READDATA.C 


Input: 
Output: 


Side Effects: 


Description: 


None. 
None. 
Readdata fills the global array gridpizel. 


Readdata opens and reads the values from the terrain elevation 
data file and stores the values in the gridpizrel array. Note that the 
elevation data file is arranged in a format as discussed in Chapter 
II. The gridpirel array is arranged in straight rows and columns 
analogous to the geographic positions of the data. 


ROAD BOUNDS.C 


Input: 


Output: 


Side Effects: 


Description: 


Road bounds takes as input the following: 

- Three arrays (ptl, pt2 and pt3) containing the X and Z 
coordinates of three points along the centerline of the road. The 
line segment from ptl to pt2 defines the first segment of the road. 
The segment from pt2 to pt3 defines the next segment of the road. 
- A float, width, which is the width of the road in feet. 


Road bounds returns the following as outputs: - Four arrays 
(left | pel, right pti, left pt2, and right _pt2) which contain the X 
and Z coordinates of the first segment’s left and right sides. The 
left side runs from left pt1 to left pt2 and the right side runs from 
right ptl to right pt2. 

- Four integers, first zgrid, first zgrid, last _xgrid and last zgrid, 
which are the indices of the bounding box surrounding the first 
road segment (see Figure 8.2). 


None. 


Given three points along the center line of the road, and the road’s 
width, road bounds computes the start and end coordinates for the 
first segment’s left and right sides. The end coordinates are 
computed as the intersection of the first segment’s left (or right) 
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side with the second segment’s left (or right) side. This insures 
that adjoining segments will meet cleanly. The second function of 
road bounds is to compute a bounding box around the first road 
segment. This box is defined as the row indices of the northern 
and southern most gridsquares that the road segment intersects, 
and the column indices of the eastern and western most gridsquares 
that the road segment intersects (See Chapter VIII for a more 
detailed discussion). 


SORT ARRAY.C 


Input: 


Output: 


Side Effects: 


Description: 


Sort array takes as inputs: 

- An array of points, pnts. 

- An integer that is the number of entries in the pnts array. 

- A boolean, which is TRUE if the array should be sorted in 
descending order, FALSE if the array should be sorted in ascending 
order. 

- The index number of the coordinate that is the sort key: 0 for the 
X coordinate, 1 for the Y coordinate, and 2 for the Z coordinate. 


Sort array returns the array pnts with the points sorted according 
to the input parameters. 


None. 


Sort array performs a simple "bubble-sort" of the input points 
according to the input parameters. 


UP LOOK POS.C 


input: 


Output: 


Side Effects: 


Description: 


Up look pos takes the following as inputs: 

- The heading of the missile in radians. 

- The pan angle of the camera in radians. 

- The tilt angle of the camera in radians. 

- The X, Y, and Z coordinates of the missile’s position. 

- The X, Y, and Z coordinates of the locked-on target (if any). 

- A boolean which is TRUE if the missile is locked-on a target, 
FALSE otherwise. 


Up look pos returns as outputs the X, Y, and Z coordinates of the 
camera’s look-at position. 


None. 


Up look position computes a point along the camera’s line of 
sight. If the missile is locked on a target, the look-at position is the 
locked-on target’s position. Otherwise it is any point along the 
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camera’s line of sight. See Chapter VI and Figure 6.2 for a more 
detailed discussion. 


UP MSL POSIT.C 


Input: Up msl posit takes as inputs: 
- The heading of the missile in radians. 
- The speed of the missile in knots. 
- The X, Y, and Z coordinates of the missile’s position. 
- The X, Y, and Z coordinates of the locked-on target (if any). 
- A boolean which is TRUE if the missile is locked-on a target, 
FALSE otherwise. 


Output: Up msl posit returns as outputs: 
- The new heading of the missile in radians, if it was changed to 
track a locked-on target. 
- The new heading of the missile in degrees measured in the 
compass convention. 
- A boolean which is TRUE if the missile is still flying (has not hit 
a target), and FALSE if the missile has hit the target. 


Side Effects: None. 


Description: Up msl posit calculates a new missile position for the next frame. 
The new position is either based on the commanded direction, 
speed, and altitude (when the missile is NOT locked onto a target), 
or the commanded speed and the direction to the target (if the 
missile is locked onto a target). For a detailed discussion of the 
routine, see Chapter VI. 


VIEW BOUNDS.C 


Input: View bounds takes as inputs the X, Y, and Z coordinates of the 
missile’s position; the X, Y, and Z coordinates of the camera’s 
look-at position; and the field of view (zoom) value. 


Output: View bounds returns as outputs the row indices of the northern 
and southern most gridsquares to be drawn, and the column 
indices of the western and eastern most gridsquares to be drawn. 


Side Effects: None. 


Description: The purpose of view bounds is to construct a bounding box around 
the gridsquares which are to be drawn. The box is constructed by 
extending the line of sight vector down until it intersects the 
minimum elevation plane. The view bounds extends 20 
gridsquares north, south, east, and west of this intersection point. 
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If the missile’s position is not within the bounds, the bounds are 
extended to include the missile’s position. For a more detailed 
discussion, see Chapter VI and Figure 6.5 
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APPENDIX B - SOURCE LISTINGS 


BUILD ROAD 


#include "stdio.h" 
#include "fogm.h" 
#include "files.h" 
#include "gl.h" 
#include "math.h" 


#define X O 
#define Y 1 
#define Z 2 


#define DIAGONAL 0 
#define HORIZONTAL 1 
#define VERTICAL 2 


#define LOWER 0 
#define UPPER 1 


build road() 
{ 
extern Object road/|99]{99); 
extern short gridpixel| 100]| 100); 
FILE *fp, *fopen(); 
float road _ width; /* road width if feet * / 
int num pts; /* number of data points 
for the road seqment */ 
int segnum = QO; 
char temp|100 ; 
MG emt, 1s 
int vertex cnt, num duplicates; 
float gnd_level(); 7 
float elev; 
float pt1![3], pt2(3], pt3/3); 
float nw corner|3], ne corner|3], sw corner 3), se_corner(3]; 
float right pt1[3), right pt2(3); 
float left pt1.3], left pt2/3); 
float north bound, south bound, east bound, west bound; 
float delta x, delta 2; - 
float seg dir; 
int ne flag, nw flag, se flag, sw flag; 
int xgrid, zgrid; - 
int first xgrid, last xgrid, first zgrid, last zgrid; 
float poly 1,10!/3]; 


frontbuffer(TRUE); 
fp = fopen(ROAD FILE,"r'); 
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while (fscanf(fp, "%e", &road_ width) != EOF) { 
fscanf(fp, "%d", &num_ pts); 
fscanf(fp, "%e %e", &ptl/X], &pt1[Z]); 
fscanf(fp, "%e %e", &pt2 X|, &pt2/Z]); 


delta x = pt2[X] - pt1[X]; 
delta z= pt2[Z] - pt1[Z]; 
seg dir = atan2(delta z, delta x); 
left ote = ptl[X] + eee dir + HALFPI)*road_ width/2.0); 
right ptl[X] = pt1[X] + (cos(seg dir - HALFPI)*road_ width/2.0); 
left pt1[Z} = pt1[Z] + (sin(seg dir + HALFPI)*road_ width/2.0); 
right pt1[Z} = pt1[Z] + (sin(seg dir - HALFPI) * road _ width/2.0); 
for (cnt = 3; cnt <= num) pts; 1, +-+cnt)} { 

if (cnt <= num pts) { 

fscanf(fp, '"%e %e", & pt3[X], &pt3[Z)); 
} 


else { 
pt3|X| = pt2 ini; 
pt3/Z] = pt2[Z); 


/* print new road segment number on title screen */ 
segnum += 1; 
pushmatrix(); 
ortho2(0.0, 1023.0, 0.0, 767.0); 
view port (0,1023,0,767); 
sprintf(temp, "Building road segment: %d%", segnum); 
color(BLUE); 
rectf(780.0, 20.0, 1010.0, 30.0); 
color(CYAN)}; 
cmov2i(780, 20); 
charstr(temp); 
popmatrix(}; 
/* determine the boundaries of this road segment */ 
road bounds(pt1, pt2, pt3, road width, left ptl, right ptl, 
left pt2, right pt2, &first xgrid, - 
&first zgrid, &last xgrid, &last zgrid); 
for (xgrid = first xgrid; xgrid <= last xgrid; ++xgrid){ 
for (zgrid = first zgrid; zgrid <= last zgrid; ++zgrid){ 
ne flag = FALSE; 
nw flag = FALSE; 
sw flag = FALSE; 
se flag = FALSE; 
vertex cnt = -1; 
east bound = (float)(xgrid + 1) * FT 100M; 
west bound = (float)(xgrid) * FT 100M; 
north bound = (float)(zgrid + 1) * FT 100M; 
south bound = (float)(zgrid) * FT 100M; 


sw _corner|X] = west bound; 

sw _corner|Z] = south bound; 

elev = gridpixel|zgrid||xgrid| & elev mask; 
sw _corner|Y| = pow(elev, ALTSCALE); 
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se corner|X| = east bound; 

se corner|Z| = south bounel 

elev = gridpixel/zgrid/ xecden & elev_mask; 
se corner|Y] = pow(elev, ALTSCALE); 


nw corner|X| = west bound; 

nw corner|Z) = north bound; 

elev = grid pixellzgrid ail (een & elev mask; 
nw corner|Y| = pow(elev,ALTSCALE); 


ne corner|X| = east bound; 

ne corner|Z| = north bound; 

elev = gridpixel|zgrid+1||xgrid+1| & elev_mask; 
ne corner{Y] = pow(elev, ALTSCALE); 


/* determine points of intersection between the left and 
right sides of the road and the eastern grid boundary 
and add these points to the polygon vertex array */ 


do boundary(VERTICAL, UPPER, xgrid, zgrid, se corner, ne corner, 
left ptl, left pt2, right ptl, right pt2, &se flag, 
&ne flag, polyl, &vertex cnt); 


/* determine points of intersection between the left and 
right sides of the road and the northern grid boundary 
and insert these points into the polygon vertex array */ 
do boundary(HORIZONTAL, UPPER, xgrid, zgrid, ne corner, 
nw corner, left ptl, left pt2, right ptl, 
right pt2, &ne flag, &nw flag, polyl, &vertex cnt); 


/* determine points of intersection between the left and 
right sides of the road and the diagonal and 
insert these pointsinto the polygon vertex array */ 


do boundary(DIAGONAL, UPPER, xgrid, zgrid, nw corner, se corner, 
left ptl, left pt2, right ptl, right pt2, &nw flag, 
&se flag, polyl, &vertex cnt); 
/* remove duplicate entries from the polygon array */ 
num duplicates = 0; 
for (i = 1; i <= vertex _ent; ++i) { 
if ((poly1[i}[0] = Spor i-1][/O]) && 
(poly 1|i](2] == poly1{i-1][2})) { 
for (} = i; } < vertex cnt - num duplicates; ++)) { 

poly 1[j]/0] = poly 1|j+1]/0); 

poly 1{j][1] = poly 1[j+1][1]; 

poly 1[j][2] = poly1[j+1]|2]; 


num duplicates += 1; 
} 
} 


vertex cnt -= num duplicates; 
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if (vertex cnt > 0) { /* add polygon to grid object */ 
if (road|zgrid]|xgrid’ != 0) { 
editobj(road|zgrid||xgrid]); 


else { 
road|zgrid||xgrid| = genob)(); 
makeobj(road|zgrid||xgrid]); 


color(ROADGREY); 
polf(vertex cnt +1, &poly1/0]/0}); 
linewidth(3); 
poly(vertex cnt + 1, &poly1{0)|0]); 
closeobj(); 

} 

vertex cnt = -1; 

ne flag = FALSE; 

nw flag = FALSE; 

sw flag = FALSE, 

se flag = FALSE; 


/* determine points of intersection between the left and 
right sides of the road and the southern grid boundary 
and insert these points into the polgon vertex array */ 
do boundary(HORIZONTAL, LOWER, xgrid, zgrid, sw corner, 
se corner, left ptl, left pt2, rmght pt1, 
right _pt2, &sw flag, &se flag, polyl, &vertex cnt); 


/* determine points of intersection between the left and 
right sides of the road and the diagonal and 
add these points to the polygon vertex array */ 


do boundary(DIAGONAL, LOWER, xgrid, zgrid, se corner, nw corner, 
left ptl, left pt2, right ptl, right pt2, &se flag, 
&nw flag, polyl, &vertex cnt); 


/* determine points of intersection between the left and 
right sides of the road and the western grid bound 
and add these points to the polygon vertex array */ 


do boundary(VERTICAL, LOWER, xgrid, zgrid, nw corner, sw corner, 
left ptl, left pt2, right ptl, right pt2, &nw flag, 
&sw flag, polyl, &vertex cnt); 


* remove duplicate entries from the polygon array * 
p 


num duplicates = 0; 
for (i = 1;i <= vertex cnt; ++i) { 
if ((poly1|i][0] == poly1[i-1][0]) && 
(poly 1{i}[2] == poly |i-1)/2})) 4 
for (j} = i; ) < vertex cnt - num duplicates; ++)) { 
poly 1]}](0) = poly1[j+1)[0); 
poly1[j][1] = poly1[j+1](1]; 
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poly 1[j}/2] = poly 1[j+1][2]; 


num duplicates += 1; 
} 
} 


vertex cnt -= num _ duplicates; 
if (vertex cnt > 0) { /* add polygon to grid object */ 
if (road/zgrid||xgrid] != 0) { - 
editobj(road|zgrid||xgrid)); 


else { 
road|zgrid||xgrid) = genobj(); 
makeobj(road zgrid||xgrid]); 


} 

color(ROADGREY); 

polf(vertex cnt +1, &poly1[0]/0)); 
linewidth(3); 

poly(vertex cnt + 1, &poly1|0][0}); 
closeobj(); 


} 


right ptl|X] = right pt2[X]; 
right ptl Z] = right pt2[Z); 
left_pt1[X] = left _pt2|X|; 
left pt1{Z] = left pt2[Z]; 
pilin) pein 
ptl|Z) = pez Ai 
pt2|Xi = pts|X]; 
pt2\Z] = pt3/Z); 
} 

} 

fclose(fp); 

frontbuffer(FALSE); 
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BUILDTERRAIN 


/* buildterrain.c - this function builds objects representing 1km grid squares 
in 3-D, with each grid square generating 4 objects, identical except for 


order of drawing 7 
#include "gl.h" /* get the graphics defs a 
#include "device.h" /* get the graphics device defs */ 
#include "fogm.h"  /* default constants a 
#include "math.h"  /* math function declarations */ 
buildterrain() 


{ 


/* array of data points to build the terrain */ 
extern short gridpixel|100 [100]; 


extern float savetriangle|99)(99|(2]/3][3]; 
extern long gridcolor|99|/99]; 

extern Object target |99]/99]; 

extern float ground plane)4|/3]; 

extern long gnd_ plane color; 


float gnd plane ht; 


Coord triangle1/3||3', triangle2(3]|3}]; /* polygon coordinates */ 
short xgrid, zgrid; /* indexes into the grid object array */ 
short endrow, endcol; /* miscellaneous indexes etc */ 


int row, col; 
float ax,ay,az; /* interior point for use in the lightpoly function */ 
float Ix,ly,lz; /* position of light source in lightpoly function */ 


/* min and max colormap indexes for lighting the poly */ 
long colormin, colormax; 


/* color index to use returned by the lightpoly function * / 
long colortouse, color1, color2; 


char temp|50];  /* character string for countdown */ 
float x,y; 
float gammacorr; 


long rampamax, rampamin, rampbmax, rampbmin; 
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int startrow, startcol, coordidx, vertex; 


Ix = 500* FT 100M; = /* direction of light source */ 
ly = 100000 * FT_100M; 
be =e 


frontbuffer( TRUE); /* write to front buffer */ 


/* compute color for ground plane polygon */ 

gnd plane ht = pow((float)MIN, ALTSCALE); 

ground plane|0|/(0] = -NUMXGRIDS * FEETPERGRID; 
ground plane/0||1| = gnd plane ht; 

ground plane(/0|/2} = NUMZGRIDS * FEETPERGRID; 


ground plane|1|/0] = 2.0 * NUMXGRIDS * FEETPERGRID; 
ground plane i|/1] = gnd_ plane ht; 
ground plane1|[2} = NUMZGRIDS * FEETPERGRID; 


ground plane 2]/0] = 2.0 * NUMXGRIDS * FEETPERGRID; 
ground plane 2|/1] = gnd_ plane ht; 
ground plane 2/2] = -2.0 * NUMZGRIDS * FEETPERGRID; 


ground plane'3|\0}] = -NUMXGRIDS * FEETPERGRID; 
ground plane 3|/1] = gnd plane ht; 
ground plane 3||2] = -2.0 * NUMZGRIDS * FEETPERGRID; 


lightorient(ground plane,4,0.0,0.0,0.0,1x,ly,lz,256,461, &gnd plane color); 


/* compute coordinates and colors for triangles and store in global 
variable savetriangle for later display ri 


for (col = 0; col < 99; ++col) { 
/* print new countdown number on title screen */ 
pushmatrix(); 
ortho2(0.0, 1023.0, 0.0, 767.0); 
view port (0,1023,0,767); 
sprintf(temp, "Countdown to launch: %d%", 98 - col); 
color(BLUE); 
rectf(780.0, 15.0, 1010.0, 30.0): 
color(CYAN); 
cmov2i(788, 20); 
charstr(temp); 
popmatrix(); 


for (row = 0; row < 99; ++row) { 


/* choose which color ramp to use so that a checker board 
effect is acheived */ 
if ((row+col)%2) { 
colormin = 256; 
colormax = 461; 
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} 


else { 
colormin = 462; 
colormax = 667; 


} 


/* build the polygon */ 

triangle1/0]|2] = (float)row * (-41.01) * 8.0; 

triangle1 0||0] = (float)col * 41.01 * 8.0; 

triangle1|0][1] = pow( (float) (gridpixel{row||col|&elev_ mask) 
_ ALTSCALE): 


triangle1({1]/2| = (float)row * (-41.01) * 8.0; 

triangle1/1]/0) = (float)(col+1) * 41.01 * 8.0; 

triangle1/1|{1] = pow((float)(gridpixeljrow)|col+1]&elev mask) 
,ALTSCALE); 


triangle1|2|/2| = (float)(row+1) * (-41.01) * 8.0; 

triangle1|2||0) = (float)col * 41.01 * 8.0; 

triangle1|2]|1] = pow((float)(gridpixel|row+1]|[col]&elev_mask) 
‘ALTSCALE); 


/* copy common vertex values for opposing triangle of grid */ 
for (vertex = 1; vertex < 3; ++vertex) { 
triangle2|vertex|/0] = triangle1|vertex||0]; 
triangle2|vertex||/1| = triangle1[vertex||1]; 
triangle2|vertex||2| = triangle1|vertex||2|; 


} 


/* change corner coordinate to form opposing triangle of grid * / 

triangle2|0|/2) = (float)(row+1) * (-41.01) * 8.0; 

triangle2|0|/0) = (float)(col+1) * 41.01 * 8.0; 

triangle2(0||1; = pow((float)(gridpixel|row+1]|col+1)&elev_mask) 
, ALTSCALE) ; 


/* compute an interior point for trianglel */ 
ax = triangle1[0]/0| + 15.0; 

ay = -10.0; 

az = triangle1|0]|2] -15.0; 


/* light and orient triangle */ 
lightorient(triangle1,3,ax,ay,az,lx,ly,lz,colormin, colormax, &color1); 


/* compute interior point for triangle2 */ 
ax = triangle2/0]|0] - 15.0; 

ay = -10.0; 

az = triangle2|0||2| +15.0; 


/* compute the light for and orient triangle2 */ 
lightorient(triangle2,3,ax,ay,az,1x,ly,lz,colormin,colormax, &color2); 


/* compute average color for the square */ 
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colortouse = (colorl + color2) / 2; 


} 


/* save this triangles color and orientation */ 
for (vertex = 0; vertex < 3; ++vertex) 
for (coordidx = 0; coordidx < 3; ++coordidx) { 
savetriangle|row |[col]/0]|vertex]|coordidx| = 
trianglel|vertex|/coordidx|; 
savetriangle|row |(col]{1|[vertex|{coordidx] = 
triangle2|vertex||coordidx]; 


gridcolor|row||col] = colortouse; 


} 


frontbuffer(F ALSE); 
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COLORRAMP 


/* constructs the color ramps to be used for displaying the terrain. 
If greyscale is true, constructs greyscale ramps, else it 
constructs green ramps. me 


#include "fogm.h" /* fogm constants */ 


colorramp(greyscale,init) 
int greyscale, init; 


{ 


int 1; 


/* build two gamma corrected color ramps with slightly offset colors */ 

if (greyscale) { 
gammaramp(1.5,256,205,255,255,255,50,50,50); /* even terrain ramp */ 
gammaramp(1.5,462,205,245,245,245,40,40,40); /* odd terrain ramp */ 
gammaramp(1.5,668,180,235,235,235,30,30,30); /* tank ramp */ 
mapcolor(SK YBLUE, 230,230,230); /* sky color */ 
mapcolor(ROADGREY,35,35,35); 


else { 
gammaramp(1.5,256,205,0,255,0,0,50,0); /* even terrain ramp */ 
gammaramp(1.5,462,205,0,245,0,0,40,0); /* odd terrain ramp */ 
gammaramp(1.5,668,180,255,165,55,75,55,0);  /* tank ramp */ 
mapcolor(SK YBLUE, 200,200,255); /* sky color */ 
mapcolor(ROADGREY,35,35,35); 

} 

if (init) { 


mapcolor(16,0,70,0); /* set up colors for contour map * / 
mapcolor(17,0,80,0); 
mapcolor(18,0,90,0); 
mapcolor(19,0,100,0) 
mapcolor(20,0,110,0) 
mapcolor(21,0,120,0); 
mapcolor(22,0,130,0); 
) 
) 
) 


’ 


’ 


} 


mapcolor(23,0,140,0 
mapcolor(24,0,150,0 
mapcolor(25,0,165,0); 
mapcolor(26,0,180,0); 
mapcolor(27,0,190,0); 
mapcolor(28,0,210,0); 
mapcolor(29,0,225,0); 
mapcolor(30,0,240,0); 
mapcolor(31,0,255,0); 
mapcolor(32,75,55,0); 

( 

( 

( 


3 


mapcolor(33,95,60,0); 
mapcolor(34,115,70,0); 
mapcolor(35,125,76,0); 
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mapcolor(36,135,83,0); 

mapcolor(37,145,90,0); 

mapcolor(38,155,97,0); 

mapcolor(39,165,105,0); 

mapcolor(40,175,110,0); 

mapcolor(41,185,113,0); 

mapcolor(42,190,118,0); 

mapcolor(43,200,127,0); 

mapcolor(44,210,135,30); 

mapcolor(45,225,145,35); 

mapcolor(46,240, 155,45); 

mapcolor(47,255, 165,55); 

for (i=64; 1<128; 1++) mapcolor(i,0,0,255); 

for (i= 128; i<256; i++) mapcolor(i,255,0,0); 
mapcolor(851,0,150,0); /* set up colors for instruction box */ 
mapcolor(852,255, 165,55); 

mapcolor(853,95,60,0); 

mapcolor(854,0,0,0); /* color for indicator box background* / 
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COMPASS 


/* compute the compass heading in degrees of the input direction. */ 
#include "fogm.h" /* fogm constants */ 


float compass(direction) 
double direction; 


{ 


float compassdir; 


compassdir = RTOD * direction; 
if (compassdir <= 90.0) 

compassdir = 90.0 - compassdir; 
else 

compassdir = 450.0 - compassdir; 


return(compassdir); 
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DISPLAY TERRAIN 


/* Compute which polygons need to be drawn to display the terrain and 
output them in an order such that the polygons farthest from the viewer 
are drawn first and those closest are drawn last. 


Note: Eventhough this seems like a long routine, it is broken into 8 
independent cases based on the direction the camera is looking. 
If you understand one case the others are merely mirror images of the 
algorighm for other octants. a 


#include "fogm.h" 
#include 'math.h" 
#include "gl.h" 


display terrain(vx, vy, vz, px, py, pz, fovy, 
firstxgrid, firstzgrid, lastxgrid, lastzgrid) 


Coord vx, vy, VZ, PX, PY, PZ} 
int fovy; 
short firstxgrid, firstzgrid, lastxgrid, lastzgrid; 


extern float ground plane/4]|{3}; 

extern long gnd plane color; 

extern Object road 99] 99); 

extern Object target 99|/99); 

extern float savetriangle 99|/99]/2] 3]{3]; 
extern long gridcolor|99]|/99]; 


double lookdir; 

int threshold, count, startx, startz; 
short xgrid, zgrid; 

float tanval; 

float y; 


if (TV) viewport (0,474 ,0,474); 
else viewport (0,767,0,767); 
pushmatrix(); 


color(SK YBLUE); 


clear(); 


ortho2(0.0,1023.0,0.0,767.0); /* outline the screen */ 
color(BLACK); 

recti(0,0,1023,767); 

popmatrix(); 


pushmatrix(); 
perspective(fovy,1.0,0.0,19500.0); 
lookat(vx,vy,vz,px,py,pz,0.0); 
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/* determine the direction of the line of sight */ 
lookdir = (double)atan2((float)(vz - pz), (float)(-(vx - px))); 
if (lookdir < 0.0) lookdir += TWOPI; 


/* lay down the ground plane */ 
color(gnd_ plane color); 
polf(4, ground plane); 


/* put the grid objects through the geometry engine in an order 


based on the lookdir. */ 


if (lookdir > SEVEN QTR PI) 


{ 


/* 8th OCTANT */ 
threshold = (int)(tan(lookdir+HALFPI) + 0.5); 


startx = lastxgrid; 

startz = firstzgrid; 

while (startz <= lastzgrid) { 
zgrid = startz; 
xgrid = startx; 


while ((xgrid <= lastxgrid) && (zgrid <= lastzgrid)) { 


color(gridcolor|zgrid]|xgrid]); 
polf(3,&savetriangle|zgrid||xgrid|{0|/0}|0)); 
polf(3,&savetriangle zgrid||xgrid ||1|/0|/0]); 


if (road|zgrid||xgrid] != 0) callobj(road{zgrid||xgrid}); 
if (target|xgrid]/zgrid| != 0) callobj(target|xgrid]|zgrid]); 
/* check if tank should be drawn now */ 


zgrid += 1; 
count += |; 


if (count >= threshold) { 
Moriae|— al. 
count = QO; 


startx -= 1; 
count = 0; 


if (startx < firstxgrid) { 


startx = firstxgrid; 
startz += threshold; 


} 
else if ((lookdir > THREE HALVES PI) && (lookdir <= SEVEN QTR PI)) 
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i tthOCT ANT aa) 
tanval = tan(lookdir+HALFPI); 
if (tanval == 0.0) 

threshold = 1000; 


else 
threshold = (int)((1.0/tanval) + 0.5); 


count = 0; 

startx = lastxgrid; 

startz = firstzgrid; 

while (startx >= firstxgrid) { 
zgrid = startz; 
xgrid = startx; 


while ((xgrid >= firstxgrid) && (zgrid >= firstzgrid)) { 


color(gridcolor zgrid]|xgrid]); 
polf(3,&savetriangle|zgrid||xgrid|/0}{0]/0)); 
polf(3,&savetriangle|zgrid||xgrid]|{1}/0]/0)); 

if (road[zgrid] xgrid] != 0) callobj(road{zgrid]|xgrid]); 

if (target /xgrid]/zgrid] != 0) callobj(target|xgrid]|zgrid}); 


xgrid -= 1; 
count += 1; 


if (count >= threshold) { 


zenid -= |; 
Counts 0: 


startz += 1; 
count = 0; 


if (startz > lastzgrid) { 


startz = lastzgrid; 
startx -= threshold; 


} 


I 
else if ((lookdir > FIVE QTR PI) && (lookdir <= THREE HALVES PI)) 


/* 6th OCTANT */ 
tanval = -tan(lookdir+HALFPI); 
if (tanval == 0.0) 

threshold = 1000; 


else 


threshold = (int)((1.0/tanval) + 0.5); 


count = 0; 
startx = firstxgrid; 
startz = firstzgrid; 
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while (startx <= lastxgrid) { 
zerid = startz; 
xgrid = startx; 


while ((xgrid <= lastxgrid) && (zgrid >= firstzgrid)) { 


color(gridcolor|zgrid||xgrid]); 
polf(3,&savetriangle|zgrid||xgrid||0}{0]/0! 
polf(3,&savetriangle|zgrid||xgrid||1](0][0] 


if (road|zgrid||xgrid| != 0) callobj(road[zgrid||xgrid}); 

if (target|xgrid||zgrid| != 0) callobj(target|xgrid||zgrid)); 
xgrid += |; 

count += 1; 


); 
Ni 


} 


if (count >= threshold) { 


zerid -= 1; 
count = 0; 
j 
} 
startz += |; 


count = QO; 


if (startz > lastzgrid) { 
startz = lastzgrid; 
startx += threshold; 


} 


} 
else if ((lookdir > PI) && (lookdir <= FIVE QTR PI)) 


/* 5th OCTANT */ 
threshold = (int)(-tan(lookdir+ HALFPI) + 0.5); 


count = 0; 

startx = firstxgrid; 

startz = firstzgrid; 

while (startz <= lastzgrid) { 
zgrid = startz; 
xgrid = startx; 


while ((xgrid >= firstxgrid) && (zgrid <= lastzgrid)) { 
color(gridcolor|zgrid||xgrid]); 
polf(3,&savetriangle|zgrid||xgrid||0){0}[0}); 
polf(3,&savetriangle(zgrid||xgrid][1|[0][0}); 


if (road|zgrid||xgrid| != 0) callobj(road|zgrid]|xgrid)); 

if (target|xgrid||zgrid] != 0) callobj(target|xgrid]|zgrid]); 
zerid += 1; 

count += 1; 
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if (count >= threshold) { 


xgrid -= 1; 
count = QO; 
Stant< ..— 1. 
count = Q; 


if (startx > lastxgrid) { 
startx = lastxgrid; 
startz += threshold; 


} 


} 
else if ((lookdir > THREE QTR PI) && (lookdir <= P])) 


/* 4th OCTANT 7 
threshold = (int)(tan(lookdir+HALFPI) + 0.5); 


count = 0; 

startx = firstxgrid; 

startz = lastzgrid; 

while (startz >= firstzgrid) { 
zerid = startz; 
xgrid = startx; 


while ((xgrid >= firstxgrid) && (zgrid >= firstzgrid)) { 


color(gridcolorzgrid||xgrid]); 
polf(3,&savetriangle|zgrid||xgrid]|{0}|0]{0}); 
polf(3,&savetriangle|zgrid|{xgrid]{1]/0]/0}); 

if (road|zgrid||xgrid| != 0) callobj(road{zgrid]|{xgrid]); 
if (target|xgrid||zgrid| != 0) callobj(target/xgrid]|zgrid)); 


7orid=— 1; 
count += |; 


if (count >= threshold) { 


XOniGie— 1; 
count = O; 
} 
startx +— I: 
count = 0; 


if (startx > lastxgrid) { 
startx = lastxgrid; 
startz -= threshold; 
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} 
else if ((lookdir > HALFPI) && (lookdir <= THREE QTR PI)) 
{ 
/* 3rd OCTANT */ 
tanval = tan(lookdir+HALFP]); 
if (tanval == 0.0) 
threshold = 1000; 


else 
threshold = (int)((1.0/tanval) + 0.5); 


count = 0: 

startx = firstxgrid; 

startz = lastzgrid; 

while (startx <= lastxgrid) { 
zerid = startz; 
xegrid = startx; 


while ((xgrid <= lastxgrid) && (zgrid <= lastzgrid)) { 


color(gridcolor|zgrid||xgrid]); 
polf(3,&savetriangle|zgrid||xgrid]|0]/0|[0)); 
polf(3,&savetriangle|zgrid||xgrid]/1]/0]/0)); 


if (road{zgrid|[xgrid] != 0) callobj(road|zgrid||xgrid]}); 

if (target|xgrid||zgrid] != 0) callobj(target|xgrid]|zgrid]); 
xerid += 1; 

count | — 1. 


if (count >= threshold) { 


zerid += 1; 
count = 0; 
startz -= 1; 
count = 0; 


if (startz < firstzgrid) { 
startz = firstzgrid; 
startx += threshold; 


} 


} 
else if ((lookdir > QTR PI) && (lookdir <= HALFPI)) 


{ 
/* 2nd OCTANT */ 


tanval = -(tan(lookdir+HALFP]I)); 
if (tanval == 0.0) 
threshold = 1000; 
else 
threshold = (int)((1.0/tanval) + 0.5); 
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count —(); 

startx = lastxgrid; 

startz = lastzgrid; 

while (startx >= firstxgrid) { 
zZerid = startz, 
xgrid = startx; 


while ((zgrid <= lastzgrid) && (xgrid >= firstxgrid)) { 


color(gridcolor{zgrid]|xgrid]); 
polf(3,&savetriangle(zgrid|{xgrid]{0}{0](0}); 
polf(3,&savetriangle|zgrid||xgrid||1]/0][0}); 


if (road|zgrid]|xgrid) != 0) callobj(road|{zgrid '||xgrid]); 

if (target xgrid|[zgrid| != 0) callobj(target|xgrid||zgrid]); 
xgrid -= |; 

count += 1; 


if (count >= threshold) { 
Zee — al 
count = 0; 


} 


Startz-— 1: 
count = 0; 


if (startz < firstzgrid) { 
startz = firstzgrid; 
startx -= threshold; 


} 
} 
else if ((lookdir >= 0.0) && (lookdir <= QTR PI)) 


/* TstrO Ga ANG >” / 
threshold = (int)(-tan(lookdir+HALFPI) + 0.5); 


count = 0; 

startx = lastxgrid; 

startz = lastzgrid; 

while (startz >= firstzgrid) { 
zgrid = startz; 
xgrid = startx; 


while ((xgrid <= lastxgrid) && (zgrid >= firstzgrid)) { 
color(gridcolor|zgrid]/xgrid!); 
polf(3,&savetriangle|zgrid]|xgrid]/0}{0]({0}); 
polf(3, &savetriangle|zgrid]|xgrid]| 1]/0]/0)); 
if (road /zgrid||xgrid] != 0) callobj(road|zgrid){xgrid]); 
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if (target|xgrid||zgrid| != 0) callobj(target|xgrid||zgrid]); 
zerid -= 1; 
count += 1; 
if (count >= threshold) { 
xgrid += 1; 
count = QO; 


} 


startx -= 1; 
ecounc — U0: 


if (startx < firstxgrid) { 


startx = firstxgrid; 
startz -= threshold; 


} 
} 


popmatrix(); 
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DIST TO LOS 


#include "gl.h" 

#include "math.h" 

float dist to los(vx,vy,vz,px,py,pz,point) 

/* compute the distance from the point "point" to the line of sight */ 


Coord vx,vy,vZ,pX,Py,p2; 
float point/3}; 


{ 
float a,b,c; /* direction numbers of line of sight */ 
float def; 
float dist; 
a = (float)(px - vx); 
b = (float)(py - vy); 
c = (float)(pz - vz); 
d = point|0! - (float)vx; 
e = point(1] - (float)vy; 
f = point/2| - (float)vz; 
dist = sqrt((up_i(e*c - f*b,2) + up i(f*a- d*c,2) + up i(d*b - e*a,2))/ 
(up i(a,2) + up i(b,2) + up i(c,2))); 
return(dist); 
} 
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DO BOUNDARY 


#include "gl.h" 

#include "math.h" 
#include "stdio.h" 
#include "fogm.h" 


#define X O 
#define Y 1 
#define Z 2 


#define DIAGONAL 0 
#define HORIZONTAL 1 
#define VERTICAL 2 


4define LOWER 0 
ddefine UPPER 1 


#define NONE 0 


#define INTERSECT 1 
#define PROPER 2 


do boundary(bound type, which triangle, xgrid, zgrid, 
bound start, bound end, left start, 

left end, right start, right end, start corner flag, 

end corner flag, polyl, vertex cnt) 


int bound type, which triangle, xgrid, zgrid; 


float bound start 3], bound end(3], left start/3!, left end{3}], 
right start(3], right end/3}; 


int *start corner flag, *end corner flag; 
float poly 1/10][3]; 


int *vertex cnt; 


int test index, cnt, index; 


float bound right[3], bound left|3], bound start edge|3}, 
bound end edge/3]; 


float vertex array|10]|3]; 
float road_poly|10)|3); 
float grid poly|10]|3}; 


int intersect cnt; 
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int intersect type, decending sort; 


float upper bound, lower bound; 


float gnd_ level(); 
int in this_poly(); 
intersect cnt = -1; 


/* compute the verticies of the road segment currently 
being worked on */ 
for (index = 0; index < 3; ++index) { 
road _poly|0||index| = left start|index|; 
road _ poly({1|{index] = left end index); 
road _ poly|2|{index] = right Pie 
road poly 3!lindex| = right start/index|; 


} 


/* compute the verticies of the grid triangle associated with 
this boundary */ 
grid aPC as (float)(xgrid*FT 100M); 
grid poly{0!/Z] = (float)((zgrid+1)*FT 100M); 
grid ~poly(1!/X] = = (float)((xgrid+1)*FT 100M); 
grid poly|1, Z| = (float)(zgrid*FT 100M); 
if (which triangle =~ UPPER) { 
grid poly 2]/X| = (float)((xgrid+1)*FT 100M); 
grid poly 2] Z) = (float)((zgrid+1)*FT 100M); 


else { 
grid poly'2! X| = (float)(xgrid*FT 100M); 
grid poly|2] Z| = (float)(zgrid*FT 100M); 
} 
if (bound type == HORIZONTAL) { 
test index = X; 
} 


else if (bound type == VERTICAL) { 
test index = Z; 
} 


else if (bound type == DIAGONAL) { 
test index = Z; 
} 


if (bound start{test index] < bound end|test index) { 
lower bound = bound start/test index]; 
upper bound = bound > _end(test index; 


} 
else { 
lower bound = bound end|test index); 
upper bound = bound start|test index]; 
} 
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/* determine points of intersection between left and right sides 
of the road and the boundary af 


line intersect2(bound start, bound end, right start, right end, 


bound right, & intersect _type); 
if (intersect type == PROPER) { 


/* intersection lies on road line segment, add intersection 
to array */ 

intersect cnt += 1; 

vertex eee erent _ent]|[X] = bound right|X]; 

vertex array|intersect cnt|/Z] = bound right[Z]; 

vertex array[intersect cnt|[/Y] = gnd_level(bound_ right|X], 

-bound right|Z]); 


} 

else if ((intersect_ type == INTERSECT) && 
(in this poly(grid poly, 3, right start)) && 
(bound _right(test _index| > > lower bound) && 
(bound right|test index] < upper bound)) { 


/* intersection point is beyond the bound of the road’s right 
line segment, but the right start point is inside the polygon so 
add the road’s right start point to the vertex array */ 


intersect cnt += 1; 

vertex “onus [inienesct _ent|[X] = right start/X; 

vertex arraylintersect cnt|/Z| = right -start|Z); 

vertex array|intersect cnt|/Y| = gnd level(right start|X], 
-right start|Z)); 


} 

else if ((intersect type == INTERSECT) && 
(in this poly(grid_ poly, 3, right end)) && 
(bound _right|test _index| > lower bound) && 
(bound right/|test index] < upper bound)) { 


/* intersection point is beyond the bound of the road’s right 
line segment, but the right end point is inside the polygon so 
add the road’s right end point to the vertex array */ 


intersect cnt += 1; 
vertex eect _cent|{X] = right end[X]; 
vertex array|intersect cntj[Z| = right end/Z]; 
vertex array|intersect cnt|/Y| = gnd level(right end|X|, 
-right_end|Z)); 
} 
line _intersect2(bound start, bound end, left start, left end, 
bound left, &intersect. ene) 
if ee cere =— FROPER) { 
/* intersection lies on road line segment, add intersection 
to array */ 
intersect cnt += 1; 
vertex arraylintersect cnt||X| = bound left/|X]; 
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vertex array|intersect cnt||Z] = bound left|Z); 
vertex arraylintersect cnt||Y| = gnd_ level(bound left|X], 
-bound left/Z)]); 


else if ((intersect type == INTERSECT) && 
(in this poly(grid poly, 3, left start)) && 
(bound left|test _index| > lower bound) && 
(bound left[test index] < upper bound)) { 


/* intersection point is beyond the bound of the road’s left 
line segment, but the left start point is inside the polygon so 
add the road’s left start point to the vertex array */ 


intersect cnt += 1; 

vertex array|intersect cnt||X] = left start|X); 

vertex array|intersect cnt|/Z| = left start(Z|; 

vertex array|intersect cnt|/Y] = gnd_ level(left start({X|], 
-left_ start|Z)); 


} 

else if ((intersect_ type == INTERSECT) && 
(in this poly(grid poly, 3, left end)) && 
(bound _left[test index} > eaves _bound) && 
(bound left|test index] < upper bound)) { 


/* intersection point is beyond the bound of the road’s left 
line segment, but the left end point is inside the polygon so 
add the road’s left end point to the vertex array * 


intersect cnt += 1; 

vertex Brett _ent]|X] = left end; X]; 

vertex arraylintersect cnt|/Z| = left _end(Z); 

vertex array|intersect cnt|/Y) = gnd iveliten _end|[X], 
-left_ end|Z}); 


} 


/* if either of the bound’s end points fall within the bounds of the 
road, add them to the array*/ 
if ((!*start corner flag) && (in this poly(road poly, 4, bound start))) { 
* put in start bound point */ 

*start corner flag = TRUE; 
intersect cnt += 1; 
vertex rates _ent}[X] = bound start[X)]; 
vertex arraylintersect cnt|/Z| = bound start|Z); 
vertex arraylintersect cnt]|/Y] = bound _start|Y]; 


if ((!*end corner flag) && (in this poly(road poly, 4, bound end))) { 
/* put in end bound point */ 
*end corner flag = TRUE; 
intersect cnt += 1; 
vertex erodes _ent][X] = bound end|X]; 
vertex array|intersect cnt|(Z] = bound end{Z]; 
vertex array|intersect cnt|{Y] = bound end|Y]; 
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j 


/* determine the point of intersection between the start and end 
bound of the road and the grid boundary */ 

line intersect2(bound start, bound end, left start, right start, 
bound start edge, & intersect _type); 
if (intersect _type == PROPER) { 

/* intersection lies on road line segment, add intersection 

to array =) 

intersect cnt += 1; 

vertex Seance cet _ent|[X] = bound start edge|X|; 

vertex arraylintersect cnt|/Z] = bound start edge[Z]; 

vertex array|intersect cnt||Y| = gnd _level(bound : start edge|X], 

-bound start edge|Z)); 
} 
line intersect2(bound start, bound end, left end, right end, 
bound end edge, &intersect type); 
if (intersect type == PROPER) { 

/* intersection lies on road line segment, add intersection 

to array tap 

intersect cnt += 1; 

vertex eerietece _ent||X) = bound end _ edge|X); 

vertex arraylintersect cnt|{Z| = bound end edge(Z]; 

vertex arraylintersect cnt|/Y] = gnd level(bound end edge(X\, 

-bound end edge[Z]); 
j 
/* put the points from the vertex array into the poly! array in 

the proper order */ 
decending sort = (bound start|test index] != lower bound); 
sort array(vertex array, intersect cnt, decending sort, test index); 


for (cnt = 0; cnt <= intersect cnt; ++cnt) { 
“vertex cnt += 1; 
poly1|*vertex cnt||X| = vertex array|cnt||X]; 
poly1|*vertex cnt|/Y| = vertex array|cnt][Y]; 
poly1|*vertex cnt]/Z| = -vertex array|cnt||Z); 
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EDIT INDBOX 


/* update the control settings of the indicator box */ 
#include "fogm.h" 
#include "gl.h" 


edit indbox(indbox, speedtag, headingtag, elevtag, altmsltag, 
zoomtag, tilttag, pantag, desigtag, speed, compassdir, 
vX, vy, vz, pan, tilt, zoom, designate) 


Object indbox; 


Tag speedtag, headingtag, elevtag, altmsltag, zoomtag, tilttag, pantag, 
desigtag; 


float speed, compassdir; 

Coord vx, vy, VZ; 

double pan, tilt; 

int designate; 

int zoom; 

{ | | 
char chspeed|5], chheading|5|, chelev{5 , chaltmsl 5]; 


float gnd_level(); 
float zoomtic, pantic, tilttic; 


sprintf{(chspeed,"%4.0f" speed); /* convert speed to string */ 
sprintf(chheading,'%3.0,compassdir); /* convert heading to str */ 
sprintf(chelev,"%%4.0f",vy - gnd_level(vx,vz)); /* convert elev AGL to str */ 
sprintf(chaltms],"%4.0f" vy); /* convert alt MSL to str */ 

/* compute new location for zoom, pan, and tilt indicators */ 


zoomtic = zoom * -0.2766 + 222.128; 
tilttic = tilt * 721.92682 + 365.0; 
pantic = pan * -721.92682 + 435.0; 


editobj(ind box); /* update the indicator display */ 
objreplace(speedtag); 
charstr(chspeed); 
objreplace(headingtag); 
charstr(chheading); 
objreplace(elevtag); 
charstr(chelev); 
objreplace(altmsltag); 
charstr(chaltms]); 
objreplace(zoomtag); 
move2(28.0,zoomtic); 
objreplace(tilttag); 
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move2(42.0,tilttic); 

objreplace(pantag); 

move2(pantic,27.0); 

objreplace(desigtag); 

cmov2i(designate ? 10 : 19,10); 

charstr(designate ? "DESIGNATE" : "REJECT"); 


closeobj(); 
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EDIT NAVBOX 


#include "fogm.h" 
#include "math.h" 
#include "gl.h" 


edit navbox(navbox, arrowtag, vx, vz, direction, firstxgrid, firstzgrid, 
lastxgrid, lastzgrid) 

Object navbox; 

Tag arrowtag; 

Coord vx, vz; 

double direction; 

short firstxgrid, firstzgrid, lastxgrid, lastzgrid; 


{ 


Coord arrowx, arrowy, larrowx, larrowy, rarrowx, rarrowy; 


/* compute coordinates of arrow line segments for nav control box */ 
arrowx = vx + cos(direction) * 2.0 * FEETPERGRID; 

arrowy = vz - sin(direction) * 2.0 * FEETPERGRID; 

larrowx = arrowx + cos(direction - 2.3561945) * FEETPERGRID; 
larrowy = arrowy - sin(direction - 2.3561945) * FEETPERGRID; 
rarrowx = arrowx + cos(direction + 2.3561954) * FEETPERGRID; 
rarrowy = arrowy - sin(direction + 2.3561945) * FEETPERGRID; 


/* update the contour map display with new info ag 
editobj(navbox); 

objreplace(arrowtag); 

move2(vx,vz); 

draw2(arrowx, arrowy); 

draw2(larrowx, larrowy); 

move2(arrowx, arrowy); 

draw2(rarrowx, rarrowy); 

rect(firstxgrid*FT 100M,-firstzgrid*FT 100M, 
(lastxgrid+1)*FT 100M, (-lastzgrid-1)*FT 100M); 
closeobj(); 
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EXPLOSION 


#include "gl.h" 


explosion() 


{. 


int 1,}; 


pushviewport(); 

view port (0,1023,0,767); 

color(BLACK); 

clear(); 

swapbuffers(); 

color(RED); 

clear(); 

swapbuffers(); 

swapbuffers(); 

color( YELLOW); 

clear(); 

swapbuffers(); 

swapbuffers(); 

color{(RED); 

clear(); 

swapbuffers(); 

swapbuffers(); 

color(YELLOW); 

clear(); 

swapbuffers(); 

swapbuffers(); 

color(RED); 

clear(); 

swapbuffers(); 

swapbuffers(); 

for (1 = 0; 1 < 100000; 1++) 
for p= 07) < 10: j++): 

popviewport{); 
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FOGM (MAIN) 


/* fogm.c -- an IRIS-2400 program by Doug Smith & Dale Streyle 
It reads in a 10km x 10km section of a terrain map, computes a lighting 


and shading model for the terrain, and allows overflight a 
#include "gl.h" /* get the graphics defs ey 
#include "device.h" /* get the graphics device defs */ 
#include "fogm.h" /* constants : 
#include "math.h" /* math function declarations a 
#include "get.h" /* monitor type include file =) 


#include "stdio.h" 

#include "sys/signal.h" /* used for screen dump utility */ 
#include <sys/types.h>  /* contains the time sturcture tms */ 
#include <sys/times.h> /* for time calls a 


short gridpixel/100]|100|]; /* DMA elevation and vegatation data oa 
float savetriangle 99|/99]|2 [3 [3]; 

long gridcolor|99!|99}; 

Object road|99]/99]; 

Object target|99}[99}; 


float ground plane/4|/3); 

long gnd plane color; 

float tgt_pos|MAX TGTS]|(3}; 

short tgt grid idx{MAX TGTS]/2]; 

short tgt dir MAX TGTS|, tgt total = 0; 


float randx, randy, randz; /* random offsets from tank reference point */ 
int framecnt; 

float min elev, max elev; 

Coord tankx, tanky, tankz; 

float frames _sec|1000]|2}; 


main() 


{ 
int greyscale = FALSE; /* FALSE = color, TRUE = greys */ 


int designate; /* boolean indicating desig/reject status ) 
int flying = TRUE; /* boolean controlling flying loop Te 
int active = TRUE: /* boolean controlling main program loop as 


int nbyte, socket, connect client(); /* networking variables & subroutine */ 
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struct tms timestruct; /* structure for real-time clock calls ey 


int tgt idx; /* index of designated target oy 
double direction; /* direction of travel in radians =o 
float speed; /* speed of travel in knots a 

float compassdir; /* desired direction of travel in compass deg */ 
int fovy = 550; /* field of view in perspective command */ 


double pan = 0.0, 
tilt = -15.0 * DTOR; /* pan and tilt angles of 


/* contour map, indicator, instruction */ 
Object contour, navbox, indbox, instrbox; 


Object tank, pre | obj 7]; 


Tag headingtag, elevtag, speedtag, zoomtag, arrowtag, tilttag, pantag; 
Tag desigtag, altmsltag, pre | tag 6]; 


Colorindex unmask; 


Coord vx, vy, vz; /* viewer x y and z coordinates */ 

Coord px, py, pz; /* reference x y z coordinates for lookat */ 
Coord tgtx, tgty, tgtz; /* targeted position on tank * / 

float randseed(); /* random number generator initialization */ 
int frames = -1; 


long seconds, lastseconds, totalseconds = 0; 
int numpolys; 

float elapsed; 

int idx; 


FILE *fopen(), *fp; 


/* first and last x and z indexes of the grid objects to draw */ 
short firstxgrid, firstzgrid, lastxgrid, lastzgrid; 


readdata(); /* read the data file into the gridpixel array */ 


/* get socket number for networking */ 


/*if (NETWORKING) socket = connect client("npscs-iris1",3); * / 


init _iris(); /* initialize the iris * / 
unmask = (1<<getplanes()) - 1; 
writemask(unmask); 


randseed(times(&timestruct)); /* seed the random # generator * / 
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init tgts(); /* define targets */ 

Screen Dump(SCREENDUMP); /* enable screen dumping */ 
billboard(); /* produce intro screen */ 
colorramp(greyscale, TRUE); /* build all color ramps */ 


makescreens(pre | obj, pre | tag); /* build objects for prelaunch */ 
makemap(&contour); /* build map object */ 
pre 1 obj}; CONTOUR] = contour; 


prelaunch(&vx, &vy, &vz, &direction, &compassdir, 
&active, pre | obj, pre | tag); 


if (active) { 
maketank(&tank); /* build object for a tank */ 


build road(); /* build the objects that comprise the roads */ 


/* process terrain data to build polygons and compute lighting */ 
buildterrain(); 


/* build object for the navigation display contour map */ 
drawnavbox(&navbox, &arrowtag); 


/* build an object for the indicator box */ 
makeindbox(&indbox, &headingtag &elevtag,&altmsltag,&speedtag, 
&zoomtag,&tilttag, &pantag, &desigtag); 


makeinstrbox(&instrbox); /* build object for control instruction box */ 
\ /* end of if (active) block */ 
while (active) { 

framecnt = 0; 


/* initialize the operator controls (mouse and dials) */ 
init controls(&pan, &tilt, &fovy, vy, greyscale, compassdir); 


pushviewport(); 
viewport(0,1023,0,767); 
color(SKYBLUE); 

clear(); 

popviewport(); 
callobj(instrbox); 
callobj(indbox); 
editobj(contour); 
objreplace(STARTTAG); 
viewport (768, 1023 ,512,767); 
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closeob)(); 
callobj(contour); 
swapbuffers(); 
callobj(instrbox); 
callobj(contour); 
editobj(contour); 
objreplace(STARTT AG); 
viewport (0,768 ,0,768); 
closeobj(); 


flying = TRUE; /* missile is flying */ 
designate = TRUE; /* a target can be designated */ 


while(flying) { /* until tgt is hit or 3-button exit */ 


/* get values from user contols (mouse and dials) */ 
read controls(&designate, &greyscale, &flying, &active, 
&speed, &direction, &compassdir, &vy, 

&pan, &tilt, &fovy); 


/* calculate which target was closest to the line of 
sight */ 

if (!'designate) { 
nearest tgt(vx,vy,vz,px,py,pz,&tgt idx); 

} 


/* update targets’ positions */ 

get tgt posit(socket, designate, tgt idx, &tgtx, &tgty, &tgtz, tank); 
/* update missile position */ 

update missile posit(&direction, &compassdir, speed, 
designate, tgtx, tgty, tgtz, 

&vx, &vy, &vz, &flying); 


/* update camera lookat position */ 
update look posit(direction, pan, tilt, vx, vy, vz, 
tgtx, tgty, tgtz, designate, &px, &py, &pz); 


/* determine which polygons need to be drawn */ 
view bounds(vx, vy, vz, px, py, pz, tilt, fovy, 
&firstxgrid, &firstzgrid, &lastxgrid, &lastzgrid); 


/* edit control display objects to reflect new values * / 

edit navbox(navbox, arrowtag, vx, vz, direction, firstxgrid, 
firstzgrid, lastxgrid, lastzgrid); 

edit _indbox(indbox, speedtag, headingtag, elevtag, altmsltag, 
zoomtag, tilttag, pantag, desigtag, speed, 

compassdir, vx, vy, vz, pan, tilt, fovy, designate); 


/* display the 3-D view of the terrain as seen by 
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the camera a, 
display terrain(vx, vy, vz, px, py, pz, fovy, 
firstxgrid, firstzgrid, lastxgrid, lastzgrid); 


/* display the control boxes */ 
writemask(SA VEMAP); 
callobj(navbox); 
writemask(unmask); 


callobj(indbox); 
swapbuffers(); 


seconds = times(&timestruct); 
numpolys = (lastxgrid - firstxgrid) *(lastzgrid-firstzgrid) *2; 
elapsed = (float)(seconds - lastseconds) /60.0; 
if ((frames >= 0) && (frames < 1000) ){ 
frames sec|frames||0| = (float)numpolys; 
frames sec|frames|(1] = 1.0/elapsed; 
} 
totalseconds += (seconds-lastseconds); 
if (totalseconds > 7200) { 
compactify(); /* do garbage collection every 2 mins */ 
totalseconds = 0.0; 
} 
lastseconds = seconds; 
frames += 1; 


\  /* end of flying loop */ 


if (active) { /* explode & restart */ 
explosion(); 
prelaunch(&vx, &vy, &vz, &direction, &compassdir, 
& active, pre 1 obj, pre | tag); 


} 


\ /* end of active loop */ 


/* write out performance stats */ 
fp = fopen("speed.data", "w"); 
if (frames > 999) frames = 999; 
for (idx = 0; idx <= frames; ++idx) { 
fprintf(fp,"%.2f %.2f0, frames sec{idx||0|, frames sec{idx||1]); 
} 
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/* gracefully exit */ 
if (NETWORKING) close(socket); 
setmonitor(HZ60); 
color(BLACK); 
clear(); 
swapbuffers(); 
clear(); 
gexit(); 
textinit(); 
exit(); 

} /* end of main */ 
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FILES.H 


/* These are the files which contain data for the terrain elevations 
and roads */ 

#define TERRAIN FILE "/work/terrain/tenkmsq.dat" 

#define ROAD FILE "/work/terrain/Road.data" 


FOGM.H 

#define elev_mask Ox 1 fff /* mask to obtain elev value from datum  */ 
#define veg mask  0x0007 /* mask to obtain vegatation value from 

shifted datum */ 
#define RD 0 /* code for reading a file in "open" i) 
#define MAX 2800 #=/* max elev (ft) in contour map of 
#define MIN 967 /* min elev (ft) in contour map a 
#define SKYBLUE 4095 /* color index for sky color ea 
#define ROADGREY 850 /* color index for the road | 


#define DELTAFOVY 50 /* field of view (zoom) increment of 5 deg */ 


# define PI 3.1415927 
#define TWOPI 6.2831853 
#define HALFPI 1.5707963 


#define THREE HALVES PI 4.7123889 

define QTR PI 0.7853982 

#define THREE QTR PI 2.3561945 

#define FIVE QTR PI 3.9269908 

#define SEVEN QTR PI 5.4977871 

#define RTOD 57.29578 /* radians to degrees conversion factor */ 
#define DTOR 0.0174533 /* degrees to radians conversion factor */ 
#define FPS TO KTS 35.525148 /* convert feet per 60th seconds to knots */ 
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#define PANSENS 30.0 —_//* scale factors (sensitivity) for 


navigaion controls (mouse and dials) */ 
#define SPEEDSENS —_20 
#define TILTSENS 50.0 
#define DIRSENS 20.0 


#define MAXLOOKDIST 32808.0 /* maximum distance that the camera can 
look ahead in feet oy 


#define FEETPERGRID 3280.8 /* number of feet in 1000 meters a 
#define ALTSCALE 1.05 __//* altitude expansion factor, altitudes are 
raised to this power to give an 


exagerrated effect “i 


#define NUMXGRIDS 10 /* number of 1K grid squares in the East- 
West direction af 


#define NUMZGRIDS 10 /* number of 1k grid squares in the North- 


South direction af 
#define FT 10K 32808 /* number FT in 10Km a | 
#define FT 100M 328.08 /* number FT in 100m a) 
#define GRID FACTOR 13.03781 /* conversion factor ) 
#define TV 0 /* 0 for SGI monitor, 1 for TV =) 
#define SCREENDUMP i /* 1 to enable screen dumping, 0 otherwise */ 
#define NETWORKING 0 /* 1 for target networking, 0 otherwise */ 
#define INIT PAN 0 /* initial, min and max pan angles in deg. */ 
#define MIN PAN -25 
#define MAX PAN 25 


#define INIT TILT -15 /* initial, min and max tilt angles in deg.*/ 
#define MIN TILT a2 
#define MAX TILT 15 


#define MAX ALT 17000 /* maximum altitude for missle  */ 
#define MIN ALT 0 /* minimum altitude for missle */ 
#define INIT SPEED — 200 /* init, min and max spd (kts) for missle */ 


#define MIN SPEED 0 
#define MAX SPEED 400 


#define INIT FOVY 550 /* initial field of view in tenth degrees */ 
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# define 
# define 
#define 
#define 
#define 
# define 
#define 


CONTOUR 
SCREEN1 
SCREEN2 
SCREENS 
INSTR 
STATS 
FLTPATH 


#define LAUNCH 


#define 
#define 
# define 
#define 
#define 


#define MAX TGT COLOR 


TARGET 
DIR 

HEAD 
TGT 4 
MISSILE 


on te Ww he — 


0 


#define MIN TGT COLOR 668 


#define 


#define 


MAX TGTS 


SA VEMAP 


100 


/* Indicies for array obj 


/* Indicies for array tag 


847 


0x00CO0 
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7 


/ 


GAMMARAMP 


/* This routine puts a gamma-corrected color ramp into the color map. = */ 
#include <math.h> 
gammaramp(gammaconst,firstcolor,ncolors, 
brightred,brightgreen,brightblue, 
darkred ,darkgreen,darkblue) 
float gammaconst; /* Strength of Gamma correction (try 1.0) */ 
long firstcolor; /* index number of the first color to set */ 
long ncolors; /* the number of colors to set */ 
long brightred,brightgreen,brightblue; /* the bright end of the ramp */ 
long darkred,darkgreen,darkblue; /* the dark end of the ramp */ 
{ 

long i; /* temp loop index */ 

float scl; /* scale factor for gamma correction */ 


long gcred,gcgreen,gcblue; /* gamma corrected colors */ 


for(i=O0; i < ncolors; i++) /* for all colors...* / 


{ 


/* compute the scale factor */ 
scl = pow((float)i/(float)(ncolors-1) , 1.0/gammaconst); 


/* compute the gamma corrected colors */ 

gcred = scl * (brightred - darkred) + darkred; 

gcgreen = scl * (brightgreen - darkgreen) + darkgreen; 
gcblue = scl * (brightblue - darkblue) + darkblue; 


mapcolor(firstcolor+i, gcred, gcgreen, gcblue); /* set the color */ 
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GET TGT POS 


/* get targets’ positions from irisl if networking. Otherwise moves 10 targets 
in straight lines, reversing when they hit an edge a 


#include "fogm.h" 

#include "gl.h" 

#include "math.h" 

#include <sys/types.h>  /* contains the time sturcture tms */ 
#include <sys/times.h> /* for time calls */ 


get tgt posit(socket,designate,tgt idx,tgtx,tgty,tgtz,tank) 


int socket, designate, tgt idx; 
float *tetx. “tety. tetz- 
Object tank; 


{ 
extern float tgt pos|MAX TGTS]|3]; 


extern float randx, randy, randz; 

extern Object target |99]/99]; 

extern short tgt grid idx|MAX TGTS]/2); 
extern short tgt total, tgt dir. MAX TGTS]; 
short 1, tgt_ num; 

int nbyte, addi(); 

float gnd_ level(), dir, dx, dz, distance; 

long dist, d2; 

static long seconds; 

static long lastsec = -999;  /* -999 is flag to indicate no value */ 
struct tms timestruct; 


seconds = times(&timestruct); 


if (lastsec == -999) /* compute distance targets move ahead */ 
distance = 0.0; 
else 
distance = (float)((15.0/FPS TO KTS) *(seconds - lastsec)); 
lastsec = seconds; /* save for next pass */ 


for (i = 0;1 < tgt total;i++)  /* delete targets from old positions */ 
if (target|tgt grid idx{i){O|||tgt grid idx|il/1]]) { 
delobj(targetitgt grid idx/i][O]|[tgt grid idx{il/1]]); 
targetitgt grid idx(il{O)|[tgt grid idx/il[1|] = 0; 


} 


if (NETWORKING) { 
nbyte = read(socket, &tgt total, sizeof(tgt total)); 
for (i = 0;1 < tgt total; i++) { 
nbyte = read(socket, &tgt grid idx|i|[0], sizeof(short)); 
nbyte = read(socket, &tgt grid idx|i]|1], sizeof(short)); 
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nbyte = read(socket, &tgt pos/i]|0), sizeof(float)); 
nbyte = read(socket, &tgt pos|il[1], sizeof(float)); 
nbyte = read(socket, &tgt pos/i][2), sizeof(float)); 
nbyte = read(socket, &tgt dir[i], sizeof(short)); 


else { 
tgt total = 10; 
fore 0 tet cotalics, 
dir = (float)(tgt dir|i] / 10) * DTOR; 
tgt pos|i][0] += cos(dir) * distance; 
tgt pos|i][2| -= sin(dir) * distance; 
tgt grid idx(i||0] = (short)(tgt pos|i][0|/FT 100M); 
tgt grid idx{i][1] = (short)(-tgt pos|i][2|/FT 100M); 
if ((tgt pos|i]/0] > FT 10K) || (tgt pos|il/0] < 0)) { 
if (tgt dir[i] > 1800) tgt dir[i] -= 1800; 
else tgt dir i] += 1800; 
tgt_pos|i]/1] = 0.0; 


else if ((tgt pos|i]/2) <-FT 10K) || (tgt pos|i][2] > 0)) { 
if (tgt dir|i] > 1800) tgt dir|i] -= 1800; 
else tgt dir'i] += 1800; 
tgt pos|i]|1] = 0.0; 


else tgt_pos|i|{1] = gnd_ level(tgt pos|i]|0|, tgt pos|i}[2)); 
} 
} 
if (!designate) { 
if (NETWORKING) { /* find which target is designated * / 
dist = up i((float)(tgt pos[0]/0] - *tgtx),2) + 
up i((float)(tgt pos|O|{2) - *tgtz),2); 
tgt idx = 0; 
for (1 = 1;1 < tgt total; i++) { 
d2 = up i((float)(tgt pos il[O] - *tgtx),2) + 
up i((float)(tgt pos/il/2] - *tgtz),2); 
if (d2 < dist) { 
dist = diz: 
tgt idx = (int)i; 


} 
} 
*tgtx = tgt pos|tgt idx][0] + randx; 
*tgty = tgt pos|tgt idx||1] + randy; 
“tgtz = tgt pos|tgt idx||2] + randz; 
} 
tgt num = tgt total; 
for (i = 0;1 < tgt num; i++) { 
dx = tgt pos|i][O| - (float)tgt grid idx|i)[0] * FT 100M; 
dz = (float)(-tgt grid idx|i]/1]) * FT 100M - tgt pos|i}[2); 
if (dx < 15.0) - 7 
if (dz < 15.0) { 
add 1(i,-1,0); 
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add1(i,-1,-1); 
add1(i,0,-1); 


else if (dz > 313.0) { 
add1(1,0,1); 
add1(i,-1,1); 
add1(1,-1,0); 
j 
else add1(1,-1,0); 
else if (dx > 313.0) 
if (dz < 15.0) { 
add1(i,0,-1); 
add1(i,1,-1); 
add1(i,1,0); 
i 
else if (dz > 313.0) { 
add1(1,1,0); 
add1(i,1,1); 
add1(i,0,1); 
} 
else add1(i,1,0); 
else if (dz < 15.0) add1(i,0,-1); 
else if (dz > 313.0) add1(i,0,1); 
} 
for (1 = 0;1 < tgt total;i++) /* add targets to new positions */ 
if (target[tgt grid idx|[i][O]|[tgt grid idx|i}[1]]) { 
editobj(target|tgt grid idx|i]/0]//tgt grid idx[i][1]]); 
pushmatrix(); 
translate(tgt pos|i) O],tgt pos|i][1],tgt pos|ij|2]); 
rotate(tgt dirii], Y’); 
callobj(tank); 
popmatrix(}; 
closeobj(); 


else { 
target|tgt grid idx|i{O]|/tgt grid idx{i/1]) = genobj(); 
makeobj(target(tgt grid idx|i|/0|| tgt grid idx{i}{1}]); 
pushmatrix(); 
translate(tgt pos/(i]{O],tgt pos{i|{1|,tgt pos|i}[2]); 
rotate(tgt dir|il, LY a\aon ae 
callobj(tank); 
popmatrix(); 
closeob)(); 


addi(tgt num,x,z) 
short tgt num,x,2Z; 
{ 
extern float tgt_ pos|/MAX TGTS|(3); 
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extern short tgt grid idx|MAX TGTS]|2|; 
extern short tgt total, tgt dir|MAX TGTS); 
short 1; 


tgt pos(tgt total][0] = tgt pos(tgt num||0]; 
tgt pos|tgt total]{1] = tgt pos|tgt num|/1]; 
tgt pos|tgt total][2] = tgt pos|tgt num|/2]; 
tgt dir[tgt total] = tgt dir{tgt num]; /* copy dir for "new" tgt */ 
tgt grid idx|tgt total]{O] = tgt grid idx|tgt num] 0] + x; /* set pos in */ 
tgt grid _idx|tet _total][1] = tgt grid idx{tgt num]/1] + 2; /* new grid sq */ 
for (i =0;i< 2;i++){  /* reset if new grid sq outside 10km square */ 

if emetic deeemcrall i < 0) tgt grid idx|tgt total]|/i] = 0; 

if (tgt grid idx|tgt total]|i] > 98) tgt grid idx|tgt total][i] = 98; 


/* copy pos. for "new" tgt */ 


} 


tgt total ++; 
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GND LEVEL 


#include "math.h" 
#include "fogm.h" 
#define X O 
#define Y 1 
#define Z 2 


float gnd_ level(vx, vz) 


float vx, vz; 


{ 


extern short gridpixel|100)|100}; 

float interp elev(); 

float grid level(); 

float point 3], nw corner|3], ne corner,3], sw corner{3], se corner|3]; 
float intersect|3/; - 

float elev; 

int xgrid, zgrid, intersect type; 


/* determine which triangle the point falls in */ 
xgrid = (int)(vx/FT 100M); 

zgrid = (int)(-vz/FT 100M); 

if (xgrid < 0) xendu—.0, 

if (xgrid > 98) xgrid = 98; 

if (zgrid < 0) 2gridi— 0. 

if (zgrid > 98) zgrid = 98; 

point|X] = vx; 

point|Z] = -vz; 

nw corner X] = (float)(xgrid*FT 100M); 

nw corner|Z]| = (float)((zgrid + 1)*FT 100M); 
elev = gridpixel|zgrid+1||xgrid] & elev mask; 
nw corner Y| = pow(elev, ALTSCALE); 

sw corner|X| = (float)(xgrid*FT 100M), 

sw corner|Z] = (float)(zgrid*FT 100M); 

elev = gridpixel(zgrid||xgrid] & elev_mask; 

sw _corner|Y| = pow(elev, ALTSCALE); 

ne corner|X| = (float)((xgrid+1)*FT 100M); 
ne corner|Z| = (float)((zgrid+1)*FT 100M); 
elev = gridpixel{zgrid+1]|[xgrid+1] & elev mask; 
ne corner|Y| = pow(elev, ALTSCALE); 

se corner|X]| = (float)((xgrid+1)*FT 100M); 
se corner|Z| = (float)(zgrid*FT 100M); 

elev = gridpixel|zgrid||xgrid + 1| & elev_mask; 
se corner/Y| = pow(elev, ALTSCALE); 


if (-vz < (nw corner|Z] - (vx - nw corner|X]}))) { 
/* point is in the lower triangle */ 


/* find the point of intersection of a line through vx,vz 
and the sw corner with the diagonal */ 
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line intersect2(sw corner, point, nw corner, se corner, intersect, 
&intersect type); 


/* find the elevation of the intersection on the diagonal */ 
intersect|Y] = interp elev(nw corner, se corner, intersect); 


/* find the elevation of the point vx, vy */ 
return(interp elev(sw corner, intersect, point}); 


else { 
/* point is in the upper triangle */ 


/* find the point of intersection of the diagonal with a line 
through th ne corner and the point */ 

line intersect2(ne corner, point, nw corner, se corner, intersect, 

&intersect_ type); 


/* find the elevation of the intersection on the diagonal */ 
intersect| Y] = interp elev(nw corner, se corner, intersect); 


/* find the elevation of the point vx, vz */ 
return(interp elev(ne corner, intersect, point)); 
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IN THIS POLY 


#include "gl.h" 


#define X 0 
#define Y 1 
#define Z 2 


#define PROPER 2 


int in this poly(polygon, num vertex, point) 
float polygon|[10]|3]; 

int num _ vertex; 

float point 3); 


{ 
int index; 
int pt in, intersect type; 
int num crossings; 
float max xX, Max z, min xX, min 2; 
float intersect/ 3]; 
float old intersect|3); 
float start test line|3]; 


max x = polygon 0]/X); 
min x = polygon(0||X); 
max z= polygon|0 |Z]; 
min z= polygon|0||Z/; 


for (index = 1; index < num vertex; ++index) { 
if (polygon[index]|X] < min x) min x = polygon|index||X]; 
if (polygon!index]/X] > max x) max x = polygon|index||X]; 
if (polygon{index||Z| < min _z) min z = polygon|index)|Z]; 
if (polygon|index||Z| > max z) max z = polygon[index|(Z]; 


} 


if ((point|X] < max_x) && (point/X) > min x) && (point|Z] < max z) && 
(point|Z] > min z)) { 


/* point may be polygon, test further by constructing a vertical line 
from the point to a point outside the polygons bounds. Count the number 
of times this line crosses a side of the polygon. If it crosses an 
odd number of times the point is in the polygon, otherwise it is 
outside the polygon */ 


start test line[X] = point(X]; 
start test line(Z| = max z + 1000.0; 


num crossings = 0; 
old intersect|X] = -999.0; 
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old intersect|Z] = -999.0; 
for (index = 0; index < num vertex -1; ++index) { 
line intersect2(start test line, point, &polygon|index||0}, 
&polygon{index+1]{0], intersect, &intersect type); 
/* if a proper intersection exists and it is not the same point 
as the previous intersection (i.e it didn’t intersect a vertex), 
then add one to the number of crossings */ 
if ((intersect type == PROPER) && ((intersect|X| != old_intersect|X}) 
|| (intersect|Z] != old intersect|Z]))) num crossings += 1; 
old intersect|X] = intersect|X]; 
old intersect(Z] = intersect|Z]; 
} 
line intersect2(start test line, point, &polygon|num_vertex-1](0], 
&polygon|0|/0], intersect, &intersect type); 
if (intersect type == PROPER) num crossings += 1; 


/* if the number of crossings is even, the point was outside 
pt in = ((num crossings % 2) != 0); 


return(pt in); 


else { 


return(FALSE); 
} 
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INIT CTRLS 


/* initialize the operator controls */ 


#include "fogm.h" /* fogm constants * / 

#include "device.h" /* graphics device definitions */ 
#include "gl.h" /* graphics routine definitions * / 
#include "math.h" /* math function definitions */ 


init _controls(pan, tilt, fovy, alt, greyscale, compassdir) 


double *pan; /* initial pan angle in radians */ 

double *tilt; /* initial tilt angle in radians */ 

int *fovy; /* initial field of view in tenths of degrees * / 
Coord alt; /* initial altitude of missile */ 

int greyscale; /* initial value of greyscale boolean */ 

float compassdir; /* initial compass direction */ 

{ 


*pan = INIT PAN * DTOR; 
“tilt = INIT STILT  BLor;: 
*fovy = INIT FOVY; 


/* set initial, min, and max values for mouse & dials */ 


setvaluator(MOUSEX,(short)(INIT PAN*PANSENS),(short)(MIN PAN*PANSENS), 
(short)(MAX PAN*PANSENS)); 


setvaluator(MOUSEY ,(short)(INIT_TILT*TILTSENS),(short)(MIN TILT*TILTSENS), 
(short)(MAX_ TILT*TILTSENS)): 


setvaluator(DIALO,(short)(compassdir* DIRSENS), (short}(-360*DIRSENS), 
(short)(720* DIRSENS)); 


setvaluator(DIAL4,(short)alt,MIN ALT,MAX ALT); 
setvaluator(DIAL2, (short)(INIT SPEED*SPEEDSENS), 
(short)(MIN SPEED*SPEEDSENS), 

(short}(MAX SPEED*SPEEDSENS)); 


setvaluator(DIAL3,greyscale,0,1); 
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INIT IRIS 


/* Initialize the graphics environment for the iris workstation */ 


#include "gl.h" /* graphics definitions */ 
#include "get.h" /* monitor type definitions */ 
#include "fogm.h" /* fogm constants */ 


init _iris() 


{ 


long chunk; /* number of bytes be which objects 
increment */ 
ginit(); /* initialize the IRIS system “) 
doublebu ffer(); /* put the IRIS into double buffer mode */ 
enunk = 126: 
chunksize(chunk); 
gconfig(); /* (means use the above command settings) */ 
if (TV) { 
setmonitor(NTSC); /* choose tv or SGI monitor my 
fontdef(1,"T'V font"); 
font(1); 


else setmonitor(HZ60); 


cursoff(); /* turn off the cursor zy 
backface(TRUE); /* turn on backface polygon removal / 
color(BLACK); 

clear(); 

swapbuffers(); 
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INIT TGTS 


#include "fogm.h" 
#include "gl.h" 


init tgts() 


} 


extern short tgt total; 
extern Object target|99]|/99]; 
short x, y; 

int init tgt(); 


for (x = 0; x < 99; x++) for (y = 0; y < 99; y++) target|x]||y] = 0; 

if ((NETWORKING) { 
tgt total = 10; 
init tgt(0,9.8,3.5,1295); 
init tgt(1,9.5,3.5,1295): 
init tgt(2,9.4,3.1,1295); 
init tgt(3,9.8,0.5, 1800); 
init tgt(4,9.5,0.0, 1355); 
init tgt(5,8.0,0.0, 1445); 
init tgt(6,4.0,0.0,1450); 
init tgt(7,0.0,0.5,450); 
init tgt(8,9.5,9.8, 2700); 
init tgt(9,9.8,8.5, 1800); 


init tgt(tgt num,xoffset ,zoffset direction) 


short tgt num, direction; 
float xoffset, zoffset; 


{ 


extern short tgt_dir MAX TGTS); 
extern float tgt pos;|MAX TGTS]|3}; 


tgt pos|tgt num|/0] = xoffset * FEETPERGRID; 
tgt pos|tgt num||2] = -zoffset * FEETPERGRID; 


tgt dir{tgt num] = direction; 
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INTERP ELEV 


#include "math.h" 


#define X O 
#define Y 1 
#define Z 2 


float interp elev(line start, line end, point) 

float line start(3], line end[3], point/3]; 
long float line deltax, line deltaz, point deltax, point deltaz; 
float line length, dist to _point; 


float interpolation; 


line deltax = (long float)(line end X| - line start/X}); 
line deltaz = (long float)(line end(Z] - line start|Z]); 


point deltax = (long float)(line start|X] - point|X)); 
point deltaz = (long float)(line start|{Z] - point|Z]); 


line length = (float)hypot(line deltax, line deltaz); 
dist to point = (float)hypot(point « deltax, point deltaz); 


interpolation = line start/Y|] + ((line end[Y] - line start/Y]) * 
(dist to point/line length)); 


return(interpolation); 
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/* this is file lightorient.c */ 


i= 


LIGHTORIENT 


lt is a routine that computes lighting for a polygon based 
upon the angle between the Normal vector of the polygon 
and the direction to the light source. 


lightorient(xyz,ncoords,ax,ay,az,1x,ly,lz,colormin,colormax,colortouse) 


xyz||/3] = floating coords of the polygon. 


ncoords = number of coordinates. 


ax,ay,aZ = interior point of the whole object. Used to determine 


outward facing normal of the polygon. This is the same 


point of reference that would be used for backface 


polygon removal. 


Ix ,ly,lz = vector pointing in direction of the light source. 


colormin, colormax = indices used for the colors assigned to this 


polygon. The user is responsible for setting 


up the color ramp. 


colortouse = returned color used to light the polygon. 


Note: the routine also puts the polygons out ordered counterclockwise 


with respect to the interior point for ease of backface polygon 


removal. 
| 
#include <math.h> 
#include <gl.h> 
#define MAXCOORDS 80 


#define PIDIV2 1.570796327 


float txyz/MAXCOORDS]/(3); 


/* temp coord hold */ 


lightorient(xyz,ncoords,ax,ay,az,lx,ly,1z,colormin,colormax,colortouse) 


float xyz(|(3]; 
long ncoords; 
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float ax,ay,az; /* interior point of the whole object. */ 


float Ix,ly,lz; /* direction to the light source oH 


long colormin,colormax; /* color min/max indices 7) 


long *colortouse; /* color used to light the polygon (return value) */ 


{ 


long i,j; /* loop temps */ 
long npoly orient(); /* direction test function */ 


float v1/3],v2/3]; /* vectors used to compute 
the polygon’s normal */ 


float normal|/3); /* the polygon’s normal */ 

float normalmag; /* normal’s magnitude */ 

float lightmag; /* magnitude of the light vector */ 
double dotprod; /* dot product of N and L */ 
float radians; /* angle between N and L */ 


/* check the number of coords in the input array */ 


if(ncoords > MAXCOORDS) 
{ 


printf("LIGHTORIENT: too many coords passed to me! = %d0,ncoords); 
exit(1); 


} 


/* orient the polygon so that its counterclockwise with respect 
to the interior point */ 


if(npoly orient(ncoords,xyz,ax,ay,az) == 1) 
/* the polygon is clockwise, reverse it. * / 
for(i=0; i < ncoords; i=i+1) 

{ 
for(j=0; j < 3; j=j+1) 
txyz|il[j] = xyz|ncoords-i-1]|j]; 


} 


for(i=0; 1 < ncoords; ++) 
for (j=0; j < 3; ++)) 
xyzi]|j] = txyzil]j]; 
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} 


/* the coordinates are ordered counterclockwise in array xyz */ 


/* compute the normal vector for the polygon using the first 
three vertices...* / 


/* compute the first vector to use in the computation */ 
v1/0} = xyz/[2][0] - xyz[1][0]; 
vil] = xyz(2][1] - xyz[1][1]; 
vi/2) = xy2[2][2] - xyz[1][2]; 


| 


/* compute the second vector to use in computing the normal */ 
v2.0: = xyz[0][0] - xyz/[1][0]; 
va(t} = xyz(ola] - xy2[a]{1} 
v2/2) = xyz[0}[2] - xyz/1][2]; 


/* the normal is vl x v2 */ 

normal/0] = v1{1]*v2[2] - v1[2]*v2/1]; 
normal(1] = v1[2]*v2[0] - v1[0]*v2/[2]; 
normal|2] = vt(0]*v2[1] - v1{1]*v2(0); 


/* compute the magnitude of the normal */ 
normalmag = sqrt((normal/0 *normal{0!)+(normal|{1]*normal|{1|)+ 
(normal|2|*normal|2])); 


/* check the magnitude of the normal */ 
if(normalmag == 0.0) 


{ 
} 


normalmag = 0.00001; /* a small number */ 


/* compute the light mag */ 
lightmag = sqrt((Ix*lx)+(ly *ly) +(lz*lz)); 
if(lightmag == 0.0) 


lightmag = 0.00001; /* a small number */ 
} 


/* compute N . L (normal dot product with the light source direction) */ 
dotprod = (normal|0] * lx) + (normal|t] * ly) + (normal{2! * Iz); 


dotprod = dotprod/(lightmag*normalmag); 
/* dotprod = cos(theta) of the angle between N and L. 
Convert to angle in radians */ 


radians = acos(dotprod); 


/* compute the color we should use */ 


if(-PIDIV2 <= radians && radians <= PIDIV2) 
{ 
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/* if the angle is negative, set to positive */ 
if(radians < 0.0) 


radians = -radians; 


} 


*colortouse = ((colormax-colormin) /PIDIV2)*(PIDIV2-radians)+colormin; 


} 


else 


{ 
} 


/* set the color */ 
color(*colortouse); 


*colortouse = colormin; 


/* draw the poly */ 
/* polf(ncoords,txyz); */ 
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LINE INTERSECT2 


#include "gl.h" 


#define X 0 

#define Z 2 

#define NONE 0 
#define INTERSECT 1 
#define PROPER = 2 


line intersect2(startl, end1, start2, end2, intersect, 
intersect type) 


float start1/3], end1|3’, start2/3], end2/3), intersect|3}; 
int “intersect type; 


{ 


/* given two lines of the form z = mx + b and z=nx+c, 
solving for x when the z’s are equal gives x = (c-b)/(m-n). 
Then solve for z using x in either of the above equations. */ 


float m,n,b,c; 
float minl x, min2 x, maxl x, max2 x, minl z, min2 z, maxl z, max2 z; 


*intersect type = PROPER; 


/* slope and z intercept of linel */ 


if (end1|X] != start1[X’) { 
m = (end1/Z) - start1{Z])/(end1,X! - start1|X)); 
b = ((start1/Z] - end1[{Z])/(end1[X] - start, X!)) * start1|X] + start1{Z]; 
if (end2/X] != start2/X}) { /* both lines are non-vertical */ 
/* slope and z intercept of line2 */ 


n = (end2 Z] - start2[Z])/(end2[X] - start2|X]); 





c = ((start2'Z) - end2{Z|)/(end2|X! - start2 X|)) * start2|X] + 
start2'Z]; 
te ne 
intersect|X] = (c-b)/(m-n); 
intersect|Z| = m*intersect|X] + b; 


else { /* both lines have equal slopes */ 
“intersect type = NONE; 
} 


else { /* linel is non-vertical, line2 is vertical */ 
intersect|X| = end2|X); 
intersect|Z| = m*intersect X 











+ b; 


else { 


184 


if (end2(X] != start2|X]) { /* linel is vertical, line2 is non-vertical* / 
/* slope and z intercept of line2 */ 
n = (end2/Z] - start2{Z])/(end2|X| - start2|X)); 
c = ((start2/Z] - end2[{Z|)/(end2|X] - start2/X])) * start2/X] + 
start2 Z'; 
intersect|X] = end1|X]; 
intersect|Z] = n*intersect|X] + c; 


else { /* both lines are vertical */ 
‘intersect type = NONE; 
} 
} 


if (*intersect type != NONE) { 
/* see if the intersection is proper, or if only the extensions of the 
line segments intersect */ 
if (start1[X] < end1|X]) { 
minl x = start1|X]; 
maxl x = end1/X]; 


else { 
minl x = end1|[X]; 
maxl x = start1|X]; 


} 

if (start1[Z| < end!{Z)) { 
minl z = start1[Z]; 
maxl z = end1[Z]; 

} 

else { 
minl z = end1'Z]; 
maxl z= startl Z); 


if (start2|X] < end2/X]) { 
min2 x = start2|X]; 
max2_ x = end2[X]; 

} 

else { 
min2 x = end2[X]; 
max2 x = start2[X]; 


} 

if (start2|Z] < end2[Z)) { 
min2 z = start2(Z]; 
max2 z = end2/Z]; 


else { 


min2 z = end2/Z]; 


max2 z = start2[Z]; 
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if ({intersect|/X] <= maxl x) && (intersect;X| <= max2 x) && 
(intersect|X] >= minl x) && (intersect|X] >= min2 x) && 
(intersect|Z] <= max z) && (intersect/Z] <= max2_z) & & 
(intersect|Z] >= minl_ z) && (intersect/Z] >= min2_2z)) { 


“intersect type = PROPER; 


else { 
*intersect type = INTERSECT; 
} 
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MAKENAVBOX 


/* drawnavbox.c - this function is called by the FOG-M missile simulator to 
build an object on top of the contour map in the upper right-hand corner 
of the screen. Navbox contains the direction arrow and view box in red. */ 


#include "gl.h" 
#include "fogm.h" 
#include "device.h" 


drawnavbox(navbox, arrowtag) 
Object *navbox; 
Tag ‘*arrowtag; 
{ 
*navbox = genobj(); /* create the navigation contol and display object */ 
makeobj(*navbox); 
if (TV) viewport(475,635,323,474); 
else viewport(768,1023,512,767); /* upper right hand corner of screen */ 
pushmatrix(); /* draw arrow in feet coordinates */ 
ortho2(-10.0,10.0 + NUMXGRIDS*FEETPERGRID, -10.0, 
-10.0 - NUMZGRIDS*FEETPERGRID); 
color(BLACK); 
clear(); 
color(128); 
*arrowtag = gentag(); 
maketag(*arrowtag); 
move2(0.0,0.0); 
draw2(0.0,0.0); 
draw2(0.0,0.0); 
move2(0.0,0.0); 
draw2(0.0, 0.0); 
rect(0.0,0.0,0.0,0.0); /* view box */ 
popmatrix(); 
closeob}(); 
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MAKEINDBOX 


/* makeindbox.c is a function that creates an object that displays the control 


idicators for the FOG-M missile simulation a 


#include "'gl.h" 
#include "fogm.h" 


makeindbox(indbox,headingtag,elevtag,altmsltag,speedtag,zoomtag,tilttag, pantag,desigtag) 
Object *indbox; 

Tag *headingtag, *elevtag, *speedtag, *zoomtag, *tilttag, *pantag, *desigtag; 

Tag *altmsltag; 


{ 


*indbox = genobj(); 

makeob)(*indbox); 

if (TV) viewport (475,635, 162,322); 

else viewport(768,1023,256,511); /* middle box on side of screen */ 
pushmatrix(); 

ortho2(0.0,255.0,0.0,255.0); /* use screen sized coordinates */ 


color(854); /* clear the window > 
clear(); 

linewidth(2); 

color(BLACK); 

recti(0,0,255,255); ,* outline box “| 
color(/ YELLOW), /* print labels for readouts */ 


cmov2i(10,240); 

charstr("SPEED"); 

cmov2i(55,225); 

charstr("kts"); 

cmov2i(90,240); 

charstr("HEADING"); 

cire( 140:0/232-0;3.0); /* "degree" symbol = 
cmov2i(180,240); 

charstr("Alt AGL"); /* AGL = above ground level */ 
cmov2i(225,225); 

charstr("ft"); 

cmov2i( 180,200); 

charstr("Alt MSL"); /* MSL = mean sea level */ 
cmov2i(225,185); 

charstr("ft"); 


cmov2i(50,130); 

charstr("ZOOM"); 

move2i(45,200); /* draw slider bar frame */ 
draw2i(25, 200); 

draw2i(25,70); 

draw2i(45,70); 

cmov2i(15,196); 
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charstr("8"); /* label slider bar values */ 
cmov2i(6,170); 
charstr("15"); 
cmov2i(6, 144); 
charstr("25"); 
cmov2i(6,118); 
charstr("35"); 
cmov2i(6,92); 
charstr("45"); 
cmov2i(6,66); 
charstr("55"); 


color( WHITE); /* readouts in white... a 
cmov2i(10,225); /* initialize to dummy values */ 
*speedtag = gentag(); 

maketag(*speedtag); 

charstr(" 200"); /* speed */ 


cmov2i({108,225); 

*headingtag = gentag(); 
maketag(*headingtag); 

ehansur|:) 0"); /* heading */ 


cmov2i(180,225); 

*elevtag = gentag(); 

maketag(*elevtag); 

charstr("1000"); /* altitude above ground level */ 


cmov2i(180,185); 
*altmsltag = gentag(); 


maketag(*altmsltag); 

charstr("1000"); /* altitude from mean sea level */ 
color(RED); 

*zoomtag= gentag(); /* indicator for zoom slider bar */ 
maketag(*zoomtag); 


move2(28.0,135.0); 
rdr2(10.0,5.0); 
rdr2(0.0,-10.0); 
rdr2(-10.0,5.0) 


° 


popmatrix(); 


if (TV) viewport (0,474,0,474); = /* reset for heads-up display */ 
else viewport (0,767,0,767); 


pushmatrix(); 
ortho2(0.0,767.0,0.0,767.0); /* use screen sized coordinates */ 
color(WHITE); 
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if (TV) linewidth(2); 
else linewidth(1); 


rect fi(365,370,370,375); /* draw center of crosshairs */ 
rect fi(396,370,401,375); 

rect fi(365,391,370,396): 

rect fi(396,391,401,396); 

move2i(0,383); 

draw2i(360,383); /* draw crosshairs sae 
move2i( 406,383): 

draw 21(767,383); 

move2i(383,0); 

draw 21(383,365); 

move2i(383,401); 

draw 2i(383,767); 


linewidth(2); 


move2i(30,50); /* draw TILT slider bar frame */ 
draw2i(40,50); 
draw2i(40,680); 
draw 2i(30,680); 
cmov2i(0,676); 
charstr("+25"); /* label slider bar values =) 
cmov2i(0,613); 
charstr("+20"); 
move2i(30,617); 
draw 2i(40,617); 
cmov2i(0,550): 
charstr("+15"); 
move2i(30,554); 
draw 2i(40,554); 
cmov2i(0,487); 
charstr("+10"); 
move2i(30,491); 
draw2i(40,491); 
cmov2i(0,424); 
charstr(" +5"): 
move2i(30,428); 
draw 2i(40,428); 
cmov2i(0,361); 
eharstr( 200). 
move2i(30,365); 
draw2i(40,365); 
cmov2i(0,298); 
charsth( -34): 
move 2i(30,302); 
draw 2i(40,302); 
cmov 2i(0,235); 
charstr("-10"); 
move2i(30,239); 
draw 2i(40,239); 
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cmov2i(0,172); 
charstr("-15"); 
move2i(30,176); 
draw 2i(40,176); 
cmov21(0, 109); 
charstr("'-20"); 
move2i(30,113); 
draw2i(40, 113); 
cmov2i(0,46); 
charstr("-25"); 


*tilttag = gentag(); 


maketag(*tilttag); 
move2(42.0,365.0); 


rdr2(10.0,-5.0); 
rdr2( 0.0,10.0); 
rdr2(-8.0,-4.0): 
rdr2( 6.0,-3.0); 
rdr2( 0.0, 4.0): 
rdr2(-2.0,-1.0); 
rdr2( 1.0,-1.0); 


move2i(120,15); 
draw2i( 120,25); 
draw 2i(750,25); 
draw 2i(750,15); 
cmov2i(107,3); 
charstr("-25"); 
cmov2i(170,3); 
charstr("-20"); 
move2i(183,15); 
draw 21(183,25); 
cmov2i(233,3); 
charstr("-15"); 
move2i(246,15); 
draw2i( 246,25); 
cmov2i(296,3); 
charstr("-10"); 
move 2i(309,15); 
draw2i(309,25); 
cmov2i(363,3); 
charstr("-5"); 
move 2i(372,15); 
draw 2i(372,25); 
cmov2i(431,3); 
charstr("0"); 
move2i(435,15); 
draw2i(435,25); 
cmov2i(494,3); 
charstr("+5"); 


move2i(498,15); 


draw 2i(498,25); 


/* draw PAN slider bar frame 


/* label slider bar values 
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ul 


/* indicator for TILT slider bar */ 


aa 


cmov2i(552,3); 

charstr("+ 10"); 
move2i(561,15); 
draw2i(561,25); 
cmov2i(615,3); 

charstr("+15"); 
move2i(624,15); 
draw2i(624,25); 
cmov2i(678,3); 

charstr("+20"); 
move2i(687,15); 
draw 2i(687,25); 
cmov2i(741,3); 

charstr("+25"); 


*pantag = gentag(); /* indicator for PAN slider bar */ 
maketag(*pantag); 

move2(435.0,27.0); 

Fan? (soe, 10-0). 

rdr2(- 10.05.00); 

rdr2( 4.0,-8.0); 

rdr2( 3.0, 6.0); 

rdr2( -4.0, 0.0); 

rdr2( 1.0,-2.0); 

rdr2( 1.0, 1.0); 


move2i(0,30); /* designate/reject box */ 
draw 2i( 100,30); 

draw2i(100,0); 

*desigtag = gentag(); 

maketag(*desigtag); 

cmov2i(10,10); 

charstr("DESIGNAT EE"): 


popmatrix(); 
closeobj(); 
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MAKEINSTRBOX 


/* makeinstrbox.c - this function builds an object that contains an instruction 


summary for the FOG-M missile simulation ai 


#include "gl.h" 
#include "fogm.h" 


makeinstrbox(instrbox) 


Object *instrbox; 


{ 


*instrbox = genobj(); 

makeobj(*instrbox); 

if (TV) viewport (475,635,0,161); 

else viewport(768,1023,0,255); = /* box is in lower right hand corner */ 
pushmatrix(); 


ortho2(0.0,255.0,0.0,255.0); /* use screen-sized coordinates ah 
color(851); /* use a medium green ah 
clear(); 

linewidth{2); 

color(852); /* use light brown af 


rectfi(10,20,110,195); /* draw the mouse control box */ 
rectfi(135,80,245,195); /* draw the dial control box */ 
color{(BLACK); /* outline controls os) 
recti(10,20,110,195); 

recti({135,80,245,195); 

rect1(0,0,255,255); 


color(BLACK); 

cmov 21{60,230); 

charstr © ON TROLS"); 
cmov2i{37,200); 
charstr("MOUSE"); 

cmov2i( 172,200); 
charstr("DIALS"); 


cmov2i(25,60); 

charstr("TILT"); 

move2i(70,62); /* draw arrow */ 
draw2i(75,55); 

draw2i(75,75); 

draw2i(70,68); 

move2i(75,75); 

draw2i(80,68); 

move2i(75,55); 

draw2i(80,62); 
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cmov2i(25,30); 

charstr("PAN"); 

move2i(67 ,40); /* draw arrow */ 
draw2i(60,35); 

draw2i(80,35); 

draw2i(73,40); 

move2i(80,35); 

draw2i(73,30); 

move2i(60,35); 


draw2i(67,30); 
color(853); /* dark brown a 
rect fi(20,85,40,185); /* draw mouse buttons a 


rect fi(50,85,70,185); 

rectfi(80,85, 100,185); 

color(BLACK), /* outline buttons oy) 
recti(20,85,40,185); 

recti(50,85,70,185); 

recti(80,85, 100,185); 


color(853); 

circfi( 160,165 ,20); /* draw dials a) 
circfi(160,110,20); 

circ fi(220,165,20); 

circfi(220,110,20); 

color(BLACK); /* outline dials 
circi(160,165,20); 

circi(160,110,20); 

circi(220, 165,20); 

circi(220,110,20); 

color( WHITE); 

cmov2i(147,160); 

charstr("SPD"); /* label dials */ 

cmov2i( 147,106); 

charstr("DIR"); 

cmov 2i(207, 106); 

charstr("ALT"); 

cemov 2i(207,160); 

charstr("CLR"); 


cmov2i(25,170); 
charstr("Z"); —/* label mouse buttons */ 
cmov2i(25,158); 
chiarstn (307). 
cmov2i(25, 146); 
charstr("O"); 
cmov2i(25,134); 
charstr("M"); 
cmov2i(25,110); 
cliarstri le); 
cmov2i(25,98); 
charstr("N"); 
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cmov2i(55,170); 
charstr("D"); 
emov2i(55, 158); 
charstr("E"); 
cmov2i(55, 146); 
Onarstry (5°): 
cmov2i(55, 134); 
charstr("I"); 
cmov2i(55,122); 
eharstr("G"); 
cmov2i(85,170); 
charstr("Z"); 
cmov21(85,158); 
charstr("O"); 
cmov2i(85, 146); 
ehlarsen( 'O""); 
cmov 21(85,134); 
charstr("M"); 
cmov2i(85,110); 
enarstr("'O'); 
cmov 2i(85,98); 
charstr("U"); 
cmov2i(85,86); 
charstr("T"); 


popmatrix(); 
closeobj(); 
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MAKEMAP 


/* makemap.c - this function is called by the FOG-M missile simulator to 


build an object containing a contour map. The map 1s used for the full 
screen display in prelaunch, and in the upper right corner of the flight 
display in fogm. a 


#include "gl.h" 
#include "fogm.h" 
#include "device.h" 


makemap(contour) 
Object *contour; 


{ 


short i, j, elev, length, lastcolor, breakpt| 15]; 
int colour; 
extern short gridpixel 100 |100]; /* terrain elevations & vegetation * / 


/* compute elevations where color changes should occur */ 


for (i = 1; 1 < 16; i++) breakptli-1] = (((MAX - MIN) / 16 ) * i) + MIN; 


*contour = genobj(); /* create the navigation contol and display object */ 
makeobj(*contour); 

view port(0,767,0,767); 

pushmatrix(); 

ortho2(0.0,100.0,0.0,100.0);  /* use array index space */ 

color(BLACK); 


clear(); 


lastcolor = BLACK; 
linewidth(8); 


for (i=0; 1 < 100; ++3) { /* draw column i p 
move 2i(i,0); /* start at bottom of column = */ 
length = 0; /* # adjacent points of the same color */ 
for (j = 0;j < 100; ++)) { /*foreachrowincolumni  */ 
elev = gridpixel|j||i] & elev mask; /* mask off veg code */ 
if (elev < breakpt|O]) colour = 16; /* assign green colors 7 
else if (elev < breakpt({1]) colour = 17; 
else if (elev < breakpt|2]) colour = 18; 
else if (elev < breakpt|3]) colour = 19; 
else if (elev < breakpt|4|) colour = 20; 
else if (elev < breakpt|5|) colour = 21; 
else if (elev < breakpt|6]) colour = 22; 
else if (elev < breakpt|7|) colour = 23; 
else if (elev < breakpt|8]) colour = 24; 
else if (elev < breakpt|9|) colour = 25; 
else if (elev < breakpt({10]) colour = 26; 
else if (elev < breakpt{11]) colour = 27; 
else if (elev < breakpt{12]) colour = 28; 
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else if (elev < breakpt 13]) colour = 29; 
else if (elev < breakpt 14]) colour = 30; 
else colour = 31; 


/* if veg-code = 0 (i.e. veg < 1 meter) shift to brown colors */ 
if (!((gridpixel|j|/i] >> 13) & veg mask)) colour += 16; 


if (colour == lastcolor) length++; /* don’t draw yet */ 


else { /* draw now that color has changed */ 
color(lastcolor); 
rdr2i(0,length); 
lastcolor = colour; /* reset for new draw */ 
length = I; 


\ /* end for j */ 


color(colour); /* draw last (top) line */ 
rdr2i(0,length); 
\ /* end fori */ 


if ('TV) { 
color(BLACK); /* draw grid on top of map */ 
linewidth(1); 
for (i = 10; 1 < 100; i+=10) { /* draw interior lines * / 
move2i(i,0); /* horizontals * / 
draw 2i(i, 100); 
move2i(0,i); /* verticals */ 


draw2i(100,1); 


} 


linewidth(2); /* draw exterior border */ 
rect (0.0,0.0,100.0,100.0); 


popmatrix(); 
closeobj(); 
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MAKESCREENS 


/* makescreens.c - builds graphical objects for prelaunch’s instructional 


screens and readout boxes. a 


#include "gl.h" 
#include "device.h" 
#include "fogm.h" 


makescreens(obj,tag) 


Object obj|7]; 


Tag 


{ 


tag 6}; 


obj|INSTR] = genobj(); /* object for pre-launch instructions */ 
makeobj(obj INSTR)); 

if (TV) viewport(475,635,239,474);: 
else viewport (767 ,1023,385,767); 
pushmatrix(); 
ortho2(0.0,255.0,0.0,384.0); 
color(CYAN); 

clear(); 

color(BLUE); 

rectfi(10, 10,245,374); 

color(WHITE); 

cmov2i(30,340); 
charstr("PRE-LAUNCH INSTRUCTIONS"); 
cmov2i(25,300); 

charstr("1. PRESS LEFT MOUSE"); 
cmov2i(52,285); 

charstr(" BUTTON TO LOCK IN"); 
cmov2i(52,270); 

charstr("LAUNCH POSITION"); 
cmov 2i(25,220); 

charstr("'2. PRESS RIGHT MOUSE"); 
cmov2i(52,205); 

charstr("BUTTON TO LOCK IN"); 
cmov 2i(52,190); 

charstr("TARGET LOCATION"); 
cmov2i(25,140); 

charstr("3. PRESS MIDDLE MOUSE"); 
emov 21/52 02a): 

charstr("BUTTON TO LAUNCH"); 
cmov2i(25, 75); 

charstr("'4. PRESS ALL THREE"); 
cmov2i(52, 60); 

charstr("BUTTONS [TO Exit: 
popmatrix(); 

closeobj(); 
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/* define object for displaying user input for missile launch 
position and target location. Also displays computed heading 
and distance to target */ 


obj[STATS] = genobj(); 
makeobj(obj[/STATS]); 

if (TV) viewport(475,635,0, 238); 
else viewport(767,1023,0,384); 
pushmatrix(); 
ortho2(0.0,255.0,0.0,384.0); 
color(CYAN); 

clear(); 

color(BLUE); 
rectfi(10,10,245,374); 

color( WHITE); 
cmov2i(30,340); 
charstr("PRE-LAUNCH STATISTICS"); 
cmov2i(25,260); 
charstr("LAUNCH POSITION: 10SFQ"); 
cmov2i(70,235); 

charstr("X COORD: "); 
cmov2i(70,220); 

charstr("Y COORD: "); 
cmov2i(170,235}; 

tag) LAUNCH] = gentag(); 
maketag(tagj/LAUNCH]); 
charstr(" "); 
cmov2i(170,220); 

charstr(" "); 

cmov2i(25, 180); 
charstr("TARGET LOCATION: 10SFQ"); 
cmov2i(70,155); 

charstr("X COORD: ny 
cmov2i(70,140); 

charstr("Y COORD: ui 
cmov2i(170,155); 

tagi| TARGET] = gentag(); 
maketag(tag|TARGET)); 
charstr(" "); 
cmov2i(170,140); 

charstr("' '); 
cmov2i(25,100); 
charstr("HEADING: "'}; 
cmov2i{ 25,60); 
charstr("DISTANCE: "); 
cmov2i( 106,100); 

tag/ HEAD] = gentag(); 
maketag(tagjHEAD)); 
charstr(" "'); 

cmov2i(115,60); 

eharstr(" "): 

popmatrix(); 
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closeobj(}; 
/* define object for lines & circles showing flightpath on contour map */ 


obj{[FLTPATH] = genobj(); 
makeobj(obj/|FLTPATH}); 
pushmatrix(); 

if (TV) viewport (0,474,0,474); 
else viewport (0,767 ,0,767); 
ortho2(0.0,100.0,0.0,100.0); 
color(BLACK); 

clear(); 

color(64); 

linewidth(3); 

tag|MISSILE] = gentag(); 
maketag(tag| MISSILE]); 
circf(0.0,0.0,0.0); 
move2(0.0, 0.0, 0.0); 
draw2(0.0, 0.0. 0.0); 
color(128); 

tagiTGT] = gentag(); 
maketag(tag|TGT)); 
circf(0.0,0.0,0.0); ; 
popmatrix(}; 

closeobj(); 


/* define object for displaying first screen of operator instructions */ 


obj[SCREEN1] = genobj(); 
makeobj(obj[SCREEN1)); 

color(BLUE); /* set background color */ 
clear(); 

color(RED); 

linewidth(10); 

recti(0,0,1023,767); 

linewidth(1); 

color(WHITE); 

cmov2i(420,500); 

charstr("WELCOME"); 

cmov2i(420,450); 

charstr("TO THE"); 

cmov2i(320,400); 

charstr("FIBER-OPTICALLY GUIDED MISSILE"); 
cmov 21{420,350}; 

charstr("(FOG-M)"); 

cmov2i(410,300); 

charstr("SIMULATION"); 

cmov2i(310,100); 

charstr("PRESS MIDDLE MOUSE BUTTON TO CONTINUE..."), 
cmov2i(315,85); 

charstr("OR PRESS ALL 3 MOUSE BUTTONS TO EXIT."); 
closeobj(); 
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/* define object for displaying second screen of operator instructions */ 


obj/[SCREEN2] = genobj(); 

makeobj(obj|[SCREEN2)); 

color(BLUE); /* set background color */ 

clear(); 

linewidth(10); 

color(RED); 

recti(0,0,1023,767); 

linewidth(1); 

color( WHITE); 

cmov2i(210,600); 

charstr("THE FOG-M PROGRAM PROVIDES A SIMULATED MISSILE LAUNCH AND"); 
cmov2i(210,575); 

charstr("OUT-THE-WINDOW VIEW OF THE TERRAIN AS SEEN FROM THE OPERATOR’S"); 
cmov2i(210,550); 

charstr("CONSOLE ON THE GROUND."); 

cmov2i(210,500); 

charstr("THE GENERAL AREA FOR THIS FLIGHT SIMULATION IS FT HUNTER LIGGETT"}; 
cmov2i(210,475); 

charstr("CALIFORNIA AND VICINITY."); 

cmov2i(210,425); 

charstr("THE SPECIFIC TEST AREA IS A 10 KILOMETER REGION DESIGNATED BY"); 
cmov 21(210,400); 

charstr("UNIVERSAL TRANSVERSE MERCATOR (UTM) GRID COORDINATES 10SFQ58."); 
cmov2i(300,100); 

charstr("PRESS MIDDLE MOUSE BUTTON TO CONTINUE,"); 

cmov2i(305,85); 

charstr("OR PRESS ALL 3 MOUSE BUTTONS TO EXIT."); 

closeobj(); 


/* define object for displaying third screen of operator instructions */ 


obj[SCREEN3] = genobj(); 

makeobj(obj[SCREEN3)); 

color(BLUE); /* set background color */ 

clear(); 

linewidth(10); 

color(RED); 

recti(0,0,1023,767); 

linewidth(1); 

color( WHITE); 

cmov 21{285,650); 

charstr("PRE-LAUNCH ORIENTATION"); 

cmov 21( 200,600); 

charstr("1. WHEN THE PRE-LAUNCH PHASE OF THE FOG-M SIMULATION BEGINS, A"); 
cmov2i(200,585); 

charstr("2-DIMENSIONAL CONTOUR MAP OF THE TEST AREA (UTM 10SFQ58) WILL BE"); 
cmov2i(200,570); | 

charstr("DISPLAYED ON THE OPERATOR CONSOLE. TWO CONTROL PANELS CONTAINING"): 
cmov2i(200,555); 

charstr("PRE-LAUNCH INSTRUCTIONS AND CURRENT LAUNCH STATISTICS WILL ALSO"); 


201 


cmov2i(200,540); 

charstr("BE DISPLAYED."); 

cmov2i(200,490); 

charstr("2. THE OPERATOR WILL BE REQUIRED TO PROVIDE TWO CRITICAL DATA"); 
cmov2i(200,475); 

charstr("ITEMS TO THE LAUNCH CONTROL SYSTEM; INITIAL LAUNCH POSITION AND"); 
cmov2i(200,460); 

charstr("TARGET LOCATION."); 

cmov2i(200,410); 

charstr("3. TO DEFINE INITIAL LAUNCH POSITION, MOVE CURSOR OVER DESIRED"); 
cmov2i(200,395); 

charstr("LOCATION (REFER TO LAUNCH STATISTICS CONTROL PANEL TO VIEW THE"); 
cmov2i( 200,380); 

charstr('CURRENT UTM GRID COORDINATES). PRESS LEFT MOUSE BUTTON TO LOCK"); 
cmov2i(200,365); 

charstr("IN LAUNCH POSITION."); 

cmov2i(200,315); 

charstr("4. TO DEFINE TARGET LOCATION, MOVE CURSOR OVER DESIRED LOCATION"); 
cmov2i(200,300); 

charstr("(REFER TO LAUNCH STATISTICS CONTROL PANEL TO VIEW CURRENT UTM"); 
cemov2i( 200,285): 

charstr("GRID COORDINATES). PRESS RIGHT MOUSE BUTTON TO LOCK IN TARGET"); 
cmov2i(200.270); 

charstr("LOCATION. THE BLUE LINE DISPLAYS THE PROJECTED FLIGHT PATH. THE"); 
cmov2i( 200,255); 
charstr("MISSILE WILL FLY AT A CONSTANT VELOCITY AND HEADING. THE LAUNCH"); 
cmov2i(200,240); 
charstr("STATISTICS CONTROL PANEL WILL DISPLAY COMPUTED MISSILE HEADING"); 
cmov2i(200,225); 

charstr("IN DEGREES (0 DEGREES DUE NORTH)."); 

cmov2i(240,100); 

charstr("PRESS MIDDLE MOUSE BUTTON TO MOVE INTO PRE-LAUNCH PHASE,"); 
cmov2i(326,85); 

charstr("OR PRESS ALL 3 MOUSE BUTTONS TO EXIT."); 

closeobj(); 
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MAKETANK 


#include "gl.h" 
#include "fogm.h" 


maketank(item) 
Object *item; 


{ 


long points = 4, bigpoints = 8; 

float parray[8]|3]; 

float Ix,ly,1z; 

long cmin = MIN _TGT COLOR, cmax = MAX _TGT COLOR, cl; 


Ix = 400.0 * 41.01; /* direction of lightsource */ 
ly = 6000.0; 
lz = 200.0 * (-41.01); 


*item=genobj(}; 
makeobj(*item); 


/* draw right side of tank CCW */ 
parray(0|[0] = -10.0; 
parray|0|[1] = 6.0; 
parray|0|[2] = -5.0; 
parray(1|[0] = -15.0; 
parray|1]{1] = 4.0; 
parray/|1|{2] = -5.0; 
parray|2|[0] = -15.0; 
parray|2|[1] = 2.0; 
parray|2]{2] = -5.0; 
parray(3]{0] = -10.0; 
parray|3|{1] = 0.0; 
parray|3]{2] = -5.0; 
parray|4|{0] = 10.0; 
parray|4|{1] = 0.0; 
parray|4|{2] = -5.0; 
parray(5|(0| = 15.0: 
parray|5|[1] = 2.0; 
Barray|5||2| = =5x0; 
parray|6|/0| = 15.0; 
parray(6|[1] = 4.0; 
parray|6|{2] = -5.0; 
parray|7|(0| = 10.0; 
parray|7|[1] = 6.0; 
parray|7|[2] = -5.0; 
lightorient(parray,bigpoints,0.0,0.0,0.0,1x,ly,lz,cmin,cmax,&cl); 
color(c1); 
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polf(bigpoints, parray); 


/* front of tank CW */ 
parray(0][0] = 15.0; 
parray(0]|1| = 5.0; 
parray|0|[2] = -5.0; 
parray|1](0] = 15.0; 
parray({1]/1} = 3.0; 
parray(1|{2] = -5.0; 
parray(2]/0] = 15.0; 
parray|2][1] = 3.0; 
parray|2|{2] = 5.0; 
parray(3][0} = 15.0; 
parray|3]{1] = 5.0; 
parray(3|(2| — 5.0; 
lightorient(parray,points,0.0,0.0,0.0,lx,ly ,lz,cmin,cmax,&cl1); 
color(c1); 
polf(points,parray); 


/* draw left side of tank CW */ 
parray(0|[0! = 10.0; 
parray|0][1] = 6.0; 
parray|0|[2] = 5.0; 
parray|1[0| = 15.0; 
parray|1/[1] = 4.0; 
parray(1][2! = 5.0; 
parray(2|/0] = 15.0; 
parray|2|[1] = 2.0; 
parray|2]|2| = 5.0; 
parray|3]/0] = 10.0; 
parray|3]{1] = 0.0; 
parray(3]/2] = 5.0; 
parray|4][0] = -10.0; 
parray(4][1| = 0.0; 
parray|4]{2] = 5.0; 
parray|5]|(0| = -15.0; 
parray(5|[1] = 2.0; 
parray|5|[2} = 5.0; 
parray|6][0) = -15.0; 
parray(6|[1] = 4.0; 
parray|6][2] = 5.0; 
parray!7]|{O] = -10.0: 
parray|7|[1| = 6.0; 
parray|7|/2} = 3.0: 
lightorient{ parray, big points,U0.0,0.0,0.0,1x,ly,1z,cmin,cmax,&c1); 
color(c1); 
polf(bigpoints, parray); 


/* back of tank CCW */ 
parray(0][0] = -15.0; 
parray|[O]{1] = 4.0; 
parray(0]|[2| = 5.0; 
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parray(1]/0] = 
parray(1]/1] = 
parray(1](2] = 
parray([2][0] = 
parray/2][1] = 
parray[2][2] = 


=~ 13.0; 
2.0: 
9.0; 
-15.0; 
Zee 


-5.0; 


parray|3](0] = -15.0; 
parray[3|[1] = 4.0; 


parray(3][2] = -5.0; 
lightorient(parray,points,0.0,0.0,0.0,1x,ly,lz,cmin,cmax,&c1); 
color(cl); 


polf(points,parray); 


/* top middle of tank body CCW */ 
parray(0|(0] = -10.0; 
parray(0]/1] = 6.0; 


parray(0]/2] = -5.0; 
parray(1][0] = -10.0; 
parray(1][1] = 6.0; 
parray|1]{2] = 5.0; 
parray|2][0] = 10.0; 
parray(2][/1] = 6.0; 
parray(2][2] = 5.0; 
parray({3|/0| = 10.0; 
parray|3|[1] = 6.0; 


patcay|3||2| — -5.0; 
lightorient(parray,points,0.0,0.0,0.0,1x,ly,lz,cmin,cmax,&c1); 
color(c1); 

polf(points,parray); 


/* top front of tank body CCW */ 
parray(0]/0] = 10.0; 
parray(0][1] = 6.0; 
parray(0][2] = -5.0; 
parray(1][0] = 10.0; 
parray|1|{1] = 6.0; 
parray(1][2] = 5.0; 
parray/|2]{0] = 15.0; 
parray|2][1] = 4.0; 
parray(2]/2] = 5.0; 
parray(3]/0] = 15.0; 


parray(s|(1] = 4.0: 
parray|3|{2| = -5.0; 


lightorient(parray, points,0.0,0.0,0.0,1x,ly,lz,cmin.cmax,&c1); 
color(cl); 
polf(points, parray); 


/* top back of tank body CCW */ 
parray(0][0] = -10.0; 

parray(0|[1] = 6.0; 

parray(0][2] = -5.0; 

parray|1][0] = -15.0; 
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parray(|1]{1] = 

parray([1]{2] = -5.0; 

parray|2|(0] = -15.0; 

parray|2][1] = 4.0; 

parray[2|[2] = 

parray(3][0] = -10.0; 

parray(3|[1] = 6.0; 

parray(3][2] = 5.0; 
lightorient(parray,points,0.0,0.0,0.0,lx,ly,lz,cmin,cmax, &c1); 
color(c1); 

polf(points,parray); . 


/* bottom middle of tank CW*/ 
parray|0][0] = -10.0; 

parray|0][1] = 0.0; 

parray|0|{2] = -5.0; 

parray|1]{0] = 10.0; 
parray|1|[1] = 0.0; 

parray|1]{2] = -5.0; 

parray|2]/0] = 10.0; 
parray(2]/1] = 0.0; 

parray(2]{2] = 5.0; 

parray|3 Ol == 1000; 
parray|3|[1] = 0.0; 
parray|3]|2] = 5.0; 
lightorient(parray ,points,0.0,0.0,0.0,1x,ly,lz,cmin,cmax,&c1); 
color(c1); 
polf(points, parray ); 


/* bottom front of tank CW */ 
parray(0||0] = 10.0; 
parray(0|{1] = 0.0; 
parray|0|[2] = -5.0; 
parray(1]{0] = 15.0; 
parray|1|{1}] = 2.0; 
parray|1|{2] = -5.0; 
parray|2]|0] = 15.0, 
parray(2]{1] = 2.0; 
parray|2|[2] = 5.0; 
parray|3](0] = 10.0; 
parray(3][1] = 0.0; 


parrav| sic) =— 20: 
lightorient(parray, points.0.0,0.0,0.0.lx.ly,Jz,cmin,cmax.&cl)}; 
color(cl); 


polf( points, parray ); 


/* bottom back of tank CW */ 
parray|0|[0] = -10.0; 
parray|0]|1| = 0.0; 

parray|0][2] = -5.0; 
parray|1](0] = -10.0; 
parray(1][1] = 0.0; 
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parray|1]|2| = 5.0; 
parray|2|{0| = -15.0; 
parray|2][1] = 2.0; 
parray|2||2] = 5.0; 


parray|3][0] = -15.0; 

parray(3][1] = 2.0; 

parray|3|(2] = -5.0; 

lightorient(parray, points,0.0,0.0,0.0,lx,ly,lz,cmin,cmax,&cl); 
color(c1); 


polf(points, parray); 


/* right side of gun barrel */ 
parray|0||0] = 1.6667; 
parray|O|{1] = 8.0; 
parray|0|{2] = -0.5; 
parray|1]/0] = 2.3333; 
parray(1|[1] = 7.0; 
parray|1|[2] = -0.5; 
parray|2//0] = 17.0; 
parray(2]/1] = 7.0; 
parray(2]/2| = -0.5; 
parray(3]{0] = 17.0; 
parray|3][1] = 8.0; 
parray/3|[2) = -0.5; 
lightorient(parray,points,5.0,2.5,0.0,1lx,ly,1z,cmin,cmax, &c1); 
color(cl); 

polf(points, parray); 
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/* top of gun barrel */ 
parray(0||0] = 1.6667; 
parray(0][1] = 8.0 
parray|0][2] = 0.5; 
parray(1]/0] = 1.6667; 
parray(1]{1] = 8.0; 
parray(1]/2] = -0.5; 
parray|2|(0] = 17.0; 
parray(2|[1] = 8.0; 
parray(2]||2] = -0.5; 
parray|3][0] = 17.0; 
parray|3]{1] = 8.0; 
parray|3]/2] = 0.5; 
lightorient(parray.points,5.0,2.5,0.0.1x,!y,1z,cmin,emax,&c1); 
color(ci); 

polf( points, parray}; 


/* left side of gun barrel */ 
parray|0][0] = 17.0; 
parray(O|[1] = 8.0; 
parray|0]/2] = 0.5; 
parray(1]/0] = 17.0; 
parray|1]/1] = 7.0; 
parray|1]/2] = 0.5; 
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parray|2|(0] = 2.3333; 

parray(2][1] = 7.0; 

parray|2|[2| = 0.5; 

parray|3]/0] = 1.6667; 

parray|3][1] = 8.0; 

parray(3][2] = 0.5; 
lightorient(parray,points.5.0,2.5,0.0.1x,ly,lz,cmin,cmax, &cl); 
color(cl); 

polf(points, parray); 


/* end of gun barrel */ 
parray |0|[0] = 17.0; 
parray|O][1] = 8.0; 
parray|0][2] = -0.5; 
parray(1]/0] = 17.0; 
parray|1|/1] = 7.0; 
parray|1|(2] = -0.5; 
parray|2]/0] = 17.0; 
parray|2][1! = 7.0; 
parray|2|(2| =* O25. 
parray[3]|0| = 17.0; 
parray(3]/1] = 8.0; 
parray/3|[2| = 0.5; 
lightortent( parray, points,5.0,2.5,0.0,!x,ly,1z,cmin,cmax,&cl); 
color(cl); 

polf( points, parray); 
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/*bottom of gun barrel */ 

parray|0][0] = 2.3333; 

parray|0|[1] = 7.0; 

parray|0|[2] = 0.5; 

parray|1](0] = 2.3333; 

parray(1][1] = 7.0; 

parray(1]{2] = -0.5; 

parray(2]|0] = 17.0; 

parray(2][1] = 7.0; 

parray(2]/2] = -0.5; 

parray(3]{0] = 17.0; 

parray|3][1] = 7.0; 

parray|3][2] = 0.5; 
lightorient(parray,points,5.0,2.5,0.0,lx,ly,lz,cmin,cmax,&cl); 
color(c1); 

ooif( points, parray); 

/* right side of turret * 
parray(0][0] = -3.0; 
parray(0]/1] = 9.0; 
parray|0|/2] = -1.0; 
parray(|1]/0] = -5.0; 
parray|1][{1] = 6.0; 
parray|1]/2| = 220; 
parray(2|/0] = 3.0; 
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parray|2|[1] = 6.0; 
parray|2|[2] = -3.0; 
parray|3|{0] = 1.0; 
parray([3][1] = 9.0; 
parray(3][2] = -1.0; 


lightorient(parray,points,-1.0,2.5,0.0,lx,ly,lz,cmin,cmax, &c1); 
color(c1); 
polf(points, parray); 


/* front side of turret */ 
parray|0|(0] = 1.6667; 
parray|0]/1] = 9.0; 
parray|0|[2] = -1.0; 
parray|1][0] = 3.0; 
parray(1|[1] = 6.0; 
parray(1][2] = -3.0; 
parray|2]{0] = 3.0; 
parray|2](1] = 6.0; 
parray|2|(2] = 3.0; 
parray|3|(0] = 1.6667; 
parray(|3|[1] = 9.0; 
parray|3][2] = 1.0; 
light orient (parray,points,-1.0,2.5,0.0,1x,ly,lz,cmin,cmax,&c1); 
color(c1); 

polf(points, parray); 


/* left side of turret */ 


parray|0|/0] = 1.6667; 
parray|0|[1] = 9.0; 
parray|0|[2] = 1.0; 
parray|1]/0] = 3.0; 
parray|1]/1] = 6.0; 
parray|1|{2] = 3.0; 


parray|2](0] = -5.0; 

parray|2|[1] = 6.0; 

parray(2][2] = 3.0; 

parray(3]|0] = -3.0; 

parray(3][1] = 9.0; 

parray|[3](2] = 1.0; 

lightorient(parray ,points,-1.0,2.5,0.0,1x,ly ,lz,cmin,cmax,&c1); 
color(cl); 

polf( points, parray); 
’* back side of turret */ 
parray|0|/0| = -3.0; 
parray|0|{1] = 9.0; 
parray|0|[2] = 1.0; 
parray|1]{0] = -5.0; 
parray/1|{1] = 6.0; 
parray|1]/2] = 3.0; 
parray|2]/0| = -5.0; 
parray|2|{1] = 6.0; 
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parray|2|{2| = -3.0; 

parray(3|(0] = -3.0; 

parray{3]{1] = 9.0; 

parray(3]{2] = -1.0; 
lightorient{parray,points,-1.0,2.5,0.0,lx,ly,lz,cmin,cmax,&c1); 
color({c1); 

polf(points, parray); 


/* top of turret */ 
parray(0|[0] = -3.0; 
parray(0|{1] = 9.0; 
parray|0]{2} = 1.0; 
parray{1|{0| = -3.0; 
parray|1][1] = 9.0; 
parray|1|(2] = -1.0; 
parray|2]|0] = 1.0; 
parray(2][1] = 9.0; 
parray|2|(2] = -1.0; 
parray(3][0] = 1.0; 
parray([3|{1] = 9.0; 
parray|3|{2] = 1.0; 
lightorient(parray,points,-1.0,2.5,0.0,lx,ly,lz,cmin,cmax,&c1); 
color(c1); 
polf(points, parray ); 


closeobj{); 


210 


NEAREST TGT 


#include "gl.h" 

#include "fogm.h" 

nearest tgt(vx,vy,vz,px,py,pz,tgt idx) 
Coord vx, vy, VZ, PX, Py, Pz; 

int *tgt idx; 


{ 
float dist, dist_to_los(); 
float min dist; 
float num_tgts; 
extern float tgt_ pos|MAX TGTS][3}; 


int index; 


num tgts = 10; 
min dist = dist to los(vx,vy,vz,px,py,pz,étgt_pos|0|[0}); 
*tgt_ idx = 0; 


for (index = 1; index < num_tgts; ++index) { 
dist = dist to los(vx,vy,vz,px,py,pz,étgt pos|index]{0]); 
if (dist < min dist) { 
min dist = dist; 
“tgt_idx = index; 
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NPOLY ORIENT 


/* npoly orient.c */ 


#include <gl.h> 
+#include <math.h> 


int npoly_orient(ncoords,xyz,xinside,yinside,zinside) 
unsigned int ncoords; 

Coord xyz||{3]; 

Coord xinside, yinside, zinside; 


{ 


register unsigned short int i,j; /* loop temps */ 
Coord center(3];_ /* center coordinate of the polygon */ 


Coord a[3], b[3]; /* vector hold locations for the vectors that run 
from the center coordinate to the points of the 
polygon */ 


Coord xn[3], xmn[3];_ /* points on line containing normal that are 
on opposite sides of the piane containing 
the polygon. 


vy 


float distton; /* distance to point n from pt inside. */ 
float disttomn; /* distance to point -n from pt inside. */ 
Coord normal(3}; /* the normal vector computed from a x b */ 


/* compute the center coordinate of the polygon */ 
center|0] = 0.0; 
center|1] = 0.0; 
center|2] = 0.0; 


for(i=0; i < ncoords; i++) 
{ 
for(j=0; j < 3; j++) 


{ = 
\ 


I 


center|j| += xyziil/j|; 


} 


/* divide out by the number of coordinates */ 
for(j=0; j < 3; j++) 


center|j] = center|j|/(float)ncoords; 
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/* check the first 2 coordinates of the polygon for their direction */ 


/* compute vector a. It runs from the center coordinate to coordinate 0 */ 
for(j=0; j < 3; j++) 


alj] = xyz[0]|j] - center|j]; 


/* compute vector b. It runs from the center coordinate to coordinate 1 */ 
for(j=0; j <3; j++) 


b{j] = xyz[1][j] - center[}]; 


/* compute a x b to get the normal vector */ 
normal|[0] = a{1]*b/2] - a{2]*b/1]; 
normal|1] = a|2|*b{0} - al0|*bl2) 
normal[2] = a[0]*b(1] - a(1]*b(0) 


/* compute point n, offset pt from center in direction of normal */ 
for(j=0; } < 3; j++) 


xn|j] = center/j] + normallj]; 


/* compute point -n, offset pt from center in opposite direction 
from normal. 
* 


for(j=0; j < 3; j++) 


xmn|j] = center|j] - normal|j]; 


/* compute the distance the inside pt is from point n */ 
distton = sqrt((xn[0] - xinside) * (xn(O] - xinside) + 
(xn|1] - yinside) * (xn[1] - yinside) + 
(xn{2] - zinside) * (xn[2] - zinside)); 


/* compute the distance the inside pt is from point -n */ 
disttomn = sqrt((xmn|0] - xinside) * (xmn{(0] - xinside) + 
(xmn|[1] - yinside) * (xmn{1] - yinside) + 
(xmn/{2! - zinside) * (xmn/{2] - zinside)); 
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/* if the dist(n) < dist(-n), then n points back towards the 
inside point and is on the same side of the plane as inside. 
a x b is then clockwise. 


* 
if(distton < disttomn) 
return(1); /* clockwise */ 


else 


return(0); /* counterclockwise */ 
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PRELAUNCH 


/* The function prelaunch is the user interface portion of the FOG-M 
flight simulation. It allows the operator to interactively enter 
critical data items necessary to simulate the missile in flight. 

The function returns the initial launch position in the x-z plane 
and also the direction of flight. */ 


#include "gl.h" 
#include "device.h" 
#include "fogm.h" 
#include "math.h" 


prelaunch(vx, vy, vz, direction, compassdir, active, obj, tag) 
Coord *vx, *vy, *vz; 

double *direction; 

float *compassdir; 

int *active; 

Object obj{7]; 

Tag tag(6l; 


{ 


float gnd_level(); 

float compass(); 

int screencnt, launchlock, targetlock; 

int xval, yval, xlaunch, ylaunch, xtarget, ytarget, utm x, utm jy; 
char xtemp[35], ytemp([35], dist{35], heading[35]; a 7 
float distance; 

double xdistance, ydistance; 

Colorindex unmask; 


xtemp|0| =’ ’; 
ytemp|0| =’ ’; 
aist}0| =” ’; 


heading|0| =’ ’; 


unmask=(1<<getplanes()) -1; 
writemask (unmask); 

af (TV) viewport(0,635.0,474); 
else viewport(0,1023,0,767)}; 
pushmatrix()}; 
ortho2{0.0,1023.0,0.0,767.0); 


*direction = 0.0; /* initialize the direction */ 
cursoff(); /* turn the cursor off */ 
callobj(obj[SCREEN1)); /* display screen 1 */ 
swapbuffers(); 
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screencnt = 1; /* initialize counter for screen displays */ 


while(TRUE) { 

frontbuffer(TRUE); 

if (getbutton(MOUSE2) && !(getbutton(MOUSE]1)) && !(getbutton(MOUSE3})) { 
ringbell(); 
while (getbutton(MOUSE2)); 
screencnt += 1; 
if (screencnt == 2) callobj(obj/[SCREEN2)); 
else if (screencnt == 3) callobj(obj[SCREENS3)); 


else break; 


if (getbutton(MOUSE1) && (getbutton(MOUSE2)) && (getbutton(MOUSE3))) { 
*active = FALSE; 
goto exit; 


} 


frontbuffer(F ALSE); 


editobj(obj|FLTPATH]); /* erase previous missile path */ 
objreplace(tag|MISSILE)]); 

circf(0.0, 0.0, 0.0); 

move2(0.0, 0.0); 

draw 2(0.0, 0.0); 

objreplace(tag/TGT)); 

circf(0.0, 0.0, 0.0); 

closeobj(); 


editobj(obj[/STATS]); /* erase previous launch statistics * / 
objreplace(tag|HEAD)); 

charstr("""); 

cmov2i(115,60); 

charstr("""); 

objreplace(tag] TARGET); 

charstr("'"); 

cmov2i(0,0); 

charstr(’™), 

closeobj(); 


setcursor(0,RED,unmask); /* set up cursor and mouse */ 
attachcursor(MOUSEX,MOUSEY); 
setvaluator{MOUSEX,384,0,767); 

setvaiuator| MOUSEY .384,0,767); 

curson(); 


launchlock = FALSE; 
targetlock = FALSE; 


callobj(obj/;CONTOUR]); /* load static displays into both buffers */ 
callobj(obj[INSTR}); 

callobj(obj[STATS]); /* included so swapped buffer doesn’t have "hole" */ 
swapbuffers(}; 
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callobj(obj[CONTOUR)); 
callobj(obj[INSTR]); 


while(TRUE) { 
if (getbutton(MOUSE]1) && (getbutton(MOUSE2)) && (getbutton(MOUSE3))) { 
*active = FALSE; 


goto exit; 


} 


xval = getvaluator(MOUSEX); /* read the x and y mouse positions */ 
yval = getvaluator(MOUSEY); 


utm_x = (50000 + (int)(xval * GRID FACTOR)); /* compute grid coordinates */ 
utm_y = (80000 + (int)(yval * GRID FACTOR)); 


sprintf(xtemp,"%4d",utm_x); /* store coordinates in temporary buffer */ 
sprintf(ytemp,"%4d",utm_y); 


/* if LEFT MOUSE selected lock in launch position and update control panel */ 


if (getbutton(MOUSES3) && (!getbutton(MOUSE2)) && (!getbutton(MOUSE1))) { 
ringbell(); 
xlaunch = xval; 
ylaunch = yval; 
launchlock = TRUE; 
*vx = ((float)((xval * FT_10K)/767)); 
*vz = -((float)((yval * FT _10K)/767)); 
*vy = gnd_level(*vx, *vz) + 200.0; 
editobj(obj[STATS]); 
objreplace(tag|LAUNCH)); 
charstr(xtemp); 
cmov2i(170,220); 
charstr(ytemp); 
closeobj(); 

}  — /* end of MOUSES hit */ 


/* As long as LEFT MOUSE not selected, keep on displaying current UTM 
grid coordinates in control panel area. */ 


if (!launchlock) { 
editobj(obj|STATS)]); 
objreplace(tag/LAUNCH)]); 
charstr(xtemp}; 
cmov2i({170,220); 
charstr(ytemp); 
closeobj(); 


} 


/* if RIGHT MOUSE selected lock in target and update control panel. */ 


if (getbutton(MOUSE1) && (!getbutton(MOUSE3)) && (!getbutton(MOUSE2))) { 
ringbell(); 
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xtarget = xval; 

ytarget = yval; 

targetlock = TRUE; 
editobj(obj[STATS]); 
objreplace(tag/TARGET)); 
charstr(xtemp); 
cmov2i(170,140); 
charstr(ytemp); 

closeobj(); 


/* As long as RIGHT MOUSE not selected keep on displaying current UTM 
grid coordinates in control panel area. */ 
if (!targetlock) { 
if (launchlock) { 
xdistance = ((double)(xval - xlaunch)); 
ydistance = ((double)(yval - ylaunch)); 
distance = sqrt((float)(xdistance * xdistance + ydistance * ydistance)); 
distance = distance * GRID FACTOR; 
sprintf(dist,"%5.0f METERS", distance); 
*direction = atan2(ydistance, xdistance); 
if (*direction < 0.0) *direction += TWOPI, 
*compassdir = compass(*direction); 


sprintf{(heading,"%d DEGREES", (int)*compassdir); 


editobj(obj/STATS]); 
objreplace(tag/ TARGET})); 
charstr(xtemp); 
cmov2i(170,140); 
charstr(ytemp); 
objreplace(tag] HEAD)); 
charstr(heading); 
cmov2i(115,60); 
charstr(dist); 

closeobj(); 


} 


/* if launch position and target location have been selected by the 
operator compute the direction of the missile and distance to target. */ 


if (launchlock && targetlock) { 
xdistance = ((double}(xtarget - xiaunch)); 
ydistance = ((double)(ytarget - ylaunch)); 
distance = sqrt((float}((xdistance * 

(ydistance * ydistance))); 

distance = distance * GRID FACTOR; 
sprintf(dist,""%5.0f METERS", distance); 
*direction = atan2(ydistance, xdistance); 
if (*direction < 0.0) *direction += TWOPI; 
*compassdir = compass(*direction); 


xdistance} — 
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sprintf(heading,"%d DEGREES", (int)*compassdir); 
editobj(obj[STATS)}); 

objreplace(tag] HEAD)); 

charstr(heading); 

cmov2i(115,60); 

charstr(dist); 

closeobj(); 


} 


/* add small red and blue circles to contour map to indicate launch 
position and target location. Connect circles to indicate missile 


flight path */ 


if (launchlock) 

if (targetlock) { 
editobj(obj[FLTPATH)}); 
objreplace(tag{/MISSILE]); 
circf((float)(xlaunch) /767.0* 100.0, (float) (ylaunch)/767.0*100.0, 0.6); 
move2((float)(xtarget) /767.0* 100.0, (float)(ytarget) /767.0* 100.0); 
draw2((float)(xlaunch) /767.0* 100.0, (float) (ylaunch) /767.0* 100.0); 
objreplace(tagiTGT)); 
circf((float)(xtarget) /767.0*100.0, (float) (ytarget) /767.0* 100.0, 0.6); 
closeob}(); 


else { 
editobj(obj[FLTPATH)}); 
objreplace(tag|MISSILE}); 
circf( (float) (xlaunch) /767.0* 100.0, (float)(ylaunch) /767.0* 100.0, 0.6); 
move2((float)(xval) /767.0* 100.0, (float) (yval) /767.0* 100.0); 
draw2((float)(xlaunch)/767.0*100.0, (float )(ylaunch) /767.0* 100.0); 
closeobj(); 


} 


/* if MIDDLE MOUSE selected, launch has occurred and control transfers 
back to main portion of FOG-M program displaying out-the-window 3-D 
view of the flight area. */ 


if (getbutton(MOUSE2) && (!getbutton(MOUSE]1)) && (!getbutton(MOUSE3)) 
&& launchlock && targetlock) { 
ring bell(); 
while (getbutton(MOUSE2)); 
break; 
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writemask(SA VEMAP); 
callobj(obj[FLTPATH)); 


writemask(unmask); 
callobj(obj[STATS}); 
swapbuffers(); 
3 
exit: 
cursoff(); 
popmatrix(); 
} 
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RANDNUM 


/* randnum.c - returns a random float between zero and one me 
static long seed = 1234567; 


randseed(newseed) 
long newseed; 


{ 
} 


seed = newseed; 


float randnum() 


{ 


long mult(); 


seed = (mult(seed,31415821) + 1) % 100000000; 
return(seed / 100000000.0); 


long mult(p,q) 
long P,q; 


long pO, pl, q0, ql; 


pl = p / 10000; 

pO = p % 10000; 

ql = q / 10000; . 

qO = q % 10000; 

return((((p0*q1 + p1*q0) % 10000) * 10000 + p0*q0) % 100000000); 


221 


READCONTROLS 


/* reads the values from the operator’s controls (mouse and dials) */ 


#include "gl.h" /* graphics lib defs */ 
#include "fogm.h" /* fogm constants */ 
#include "device.h" /* device definitions */ 


read controls(designate, greyscale, flying, active, speed, direction, 
compassdir, alt, pan, tilt, fovy) 


int *designate, *greyscale, “flying, *active, *fovy; 
float *speed, *compassdir; 

double *direction, *pan, *tilt; 

Coord *alt; 


{ 


extern float randx, randy, randz; 
float randnum(); 
Colorindex colors] 1]; 


/* quit if all three mouse buttons are pushed */ 

if(getbutton(MOUSE1) && getbutton(MOUSE2) && getbutton(MOUSES3)) { 
*flying = FALSE; 
*active = FALSE; 


else { 
if (getbutton(MOUSE3) && !(getbutton(MOUSE2))) { /* Zoom In */ 
*fovy = (*fovy < (80 + DELTAFOVY)) ? 80: *fovy - DELTAFOVY; 
a 


if (getbutton(MOUSE]1) && !(getbutton(MOUSE2))) { /* Zoom Out */ 
*fovy = (*fovy > (550 - DELTAFOVY)) ? 550: *fovy + DELTAFOVY; 


} 
if (getbutton(MOUSE2)) { /* designate/reject target */ 
if (*designate) { /* see if target in sights */ 
/*pushmatrix(); 
pushviewport(); 
pushattributes(); 


viewport(0, 1023, 0, 767); 
ortho2(0.0, 1022 070 0 16g i). 
cmov 2s((Scoord)(768/2), (Scoord)(768/2)); 
readpixels(1,colors); 
if ((colors[0| >= MIN _TGT COLOR) && (colors|0] <= MAX TGT COLOR)) { 
*designate = FALSE; 
ringbell(); 
randx = 30.0 * randnum() - 15.0; 
randy = 10.0 * randnum() - 5.0; 
randz = 10.0 * randnum(); 


while (getbutton( MOUSE2)); 
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/*} 
popattributes(); 
popviewport(); 
popmatrix(); */ 


else {  /* reject currently designated target */ 
ringbell(); 
*designate = TRUE; 
/* re-adjust tilt and pan values appropriately */ ; 


} 


if (*greyscale != getvaluator(DIAL3)) { /* DIAL3 indicates color change */ 
*oreyscale = !*greyscale; 
setvaluator(DIAL3, *greyscale,0,1); 
colorramp(*greyscale,F ALSE); 


} 


*speed = (float)(getvaluator(DIAL2) / SPEEDSENS); /* get desired speed */ 
*alt = (Coord) (getvaluator(DIAL4)); 


*pan = DTOR * (double)(-getvaluator(MOUSEX)) / PANSENS; 
*tilt = DTOR * (double)(getvaluator(MOUSEY)) / TILTSENS: 


*compassdir = (float)getvaluator(DIALO) / DIRSENS; 

/* keep *direction between 0 and 360, update valuator if changed */ 

if (*compassdir >= 360.0) { 
*compassdir -= 360.0; 
setvaluator(DIALO,(int)(*compassdir* DIRSENS), (int)(-360*DIRSENS), 
(int)(720*DIRSENS)); 


if (*compassdir < 0.0) { 
*compassdir += 360.0; 
setvaluator(DIALO,(int)(*compassdir* DIRSENS), (int) (-360*DIRSENS), 
(int)(720*DIRSENS)); 


/*convert *direction from compass degrees to trigonometric radians */ 
*direction = (*compassdir <= 90.0) ? DTOR * (90.0 - *compassdir) : 
DTOR * (450.0 - *compassdir); 
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READDATA 


/* reads the raw 16 bit elevation and vegetation code data 
from the DMA data file and inserts it into the global 
gridpixel array */ 


#include "fogm.h" 
#include "files.h" 


readdata() 
{ 
int fd; /* file descriptor for the data file */ 
short row, col, rowoffset, coloffset; /* loop indicies */ 
extern short gridpixel{100]/ 100); /* DMA elev and veg. data */ 


/* read the data from the data file into the gridpixel array */ 
fd = open(TERRAIN FILE,RD); 
lseek(fd,0,0); 
for (coloffset = 0; coloffset < NUMXGRIDS * 10; coloffset += 10) { 
for (rowoffset = 0; rowoffset < NUMZGRIDS*10; rowoffset += 10) { 
for (col = 0; col < 10; ++col) { 
for (row = 0; row < 10; ++row) { 
read (fd,&gridpixel|rowoffset--row]||coloffset+col],2); 
} 
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ROAD BOUNDS 


#include "math.h" 
#include "fogm.h" 


#define X 0 
#define Y 1 
#define Z 2 


#define NONE 0 


road _bounds(pt1, pt2, pt3, road width, left ptl, right pt1, left _pt2, 
right pt2, first xgrid, first zgrid, last _xgrid, last_zgrid) 


float pt1(3], pt2[3], pt3[3], road_width; 

float left pt1[3], right pt1[3], left -pt2/3}, right pt2([3}; 

int *first xgrid, *last_xgrid, *first zgrid, *last zgrid; 

{ 
float delta x, delta z, seg dir, min_x, max x, min z, max 2; 
float left “end1(3}, right end1(3}, left "start2[3], right "start2/3}, 
left end2(3], right end2(3}; 
int intersection type; 


/* determine the corner points of the segment */ 

delta_x = pt2[X] - pt1[X]; 

delta z = pt2[Z] - pt1[Z]; 

seg dir = atan2(delta z, delta x); 

left _end1(X] = = pt2[X] + (cos(seg_dir + HALFPI)*road_ width/2.0); 
right_end1{X] = pt2[X] + (cos(seg dir - HALFPI)*road _width/2.0); 
left end1(Z| = pt2{Z] + (sin(seg dir + HALFPI)*road _width/2.0); 

right end1(Z] = pt2[Z] + (sin(seg dir - HALFPI)*road_width/2.0); 


if ((pt2[X] != pt3[X]) |] (pt2[Z] != pt3[Z})) { 
* we are not working with the final segment of this road, find 
the intersection of this segment with the next one */ 
delta _x = pt3[X] - pt2[X]; 
delta _z = pt3(Z] - pt2(Z]; 
seg dir = = atan2(delta z, delta x); 
left_start2/X] = pt2[X] + (cos(seg dir + HALFPI)*road_width/2.0); 
right start2|X] = pea + (cos(seg_dir- HALFPI)*road_width/2.0); 
left start2/Z| = pt2'Z] + (sin(seg dir + HALFPI)*road _width/2.0); 
right _start2|Z| = pt2!Z| + (sin(seg _dir - HALFPI)*road __ width/ 2.0); 
left _end2iX = pt3iX! + (cos(seg dir + HALFP!)*road _width/2 10) 
right _end2(X] = = pt3(X} + (cos(seg_ dir - HALFPI)*road _width/2.0); 
left_ end2[Z] = pt3[Z] + (sin(seg dir + HALFPI)*road _width/2.0); 
right end2[Z] = pt3[Z] + (sin(seg_ dir - HALFPI)*road_width/2.0); 





/* find the intersection point of the left hand sides of the 

first and second road segments */ 
line_intersect2(left pti, left endl, left start2, left end2, 
left pt2, &intersection _type); 


225 


if (intersection type == NONE) { 
left _pt2(X] = left_end1[X]; 
left pt2(Z] = left_end1{Z]; 


/* find the intersection point of the right hand sides of the 
first and second road segments */ 

line intersect2(right ptl, right endl, right start2, right end2, 

right pt2, &intersection type); 7 

if (intersection type == NONE) { 
right pt2(X] = right end1[X]; 
right pt2[Z] = right end1{Z]; 


else { 
/* this is the final segment of this road */ 
left_pt2{X] = left end1[X]; 
left_ pt2{Z] = left_end1{[Z]; 
right pt2/X] = right end1[X]; 
right pt2/Z| = right end1{Z]; 


} 


/* determine the min and max x and z values */ 
min x = left pt1[X]; 
max x = left _pt1[X]; 
min z = left_pt1{Z]; 
max z= left _pt1[Z]; 
if (right _pt1[X] < min_x) min_x = right_pt1[X]; 
if (right pt1[{X] > max x) max x = right pt1|X]; 
if (right pt1[Z] < min_z) min_z= right pt1[Z]; 
if (right_pt1[Z] > max_z) max _z = right _pt1[Z]; 
if (left pt2/X] < min_x) min_x = left _pt2[X]; 

if (left pt2{X] > max_x) max x = left_pt2[X]; 
if (left_pt2[Z] < min_z) min _z = left_pt2[Z]; 
if (left pt2[Z] > max_z) max _z = left_pt2[Z]; 
if (right pt2[X] < min_x) min_x = right_pt2[X]; 
if (right pt2|X] > max_x) max x = right pt2/X]; 
if (right pt2|Z] < min_z) min_z =right pt2([Z]; 
if (right pt2[Z] > max_z) max z= right pt2([Z]; 
*first_ xgrid = (int)(min _x/FT_ _ 100M); 
*first zgrid = (int)(min_z/FT 100M); 
*last xgrid = (int)(max _x/FT_ _ 100M); 
“last zgrid = (int)(max 2; FT 100M); 
if (“firse xgrid < O) “first xgrid = 0; 
if (“first zgrid < 0) “first zgrid = 0; 
if (*last_xgrid > 98) *last _xgrid = 98; 
if (*last zgrid > 98) *last zgrid = 98; 


SORT ARRAY 


sort array(array, num _ entries, decending, test_index) 
float array{10][3}; 
int num _ entries, decending, test index; 
int 1,J; 
* float temp/3}; 


for (i = 0; i < num entries; ++i) { 
for (j = i+ 1;j <= num entries; ++)) { 


if (((decending) && (array|j|[test index] > array|i|[test index])) || 
((!decending) && (array{j|[test_ index] < array|i|{test index]))) { 
temp(0] = array(i]/0]; 
temp(1] = array/i]/1]; 
temp|2] = array(i|/2]; 
array|i/[0] = array/}/(0); 
array|i|(1| = array|j](1]; 
array|i|(2] = array(j]|2]; 
array(j|(0]| = temp|0}; 
array|j}{1] = temp/1]; 
array (j|/2! = temp/2]; 
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UP LOOK POS 


/* compute the camera’s lookat position */ 


#include "fogm.h" /* fogm constants */ 
#include 'math.h" /* math routine definitions */ 
#include "gl.h" /* graphics definitions */ 


update look posit(direction, pan, tilt, vx, vy, vz, 
tgtx, tgty, tgtz, designate, px, py, pz) 


double direction, pan, tilt; 
Coord vx, vy, vz, tgtx, tgty, tgtz, *px, *py, *pz; 
int designate; 


{ 


extern int framecnt; 


double lookdir; 
if (designate) { /* missile is not locked on to a target */ 


/* compute direction camera is looking */ 
lookdir = direction + pan; 


/* compute a coordinate along camera’s line of sight */ 
*px = vx + cos(lookdir) * MAXLOOKDIST; 
*pz = vz - sin(lookdir) * MAXLOOKDIST; 


if (framecnt < 15) { 
*py = 4.0 * vy * (14 - framecnt) / 14.0; 
framecnt+-+; 


else { 
*py = vy + MAXLOOKDIST * tan(tilt); 


else { 
"Dx = teuK: 
“py = tety; 
*pz = teen 
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UP MSL POSIT 


/* Compute new missile position */ 


#include "gl.h" /* graphics definitions */ 

#include "device.h" /* graphics device definitions */ 
#include "fogm.h" /* fogm constants */ 

#include "math.h" /* math function declarations */ 


#include <sys/types.h> /* contains the time sturcture tms */ 
#include <sys/times.h> /* for time calls */ 


update missile posit(direction, compassdir, speed, designate, 
tgtx, tgty, tgtz, vx, vy, vz, flying) 


double *direction; 
float *compassdir; 
float speed; 

int designate; 

Coord tgtx, tgty, tgtz; 
int *flying; 

Coord *vx, *vy, *vz; 


{ 


static long seconds; 

static long lastsec = -999; /* -999 is flag to indicate no value */ 
struct tms timestruct; 

float deltadist, gndlevel, gnd_level(), compass(), ht above tank; 
long float deltax, deltaz, dist to tank; 


seconds = times(&timestruct); 


/* compute distance missile must move ahead to maintain speed */ 
if (lastsec == -999) 

deltadist = 0.0; 
else 


deltadist = (speed/FPS TO KTS) *(seconds - lastsec); 


lastsec = seconds; /* save for next pass */ 

if (designate) { /* missile under operator contol, not locked on tgt */ 
*vx += deltadist * cos(*direction); 
*vz -= deltadist * sin(*direction); 


/* keep missile at least 50 ft above ground level ~ 
gndlevel = gnd _level(*“vx, *vz); 
if (*vy < (gndlevel + 50.0)) *vy = gndlevel + 50.0; 


} 
else { 


deltax = *vx - tgtx; 
deltaz = *vz - tgtz; 


dist to tank = hypot(deltax, deltaz); 
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if (deltadist > (float)dist to tank) { /* hit on target */ 
deltadist = (float)dist to tank - 5.0; 
*flying = FALSE; 
lastsec = -999; /* no value flag for next launch */ 


} 


*direction = (double)atan2((float)deltaz, (float)-deltax); 
if (*direction < 0.0) *direction += TWOPI; 
*compassdir = compass(*direction); 


setvaluator(DIALO, (int) (*compassdir* DIRSENS), (int)(-360* DIRSENS), 
(int)(720* DIRSENS)); 


*vx += (deltadist * cos(*direction)); 

*vz -= (deltadist * sin(*direction)); 

ht above tank = (float)*vy - gnd_level(tgtx,tgtz); 

*vy -= (Coord)((ht above tank * deltadist) / (float)dist to tank); 
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VIEW BOUNDS 


#include "fogm.h" 
#include "gl.h" 
#include "math.h" 


view bounds(vx, vy, vz, px, py, pz, tilt, fovy, 
firstxerid, firstzgrid, lastxgrid, lastzgrid) 
Coord vx,vy,vZ; 

double tilt; 


int fovy; 
short *firstxgrid, *firstzgrid, *lastxgrid, *lastzgrid; 
float ix, iz; | /* the intersection points */ 


float lookdir; 

float deltax, deltay, deltaz, delta alt, fx, fy, fz; 
float half fovy; 

float lower edge angle; 


/* compute the direction the camera is looking */ 
lookdir = atan2((float)(vz - pz), (float) (-(vx-px))); 
if (lookdir < 0.0) lookdir += TWOPIT; 


if (vy > py) { 
/* tilt angle is negative */ 
deltax = px - vx; 
deltay = py - vy; 
deltaz = pz - vz; 


delta alt = pow((float)MIN, ALTSCALE) - vy; 


else { 
/* tilt angle is positive, use the lower fustrum edge instead 
of the line of sight to compute the view bounds */ 
/* compute a coordinate along the lower fustrum edge */ 
half fovy = ((float)fovy /20.0*DTOR); 
lower edge angle = tilt - half fovy; 
fx = vx + cos{lookdir)*MAXLOOKDIST; 
fz = vz - sin(lookdir)*MAXLOOKDIST; 
fy = vy + tan(lower edge angle)*MAXLOOKDIST; 
deltax = fx - vx; 
deltay = fy - vy; 
deltaz — 17 = vz; 
delta alt = pow((float}MIN, ALTSCALE) - vy; 
} 


ix = vx + ((deltax/deltay)*delta alt); 
iz = vz + ((deltaz/deltay)*delta alt); 


/* compute which grid objects should be sent through the geometry 
pipeline */ 
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if (deltay > 0.0) { 
/* the fustrum is looking totally skyward, don’t bother doing 
any terrain */ 
*firstxgrid = 0; 
*firstzgrid = 0; 
*lastxgrid = 0; 
*lastzgrid = 0; 


else { 
/* display 20 grid squares on all sides of the intersection point */ 
*firstxgrid = (int)(ix/FT_100M) - 20; 
*lastxgrid = (int)(ix/FT 100M) + 20; 
*firstzgrid = (int)(-iz/FT 100M) - 20; 
*lastzgrid = (int)(-iz/FT 100M) + 20; 


/* insure that objects drawn include the current missile position */ 

if ((int)(vx/FT 100M) < “*firstxgrid) 
*firstxgrid = (int)(vx/FT 100M); 

if ((int)(vx/FT 100M) > *lastxgrid) 
*lastxgrid = (int)(vx/FT_ 100M); 

if ((int)(-vz/FT 100M) < *firstzgrid) 
*firstzgrid = (int)(-vz/FT_ 100M); 

if ((int)(-vz/FT 100M) > *lastzgrid) 
*lastzgrid = (int)(-vz/FT 100M); 

if (*firstzgrid < 0) *firstzgrid = 0; 

if (*firstxgrid < 0) *firstxgrid = 0; 

if (*lastzgrid > 98) *lastzgrid = 98; 

if (*lastxgrid > 98) *lastxgrid = 98; 
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